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PREFACE 


The  Atmospheric  Physics  Branch  of^tW  Ocean  Sciences  Division, of  the 
Naval  Research  Laboratory  organized  andcexecuted  a multi  disciplined  cruise 
of  the  USNS  Hayes  during  May  and  June  1977.  The  cruise  was  'devoted  to  ob- 
serving the  electro-optical  and  meteorological  (EOMET)  properties  of  ma- 
rine atmosphere  using  state  of  the  art  instrumentation.  , 

The  aim  of  this  reportv^s^to  organize  in  one  place  the  results  of 
this  cruise.  The  data  of  most  of  theN$nvest1 gators  were  sent  to  the  edi- 
tor to  be  incorporated  into  this  report  and  thus~it^!encompasses  a fairly 
complete  description  of  the  results  of  this  cruise.  In  an  effort  to  or- 
ganize the  data  archive  this  data  was  transferred  by  various  means  to  the 
disc  files  of  the  NRL  PDP-10  computer.  Although  considerable  care  was 
exercised  in  this  process,  there  is  no  doubt  that  in  an  effort  of  thjs 
magnitude  that  errors  were  made  in  this  transfer  process.  Consequently, 
the  credit  for  the  measurements  belong  to  the  individual  investigators 
but  the  blame  for  any  erroneous  data  must  be  on  the  editor's  shoulders. 


This  report  is  organized  with  the  view  of  presenting  the  data  to  be 
of  particular  use  to  the  atmospheric  modeler  who  wishes  to  relate  electro- 
optical  transmission  parameters  produced  by  atmospheric  aerosols  to  the 
meteorological  parameters.  Consequently,  after  the  introduction  in  Sec- 
tion A,  the  data  portion  of  this  report.  Section  B,  is  divided  into  four 
subsections.  The  first  describes  the  instrumentations  and  locations  of 
the  instrumentations  reported  in  a continuous  manner  in  this  report  and 
its  data.  The  second  describes  the  geographical  aspects  of  the  cruise. 
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The  third  part  Is  devoted  to  the  key  meteorological  parameters  while  the 
fourth  part  of  the  section  is  devoted  to  describing  the  various  aspects 
of  the  aerosol  measurements. 

Section  C of  this  report  contains  various  signed  papers  describing 
special  experiments  which  were  carried  out  throughout  the  cruise. 

Section  0 describes  the  data  archive  system  and  what  products  are 
available  from  this  system. 

Section  E is  a list  of  participants. 
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A.  INTRODUCTION 


This  report  Is  the  compilation  and  condensation  of  the  data  obtained 
on  a 21  day  research  cruise  aboard  the  USNS  Hayes  which  was  designed  to 
investigate  the  relationships  which  are  thought  to  exist  between  electro- 
optical  phenomenon  and  marine  meteorology. 

Atmospheric  boundary  layer  research  on  the  high  seas  is  always  dif- 
ficult to  carry  out.  There  are  relatively  few  alternative  platforms  which 
are  available  to  investigators  who  wish  to  make  measurements  on  the  marine 
atmosphere  in  order  to  have  real  data  on  which  to  base  their  theoretical 
deductions  and  with  which  to  test  the  results  of  their  mathematical  models. 

Platforms  which  have  been  used  to  obtain  representative  data  samples 
can  be  roughly  placed  in  the  following  categories: 

a)  aircraft 

b)  ships 

c)  buoys 

d)  islands 

e)  shore  based  instrumentation  sites 

f)  offshore  man-made  platforms 

g)  satellites 

Both  in  situ  and  remote  sensing  of  the  marine  boundary  layer  can  be 
carried  out  by  the  first  six  of  the  vehicle  types  whereas  satellite  mea- 
surements are  by  necessity  strictly  done  by  remote  sensing. 

Note:  Manuscript  submitted  December  14,  1978. 
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There  are  pros  and  cons  to  all  of  these  platform  types.  Instrumen- 
tation of  fixed  locations  such  as  islands,  buoys,  or  platforms  offer  the 
opportunity  for  long  term  measurements  but  are  very  inflexible  in  obtain- 
ing geographical  variations  or  physical  variations  based  on  geographical 
location.  On  the  other  hand,  measurements  made  aboard  moving  platforms 
such  as  ships,  aircraft  offer  great  flexibility  in  geographical  location 
but  the  length  of  endurance  of  measurements  is  limited  by  economic  factors. 

The  cruise  of  the  USNS  Hayes  was  designed  to  skirt  the  east  coast  of 
the  continental  United  States  to  obtain  representative  aerosol  measure- 
ments and  to  monitor  their  evolution  in  time  as  the  air  mass  proceeds  to 
change  from  continental  to  oceanic. 

The  cruise  obtained  EOMET  data  with  extremes  in  stability  not  often 
encountered  in  the  open  ocean  but  which  is  very  much  present  off  shore  in 
the  areas  of  Nova  Scotia  and  Newfoundland  where  cold  coastal  waters  and 
warm  continental  air  mass  combine  to  produce  very  stable  boundary  layers. 

A hundred  kilometers  south  and  east  of  these  stable  areas,  the  warm  Gulf 
Stream  produces  a high  surface  instability  on  the  air  masses  which  have 
started  to  be  accommodated  to  cold  areas  closer  to  the  shore. 

Finally,  the  cruise  offers  a geographical  comparison  in  which  the 
same  instrumentation  was  used  to  look  for  differences  between  mid-Atlantic 
Ocean  and  Mediterranean  Sea  situations. 

The  overall  mission  of  this  cruise  was  to  obtain  from  a wide  data 
base  various  simultaneous  measurements  of  electro-optics  parameters  made 
simultaneously  with  marine  meteorology  measurements  in  order  to  provide 


-2- 


both  a data  base  on  which  to  base  the  construction  of  needed  EOMET  models 
and  a test  set  of  data  on  which  to  test  existing  models.  This  latter  use 
of  the  data  is  such  that  both  a measured  EOMET  parameter  and  the  simulta- 
neously measured  input  meteorological  parameters  are  recorded  for  the  same 
period  of  time.  The  EOMET  model  prediction  can  then  be  calculated  from 
the  measured  meteorological  input  parameters  and  compared  with  the  mea- 
sured predicted  parameter. 

Several  similar  sets  of  data  made  simultaneously  by  different  inves- 
tigators using  different  techniques  were  routinely  taken  of  many  critical 
parameters  throughout  the  cruise.  This  apparent  redundancy  serves  several 
valuable  purposes. 

1)  It  reduces  the  number  of  times  that  a particular  parameter  is  not 
being  measured  because  of  instrumentation  malfunction  (a  not  too 
rare  phenomenon  for  shipboard  instrumentation  in  a hostile  salt 
environment) . 

2)  The  observation  of  several  independent  instrumentation  systems 
reading  nearly  the  same  values  simultaneously  adds  credance  to 
the  validity  of  the  measurement. 

3)  It  provides  a set  of  data  for  the  detection  of  systematic  errors 
in  particular  instrumentation  systems  or  instrumentation 
locations. 

As  with  all  experiments  which  have  been  undertaken  with  many  investi- 
gators and  with  automatic  data  processing  machines,  the  total  amount  of 
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data  collected  and  stored  during  this  research  cruise  is  very  large.  The 
printing  of  reams  of  paper  filled  with  many  columns  of  numbers,  each  num- 
ber representing  the  data  of  a particular  instrument  at  a particular  sec- 
ond  of  the  cruise  would  be  completely  out  of  the  question.  Not  only  would 
the  cost  of  publishing  be  prohibitive  but  the  results  would  be  difficult 
for  any  potential  data  user  to  handle  because  of  the  great  bulk  of  printed 
matter.  The  theoretician  who  it  is  hoped  will  utilize  this  data  in  both 
the  construction  and  testing  of  his  models  is  in  danger  of  losing  sight 
of  the  proverbial  forest  of  meaning  because  of  the  trees  of  numbers  in 
his  way. 

There  have  been  attempts  at  obtaining  a solution  to  this  modern  day 
dilemma  of  seemingly  too  much  data  but  in  our  opinion  none  of  these  solu- 
tions is  entirely  satisfactory.  Our  solution  to  this  problem  is  based  on 
the  access  of  interested  investigators  to  the  data  on  three  levels. 

Our  philosophy  for  the  presentation  of  the  data  for  the  scientific 
user  is  as  follows.  The  first  level  is  a general  overview  of  the  best  es- 
timates of  the  various  physical  parameters  measured  throughout  the  cruise 
which  are  present  in  terms  of  graphical  displays  found  in  Section  B. 

These  curves  will  allow  individual  investigators  to  select  portions  of  the 
cruise  of  interest  to  them.  The  graphical  presentations  however  do  not 
give  the  precision  which  exists  in  the  original  data  either  in  the  values 
themselves  or  in  the  time  resolution.  Therefore  a second  level  of  data 
presentation  or  access  is  necessary  to  enable  the  data  to  be  utilized  to 
the  extent  envisioned  in  the  original  plans  of  the  cruise.  Listings  of 
key  parameters  are  given  in  the  tables  of  Section  B.  In  addition  all  of 
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the  cruise  data  is  available  on  a second  level  of  data  access  which  is  de- 
signed to  give  the  investigator  individualized  access  to  any  portions  of 
the  data  which  he  might  need.  This  individualized  approach  is  described 
in  Section  0 of  this  report.  The  data  in  the  data  archive  of  the  cruise 
contain  hourly  values  recorded  by  each  of  the  instruments  when  they  were 
in  apparently  working  operation.  In  addition  there  are  a number  of  files 
which  contain  the  calculations  of  certain  parameters  which  are  based  on 
measured  data.  Software  exists  for  putting  this  data  together  in  any  de- 
sired order  so  that  customized  data  sorting  is  easily  achieved. 

Finally  most  of  the  original  investigators  have  copies  of  the  raw 
data  which  they  will  probably  make  available  to  others  on  an  individua1 
basis. 


Data  files  are  referred  to  often  throughout  this  report.  The  mean- 
ings of  these  names  and  how  to  find  them  are  described  in  Section  D of 
this  report. 
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B.  CRUISE  DATA:  GENERAL  RESULTS 

I.  Cruise  Map: 

Figure  B-I-I  shows  the  ships  track  of  the  USNS  Hayes  while  carry- 
ing out  the  EOMET  77  experiments.  The  cruise  was  divided  into  two  non 
equal  segments  separated  by  a 2 day  period  in  port  at  Rota,  Spain.  The 
North  Atlantic  leg  was  a slight  modification  to  a great  circle  route.  The 
modification  was  a short  deviation  into  Nova  Scotian  waters  in  an  endeavor 
to  repeat  certain  of  the  fog  studies  that  were  obtained  during  the  1975 
fog  cruise  of  the  USNS  Hayes.  Unfortunately  conditions  were  such  cnai  no 
appreciable  fogs  were  encountered  during  the  1977  EOMET  cruise. 

The  cruise  was  planned  to  provide  a wide  variety  of  marine  atmospher- 
ic conditions.  The  ships  track  lead  up  the  eastern  coast  of  the  North 
American  continent  where  a strongly  continental  type  of  aerosol  was  en- 
countered. Very  stable  boundary  layer  conditions  were  encountered  in  the 
waters  off  of  Nova  Scotia  and  Newfoundland  during  May  18,  19,  and  20,  1977. 
Very  unstable  conditions  were  experienced  as  the  ship  passed  rapidly  from 
the  cold  waters  of  the  Grand  Banks  area  into  the  warm  waters  of  the  Gulf 
Stream  on  23  May  1977. 

One  unfortunate  feature  of  the  track  was  that  the  strongly  eastward 
component  of  the  ships  direction  paralleled  the  main  westerly  winds  of 
this  area.  Thus  during  periods  of  strong  real  wind  conditions,  tne  rela- 
tive wind  direction  experienced  at  the  ship  was  from  time  to  time  from  the 
stern  and  thus  the  aerosol  sampling  apparatus  located  for  optimum  sampling 
for  normal  head  winds  would  give  erroneous  results  unless  the  ship  was 
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turned  into  the  wind  periodically  to  produce  the  desired  sampling  orienta- 
tion of  the  ship  heading  into  the  wind. 

The  remainder  of  the  cruise  gave  samples  of  mid-Atlantic  Ocean  condi- 
tions, coastal  conditions  off  of  the  Iberian  Peninsula  and  conditions  in 
the  Mediterranean  Sea. 

Table  B-I-l  relates  Julian  dates  used  throughout  this  report  with 
calendar  month  and  day. 

TABLE  B-I-l 


Julian  Date  vs.  Calendar  Month,  Day 


Julian 

Month 

Date 

Day 

134 

May  14 

135 

May  15 

136 

May  16 

137 

May  17 

138 

May  18 

139 

May  19 

140 

May  20 

141 

May  21 

142 

May  22 

143 

May  23 

144 

May  24 

145 

May  25 

146 

May  26 

147 

May  27 

148 

May  28 

149 

May  29 

150 

May  30 

151 

May  31 

152 

June  1 

153 

June  2 

154 

June  3 

155 

June  4 

156 

June  5 

157 

June  6 

158 

June  7 

II.  Instrumentation  Description 


a)  NRL  Instrumentation 

USNS  Hayes  Instrumentation  for  EOMET  Cruise— General  Comments 

The  USNS  Hayes  (T-ACOR  16)  is  a large  (3180  tons)  oceanographic 
research  ship  configured  in  the  form  of  a catamaran  246  feet  (75  meters) 
in  length.  A detailed  description  of  the  ship  and  its  facilities  is  given 
by  A.  T.  McCl inton  in  NRL  Report  #7370.  In  addition  to  its  11  officers 
and  34  crew  members,  the  EOMET  77  cruise  berthed  31  scientists  and  tech- 
nicians in  two  legs  of  the  cruise.  Dr.  L.  H.  Ruhnke  of  NRL  was  Senior 
Scientist  on  board  (SSOB).  Stuart  Gathman,  NRL,  was  Assistant  SSOB. 

A list  of  participants  and  their  affiliations  is  in  Section  E. 

Figure  B-II-1  shows  in  a schematic  manner,  the  general  location  of 
the  instrumentation  installed  on  the  ship  for  the  EOMET  77  cruise. 

Care  was  taken  in  placing  the  various  sensors  which  are  sensitive  to 
sampling  techniques  to  the  most  favorable  locations.  These  locations  were 
determined  by  scale  model  tests  of  the  USNS  Hayes  in  a wind  tunnel.  The 
least  disturbed  sampling  areas  for  head  winds  on  the  USNS  Hayes  are  loca- 
tions on  the  forward  mast,  on  towers  located  on  the  bows  of  the  ship,  and 

1 1 

on  towers  mounted  on  top  of  the  pilot  house.  These  locations  are  shown 
In  Figure  B-II-1,  as  locations  A,  C and  D,  and  B,  respectively. 

Because  of  the  potential  problem  of  aerosol  diffusion  to  sampling 
tube  walls,  the  various  aerosol  analyzers  were  located  either  in  situ  or 
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on  the  areas  described  above  or  just  below  deck  in  the  forward  part  of 
the  ship,  location  L,  with  relatively  short  sampling  tubes  going  to  input 
ports  at  location  E. 

Sea  surface  temperatures  were  obtained  from  probes  pulled  through 
the  surface  waters  between  the  hulls  at  locations  R i S,  from  probes  lo- 
cated in  a well  at  the  bottom  of  the  hull  at  location  G and  from  infrared 
thermometers  located  at  point  B but  focused  at  *.ne  undisturbed  sea  surface 
near  the  ship.  The  remote  sensing  systems  were  located  in  the  center  and 
towards  the  stern  of  the  ship  at  locations  I,  H,  and  K. 

Finally  the  kite  balloon  sounding  system,  J,  was  operated  from  point 
T where  its  hangar  and  winch  were  located.  This  device  could  be  operated 
only  when  the  relative  wind  vector  was  within  certain  limits  in  both  speed 
and  direction.  Location  M was  main  recorder  room;  P was  computer  room. 

The  NRL  Meteorological  Instrumentation  Descriptions 
This  section  of  the  report  will  describe  the  instrumentation  used 
to  obtain  the  NRL  Code  8320  meteorological  data  contained  in  the  data 
archive. 


1)  A Y.S.l.  Scanning  Tele-Thermometer  model  47  was  used  to  moni- 
tor several  temperatures  at  locations  throughout  the  ship.  Temperatures 
are  measured  to  1/10°C  with  this  device.  This  device  was  developed  to 
sample  sequentially  seven  matched  thermometers  every  hour.  Included  in 
the  7 s«mples  were  two  fixed  resistors  used  as  calibration  monitor  points 
for  the  device.  Two  of  the  sensors  were  mounted  on  a carriage  inside  the 
starboard  midship  instrumentation  well  of  the  ship  and  lowered  to  keel 
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level,  approximately  2 centimeters  from  the  water  flow  lines.  Their  loca- 
tion in  Figure  B-II-1  is  point  G.  These  devices  thus  sample  sea  water 
temperature  at  a depth  of  approximately  seven  meters.  This  data  is  found 
in  data  file  THERM.DAT  in  column  2. 

2)  One  thermometer  was  located  inside  of  a ventilated  sun  shield 
on  a tower  at  location  B in  Figure  B-II-1.  Its  data  thus  represents  the 
air  temperature  at  an  altitude  of  19  meters  above  sea  level.  This  data  is 
found  in  data  file  THERM.DAT,  column  3. 

3)  A Cambridge  Systems,  model  110S-M  automatic  meteorological 
temperature  and  dew  point  measuring  system  was  installed  on  a tower  at 
point  B at  an  altitude  of  approximately  19  meters  above  sea  level.  Both 
temperatures  and  dewpoints  with  this  device  are  accurate  to  0.1°C.  Tem- 
perature data  is  stored  in  FILE.DAT,  column  2 and  dewpoint  data  is  in 
FILE.DAT,  column  3. 

4)  Relative  wind  speed  and  directions  were  obtained  with  a Ben- 
dix  type  F/60  windsoeed  indicator  mounted  on  the  forward  mast  (radar  mast) 
giving  wind  speed  to  an  accuracy  of  plus  or  minus  1 knot  and  a direction 
accuracy  of  + 2°  at  wind  speeds  above  5 knots.  This  data  is  found  in  file 
ASC.DAT.  Column  2 of  this  file  contains  the  wind  speed  in  knots  and  col- 

1 1 

umn  3 is  the  relative  direction  in  degrees. 

5)  Atmospheric  radon  concentration  was  obtained  from  the  NRL  de- 
signed Automatic  Radon  Counter  (R.  E.  Larson,  B.  J.  Bressan,  "ARC  for  Con- 
tinental Unattended  Operations"  R.S.I.,  49  (7),  p 965-969,  July).  Hourly 
data  from  this  device  is  found  in  RAD0N2.DAT,  column  2. 
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6)  Turbidity  measurements  were  made  with  the  Voltz  sun  photom- 
eter at  wavelengths  of  0.88p,  0.5u,  and  0.94m.  Values  of  the  air  mass 
were  calculated  at  each  observational  point  from  knowledge  of  the  ship's 
position  and  time  using  various  programs  in  the  Hewlett-Packard-65  naviga- 
tional package. 


Air  mass  data  is  in  file  V0LTZ.DAT,  column  2.  The  air  mass  has  a 
minimum  value  of  1.0  for  situations  in  which  the  sun  is  at  the  zenith. 
Turbidity  measurements  made  with  calibrated  Voltz  sun  photometers  were 
made  throughout  the  cruise  at  three  different  wave  lengths.  Calculations 
of  the  optical  thickness  due  to  aerosols  were  made  using: 


Taerosol  m *n  T 


0.157 


where  the  air  mass  is  m,  I the  calibration  of  the  device,  and  I the  mea- 

o 

sured  intensity.  The  other  effects,  i.e.,  ozone  and  Rayleigh  scattering, 
are  assumed  to  be  constant  and  are  represented  by  a thickness  of  0.157. 


File:  Turbid. DAT ( 2 ) is  the  calculated  thickness  at  x = 0.88u 

File:  Turbid. DAT ( 3 ) is  the  calculated  thickness  at  x = 0.5m 

and  File:  Turbid. DAT(4)  is  the  calculated  thickness  at  x = 0.94m 


7)  The  ship's  positions  were  constantly  being  pinpointed  by  a 
professional  scientific  navigator  Mr.  Tony  Zuccaro  of  NRL  who  relied  in 
part  on  satellite  observations  using  the  Magnavox  706C  systems.  A data 
file  NAV.DAT  relates  the  ship's  positions  for  each  hour  of  the  cruise. 


Column  2 contains  the  latitude  in  whole  degrees. 
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Column  3 contains  the  additional  fraction  of  degrees  expressed  in 


decimal  minutes 

Column  4 contains  the  longitude  in  whole  degrees 

Column  5 contains  the  additional  fraction  of  degrees  expressed  in 
decimal  minutes 

8)  Ship's  speed  and  heading  were  calculated  from  the  scientific 
navigational  log  which  was  oased  on  satellite  observations  using  the  Mag- 
navox  706C  system.  Ship's  course  made  good  data  is  stored  in  NAV . DAT ( 6 ) 
and  the  ship's  speed  data  is  stored  in  NAV . DAT ( 7 ) . 

9)  Remote  sensing  of  the  skin  temperature  of  the  sea  surface 
was  done  by  a Barnes  Engineering  Co.  Precision  Radiation  Thermometer 
model  PRT-5.  This  data  is  found  in  FILE.DAT,  column  4. 

10)  Hourly  radiational  observations  of  the  marine  boundary  were 
obtained  by  means  of  the  Barnes  Engineering  Company's  Field  Radiation 
Thermometer,  model  PRT-10.  Files  which  contain  this  data  are: 

Sea  Surface  Temperature--PRT10.DAT(2) 

Sky  Temperature--PRT10.DAT(3) 

PRT-10  aimed  10  degrees  below  the  horizon--PRT10.DAT(4) 

PRT-10  aimed  10  degrees  above  the  horizon--PRT10.DAT(5) 

PRT-10  aimed  30  degrees  above  the  horizon--PRT10.DAT(6) 

PRT-10  aimed  60  degrees  above  the  horizon--PRT10.DAT(7) 

PRT-10  aimed  at  the  cloud  base--PRT10.DAT(8) 


11)  There  are  several  subjective  types  of  observations  of  the 
general  weather  conditions  which  have  been  coded  to  numerical  values  to 
depict  various  important  atmospheric  characteristics.  These  codes  are 
all  found  in  the  file  CL0UD.DAT.  Column  2 of  this  file  contains  the 
coded  seastate  conditions  at  the  time  of  observation.  The  numerical 
code  of  this  data  is  given  in  Table  B-II-1. 

Column  3 contains  estimates  of  foam  and  whitecap  cover  as  a percent- 
age of  the  sea  surface.  Column  4 includes  approximate  visibility  in 
miles.  In  instances  where  no  haze  was  distinguishable  at  the  horizon, 
making  subjective  estimates  difficult,  this  clarity  is  assigned  the 
value  of  100. 

Column  5 corresponds  to  the  percentage  to  total  cloud  cover.  Again, 
this  number  represents  a visual  approximation.  The  following  columns 
(6-8)  break  down  the  total  cloud  cover  into  its  component  types:  low 
(column  6),  middle  (column  7)  and  high  clouds  (column  8).  These  types 
have  been  coded  numerically  based  upon  WMO  standards,  as  seen  in  Tables 
B-II-2,  B-II-3,  and  B-II-4. 

The  CL0UD.DAT  file  is  comprised  entirely  of  eye  observations  made  at 
regular  daylight  intervals.  This  meteorological  information  can  be  use- 
ful in  conveying  a general  picture  of  weather  conditions  as  they  existed 
throughout  the  EOMET  cruise. 
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Table  B-II-1 
SEASTATE  CONDITIONS 


Code 

Figure 

Description  of  sea 

Height  of 
waves  in  feet 

Height  of 
waves  in  meters 

00 

Calm  (glassy)  ----------- 

0 

0 

01 

Calm  (rippled)  ----------- 

0 - 1/3 

0 - 0.1 

02 

Smooth  (wavelets)  --------- 

1/3  - 1-2/3 

0.1  - 0.5 

03 

Slight  --------------- 

1-2/3  - 4 

0.5  - 1.25 

04 

Moderate  -------------- 

4-8 

1.25  - 2.5 

05 

Bough  --------------- 

8-13 

2.5  - 4 

06 

Very  rough  ------------- 

13  - 20 

4 - 6 

07 

High 

20  - 30 

6 - 9 

08 

Very  high  ------------- 

35  - 45 

9 - 14 

09 

Phenomenal  ------------- 

Over  45 

Over  14 

Notes : 

(1)  The  average  wave  height  as  obtained  from  the  larger 
well-formed  waves  of  the  wave  system  being  observed  is 
reported. 

(2)  If  an  exact  boundary  height  could  be  reported  by  two 
code  figures  the  lower  code  figure  will  be  reported;  e.g.f 

a height  of  13  feet  would  be  reported  by  code  figure  5 or  05. 
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Table  B-II-2.  Low  Cloud-WMO  Code 


Symbol  Cl=C1oucIs  of  Genera  Sc,  St,  Cu,  Cb 


Code 

Fig- 

ure 

Technical  Specifica- 
tions 

Nontechnical  Specifi- 
cations 

0 

No  CL  clouds 

No  Stratocumulus, 
Stratus,  Cumulus, 
or  Cumulonimbus. 

1 

Cumulus  humilis  or 

Cumulus  with  little 

Cumulus  fractus 

vertical  extent  and 

other  than  of  bad 

seemingly  flattened, 

weather,  or 

or  ragged  Cumulus 

both. 

other  than  of  bad 
weather,1  or  both. 

2 

Cumulus  mediocris 

Cumulus  of  moderate 

or  congestus, 

or  strong  vertical 

with  or  without 

extent  generally 

Cumulus  of 

with  protuberances 

species  fractus  or 

in  the  form  of 

humilis,  or  Stra- 

domes  or  towers, 

tocumulus,  all 

either  accompanied 

having  their 

or  not  by  other 

bases  at  the 

Cumulus  or  by 

same  level. 

Stratocumulus;  all 
having  their  bases 
at  the  same  level. 

3 

Cumulonimbus 

Cumulonimbus,  the 

calvus,  with  or 

summits  of  which, 

without  Cumulus, 

at  least  partially, 

Stratocumulus  or 

lack  sharp  outlines, 

Stratus. 

but  are  neither 
clearly  fibrous  (cir- 
riform)  nor  in  the 
form  of  an  anvil; 
Cumulus,  Stratocu- 
mulus or  Stratus 
may  also  be  present. 

4 

Stratocumulus  cu- 

Stratocumulus  formed 

mulogenitus. 

by  the  spreading  out 
of  Cumulus;  Cumu- 
lus may  also  be 
present. 

5 

Stratocumulus 

Stratocumulus  not  re- 

other  than  Stra- 

suiting  from  the 

tocumulus  cumu- 

spreading  out  of 

logenitus. 

Cumulus. 

Code 

Fig- 

ure 

Technical  Specifica- 
tions 

Nontechnical  Specifi- 
cations 

6 

Stratus  nebulosus 
or  Stratus  fractus 
other  than  of  bad 
weather,1  or  both. 

Stratus  in  a more  or 
less  continuous  sheet 
or  layer,  or  in  rag- 
ged shreds,  or  both, 
but  no  Stratus 
fractus  of  bad 
weather.1 

7 

Stratus  fractus  or 
Cumulus  fractus 
of  bad  weather,1 
or  both  (pannus), 
usually  below 
Altostratus  or 
Nimbostratus. 

Stratus  fractus  of  bad 
weather  1 or  Cumu- 
lus fractus  of  bad 
weather,1  or  both 
(pannus),  usually 
below  Altostratus  or 
Nimbostratus. 

8 

Cumulus  and  Stra- 
tocumulus other 
than  Stratocumu- 
lus cumulogeni- 
tus,  with  bases 
at  different  levels. 

Cumulus  and  Strato- 
cumulus other  than 
that  formed  from  the 
spreading  out  of 
Cumulus;  the  base 
of  the  Cumulus  is  at 
a different  level 
from  that  of  the 
Stratocumulus. 

9 

Cumulonimbus 
capillatus  (often 
with  an  anvil), 
wi*h  or  without 
Cumulonimbus 
calvus,  Cumulus, 
Stratocumulus, 
Stratus  or 
pannus. 

Cumulonimbus,  the 
upper  part  of  which 
is  clearly  fibrous 
(eirriform),  often  in 
the  form  of  an  anvil, 
either  accompanied 
or  not  by  Cumulo- 
nimbus without  anvil 
or  fibrous  upper 
part,  by  Cumulus, 
Stratocumulus, 

Stratus  or  pannus 
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Table  B-II-3.  Middle  Cloud-WMO  Code 


Symbol  Cm=C1ouc1s  of  Genera  Ac,  As,  Ns 


Code 

Technical  Specifics- 

Nontechnical  Specifi-  Code 

Technical  Specifics 

Nontechnical  Specifi- 

Ffg- 

tions 

cations  Fig- 

tions 

cations 

ure 

ure 

5 

Altocumulus  trans- 

Semitransparent  Alto- 

0 

No  Cm  clouds 

No  Altocumulus,  Alto- 

lucidus  in  bands, 

cumulus  in  bands,  or 

stratus  or  Nimbo- 

or  one  or  more 

Altocumulus  in  one 

stratus. 

layers  of  Alto- 

or  more  fairly  con- 

1 

Altostratus  trans- 

Altostratus,  the  greater 

cumulus  trans 

tinuous  layers  (semi- 

lucidus. 

part  of  which  is  semi- 

lucidus  or  opacus, 

transparent  or 

transparent;  through 

progressively  in- 

opaque),  progres- 

this  part  the  sun  or 

vading  the  sky; 

sively  invading  the 

moon  may  be  weakly 

these  Altocumu- 

sky;  these  Altocumu- 

visible  as  through 

lus  clouds  gen- 

lus  clouds  generally 

ground  glass. 

erally  thicken  as 

thicken  as  a whole. 

a whole. 

2 

Altostratus  opacus 

Altostratus,  the  greater 

or  Nimbostratus. 

part  of  which  is  suffi-  ® 

Altocumulus  cumu- 

Altocumulus  resulting 

ciently  dense  to  hide 

logenitus  (or 

from  the  spreading 

the  sun  or  moon,  or 

cumulonrnibo- 

out  of  Cumulus  (or 

Nimbostratus. 

genitus). 

Cumulonimbus) . 

3 

Altocumulus  trans- 

Altocumulus,  the  ^ 

Altocumulus  trans- 

Altocumulus  in  two  or 

lucidus  at  a 

greater  part  of  which 

lucidus  or  opacus 

more  layers,  usually 

single  level. 

is  semitransparent, 

in  two  or  more 

opaque  in  places,  and 

the  various  elements 

layers,  or  Alto- 

not  progressively 

of  the  cloud  change 

cumulus  opacus 

invading  the  sky ; or 

only  slowly  and  are 

in  a single  layer 

opaque  layer  of  Alto- 

all  at  a '.ingle  level. 

not  progressively 

cumulus,  not  pro- 

t 

invading  the  sky, 

gressively  invading 

or  Altocumulus 

the  sky;  or  Alto- 

4 

Patches  (often  len- 

Patches  (often  in  the 

with  Altostratus 

cumulus  together 

ticular)  of  Alto- 

form  of  almonds  or 

or  Nimbostratus. 

with  Altostratus  oi 

cumulus  translu- 

fishes)  of  Altocumu- 

Nimbostratus. 

cidus,  continually 

lus,  the  greater  part 

changing  and 

of  which  is  semi-  8 

Altocumulus  castel 

Altocumulus  with 

occurring  at  one 

transparent;  the 

lanus  or  floccus. 

sproutings  in  the  form 

or  more  levels. 

clouds  occur  at  one 

of  small  towers  or 

or  more  levels  and  the 

battlements,  or  Alto- 

elements  are  continu- 

cumulus  having  the 

ally  changing  in 

appearance  of  cumu 

1 

appearance. 

liform  tufts. 

9 

Altocumulus  of  a 

Altocumulus  of  a 

chaotic  sky,  gen- 

chaotic  sky,  gen- 

erally  at  several 

erally  at  several 

levels. 

levels. 
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Table  B-II-4.  High  Cloud-WMO  Code 


1 


Symbol  CH=Clouds  of  Genera  Ci,  Cc,  Cs 


Code 

Fig- 

ure 

Technical  Specifica- 
tions 

Nontechnical  Specifi- 
cations 

0 

No  Ca  clouds 

No  Cirrus,  Cirrocumu- 
lus,  or  Cirrostratus. 

1 

Cirrus  fibratus, 

Cirrus  in  the  form  of 

sometimes  unci- 

filaments,  strands  oi 

nus,  not  progres- 

hooks,  not  progres- 

sively  invading 

sively  invading  the 

the  sky. 

sky. 

2 

Cirrus  spissatus,  in 

Dense  Cirrus  in  patches 

patches  or  entan- 

or  entangled  sheaves, 

gled  sheaves, 

which  usually  do  not 

which  usually  do 

increase  and  some- 

not  increase  and 

times  seem  to  be  the 

sometimes  seem 

remains  of  the  upper 

to  be  the  remains 

part  of  a Cumulo- 

of  the  upper  part 

nimbus;  or  Cirrus 

of  a Cumulouiin- 

with  sproutings  in 

| bus;  or  Cirrus 

in  the^form  of  small 

castellanus  or 

turrets  or  battle- 

! floccus. 

ments,  or  Cirrus 
having  the  appear- 
ance of  cumuliform 

tufts. 

3 

Cirrus  spissatus  cu 

Dense  Cirrus,  often  in 

mulonimbogeni- 

the  form  of  an  anvil, 

tus. 

being  the  remains  of 
the  upper  parts  of 
‘ Cumulonimbus. 

4 

Cirrus  uncinus  or  fi- 

! Cirrus  in  the  form  of 

bratus,  or  both, 

hooks  or  of  filaments 

progressively  in- 

or  both,  progressively 

vading  the  sky; 

invading  the  sky; 

they  generally 

they  generally  be 

thicken  as  a 

come  denser  as  a 

whole. 

whole. 

5 

Cirrus  (often  in 

Jirrus  (often  in  bands 

bands)  and  Cirro- 

converging  towards 

stratus,  or  Cirro- 

one  or  two  opposite. 

stratus  alone,  pro- 

points  of  the  hori- 

gressively  invad- 

zon)  and  Cirrostra- 

ing  the  sky;  they 

tus.  or  Cirrostratus 

generally  thicken 

alone;  in  either  case. 

as  a whole,  but 

they  are  progres- 

the  continuous 

sively  invading  the 

veil  does  not 

sky,  and  generally 

I 

reach  45°  above 

growing  denser  as  a 

the  horizon. 

whole,  but  the  con- 

tinuous  veil  does 
not  reach  45s  above 
the  horizon. 


Code 

Fig- 

ure 

Technical  Specifica- 
tions 

| Nontechnical  Specifi- 
cations 

| 

i 

6 

Cirrus  (often  in 

Cirrus  (often  in  bands 

bands)  and  Cirro- 

converging  towards 

stratus,  or  Cirro- 

one  or  two  opposite 

stratus  alone,  pro- 

points  of  the  hori- 

1 

gressively  invad- 

zon)  and  Cirro- 

ing  the  sky;  they 

stratus,  or  Cirro- 

generally  thicken 

stratus  alone;  in 

as  a whole;  the 

either  case,  they  are 

continuous  veil 

propressively  invad- 

extends  more 

ing  the  sky,  and 

than  45s  above 

generally  growing 

the  horizon,  with- 

denser  as  a whole; 

out  the  sky  being 

the  continuous  veil 

totally  covered. 

extends  more  than 

45*  above  the  hori- 
zon, without  the  sky 

being  totally  covered. 

7 

Cirrostratus  cover-  | 

Veil  of  Cirrostratus 

ing  the  whole  sky  ! 

covering  the  celes- 
tial dome. 

8 

Cirrostratus  not 

Cirrostratus  not  pro- 

progressively  in- 

gressively  invading 

vading  the  sky 

the  sky  and  not 

and  not  entirely  | 

completely  covering 

covering  it. 

the  celestial  dome. 

Cirrocumulus  alone, 

Cirrocumulus  atone,  or 

or  Cirrocumulus 

Cirrocumulus  accom- 

predominant 

panied  by  Cirrus  or 

among  the  Ch 

Cirrostratus,  or  both. 

clouds. 

but  Cirrocumulus  is 

predominant. 
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NRL  Aerosol  Related  Instrumentation  Descriptions 


The  suite  of  instruments  which  was  assembled  by  NRL  to  determine 
the  characteristics  of  the  atmosphere  by  means  of  measurement  of  various 
physical  parameters  of  the  natural  atmospheric  aerosol  constituency  is 
described  below: 


1)  Condensation  nuclei  (CN)  density  were  precisely  obtained  by 
means  of  a Poliak  Counter  which  has  its  calibration  traceable  back  to  the 
Poliak  original  calibration.  This  data  is  found  in  file  CN.DAT  column  2. 

The  condensation  nuclei  density  distribution  about  the  ship  was  ob- 
tained with  a portable  G.E.  condensation  counter.  This  work  is  presented 
in  section  C-IV  of  this  report. 

2)  Continuous  observations  of  condensation  nuclei  were  also  ob- 
tained by  means  of  the  Environment  One  Condensation  Nuclei  Monitor,  model 
Rich  100.  Although  there  is  some  question  as  to  whether  this  instrument 
was  operating  properly,  its  data  is  stored  in  file  ASC.DAT(7). 

3)  Cloud  Condensation  Nuclei  (CCN)  counts  were  obtained  through- 
out the  cruise  by  means  of  the  NRL  CCN  equipment  (see  R.  E.  Ruskin  and 

J.  E.  Dinger,  "Descriptions  of  NRL  CCN  Equipment  in  The  Second  Interna- 
tional Workshop  on  Condensation  and  Ice  Nuclei,"  compiled  by  Lewis  0. 
Grant,  Dept,  of  Atmos.  Science,  Colorado  State  University,  Fort  Collins, 
CO,  May  1971).  This  data  is  stored  in  file  CCN. DAT  where: 

Column  2 is  the  CCN/cc  at  AT  = 2.0°C 


Column  4 Is  the  CCN/cc  at  aT  = 4.3°C 
Column  5 is  the  CCN/cc  at  aT  = 5.0°C 

Column  6 is  the  percentage  of  involatile  particles  in  the  sample  at 


a AT  * 4.3°C. 

Values  of  88888  mean  that  data  was  not  available  for  this  number  and 
values  of  99999  mean  that  very  high  numbers  of  particles  were  observed, 
i.e.,  numbers  too  high  to  accurately  count. 


4)  Scattering  coefficients  were  measured  with  the  MRI  Model 
1567  Integrating  Nephelometer.  This  device  was  located  on  a tower  at 
location  B in  Figure  B-II-1.  The  units  of  this  data  are  expressed  as 
10"4  iif1.  This  data  is  in  ASC.DAT(4). 

5)  A Royco  model  225  aerosol  sizing  device  was  used  by  NRL 
throughout  the  cruise.  This  device  divides  the  aerosol  size  distribution 
into  5 parts.  Thus,  for  every  hour  we  have  represented  in  these  data 
files  the  10  minute  mean  concentrations  of  particles  which  fall  within 

3 

the  five  siie  classes.  Numbers  used  are  the  number  of  particles  per  cm 
in  the  following  ranges  of  particle  diameter,  d: 


Size  of  Channel 
0.45  nn  £ d < 0.60  pm 
0.60  un  £ d < 1.50  un 

1.50  un  £ d < 2.50  un 

2.50  pm  £ d < 4.0  jin 
d > 4.0  ym 


Location 
R0Y.DAT(2) 
ROY .DAT ( 3 ) 
R0Y.DAT(4) 
R0Y.DAT(5) 
ROY . DAT ( 6 ) 
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These  values  are  the  "true"  ambient  particle  concentrations  and  were  ob- 
tained from  the  measured  particle  concentrations  using  the  formula  of 
Balyaev  and  Levin  (1972)  which  gives  the  particle  sampling  efficiency  of 
a tube  facing  into  the  ambient  air  flow  as  a function  of  particle  size 
and  the  ratio  of  the  ambient  wind  speed  to  the  airspeed  in  the  sampling 
tube. 

6)  There  were  two  new  experimental  aerosol  instruments  used  on 
this  cruise.  These  principles  are  not  widely  known  and  therefore  a more 
detailed  description  of  these  devices  and  the  data  obtained  will  be  found 
in  separate  reports  in  Sections  C-I  and  C-II  of  this  report. 


TABLE  B— 1 1-5 


Meteorological  Data  Files 


Parameter 

Units 

Location 

1) 

Sea  Water  Temperature  at  7m 

(°C) 

THERM.DAT(2) 

2) 

Air  Temperature  at  19m 

(°C) 

THERM.DAT(3) 

3) 

Air  Temperature  at  19m 

(°C) 

F1LE.DAT(2) 

4) 

Dewpoint  Temperature  at  19m 

(°C) 

FILE.DAT(3) 

5) 

Relative  Wind  Speed 

(knots) 

ASC.DAT(2) 

6) 

Relative  Wind  Direction 

(degrees) 

ASC.DAT(3) 

7) 

Atmospheric  Radon 

Concentration 

(pCi/cc) 

RAD0N2(2) 

8) 

Calculated  Air  Mass 

V0LTZ.DAT(2) 

9) 

Aerosol  Thickness  (x  = 0.88u) 

TURBID. DAT(2) 

10) 

Aerosol  Thickness  (x  = 0.50u) 

TURBID. DAT(3) 

11) 

Aerosol  Thickness  (x  = 0.94y) 

TURBID. DAT(4) 

12) 

Ship's  Latitude 

(whole 
degrees ) 

NAV.DAT(2) 

13) 

Ship's  Latitude 

(decimal 
fraction 
of  degree) 

NAV.DAT(3) 

14) 

Ship's  Longitude 

(whole 

degree) 

NAV.DAT(4) 

15) 

Ship's  Longitude 

(decimal 
fraction 
of  degree) 

NAV.DAT(5) 

16) 

Ship's  Course  Made  Good 

(degrees) 

NAV.DAT(6) 

17) 

Ship's  Speed 

(knots) 

NAV.DAT(7) 

18) 

Sea  Skin  Temperature  PRT-5 

(°C) 

FILE.DAT(4) 

19) 

Sea  Surface  Temperature  PRT-10 

(°C) 

PRT10.DAT(2 ) 

20) 

Sky  Temperature  PRT-10 

(°C) 

PRT10.DAT(3) 
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TABLE  B-II-5  (Contd) 


Parameter 

Units 

Location 

21) 

PRT-10  (-10°) 

(°C) 

PRT10.DAT{4) 

22) 

PRT-1Q  (+10°) 

(°C) 

PRT10.DAT(5) 

23) 

PRT-10  (+30°) 

(°C) 

PRT10.DAT(6) 

24) 

PRT-10  (+60°) 

(°C) 

PRT10.DAT(7) 

25) 

Cloud  Base  Temperature 

(°C) 

PRT10.DAT(8) 

26) 

Sea  State 

CLOUD. DAT(2) 

27) 

White  Cap/Foam  Cover 

(%) 

CL0UD.DAT(3) 

28) 

Visibility 

(mi ) 

CL0UD.DAT(4) 

29) 

Cloud  Cover 

CLOUD. DAT( 5) 

30) 

Low  Cloud  Types 

CL0UD.DAT(6) 

31) 

Middle  Cloud  Types 

CL0UD.0AT(7) 

32) 

High  Cloud  Types 

CL0UD.DAT(8) 

33) 

Poliak  Counter  Condensation 
Nuclei  Density 

(#/cc) 

CN.0AT(2) 

34) 

E.O.C.N.M.-CN  Data 

ASC.DAT{7) 

35) 

CCN/cc  at  AT  2.0°C 

CCN.DAT(2) 

36) 

CCN/cc  at  AT  3.0°C 

CCN .DAT ( 3 ) 

37) 

CCN/cc  at  AT  4.3°C 

CCN.DAT(4) 

38) 

CCN/cc  at  AT  5.0°C 

CCN.DAT(5) 

39) 

(%)  Involatile  Particles  in 
Sample  at  AT  = 4.3°C 

CCN.DAT(6) 

NRL 

ROYCO  DATA 

40) 

Channel  Size  = d, 

0.45  ini  < d < 0.60  ym 

(#/cc) 

R0Y.DAT{2) 

41) 

0.60  ym  < d < 1.50  ym 

(#/cc) 

ROY .DAT ( 3 ) 
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TABLE  B-II-5  (Contd) 


Parameter 

Units 

Locati on 

42) 

1.50  ym  _<  d < 2.50  yin 

(#/cc) 

R0Y.DAT(4) 

43) 

2.50  ym  < d < 4.0  ym 

(#/cc) 

R0Y.DAT{5) 

44) 

d > 4.0  ym 

(#/ cc) 

R0Y.DAT(6) 

MEASUREMENTS  OF  AEROSOL  AND  MICROMETEOROLOGICAL 


CHARACTERISTICS  OF  THE  MARINE  BOUNDARY  LAYER 
DURING  THE  TRANSATLANTIC/MEDITERRANEAN  CRUISE 
OF  THE  USNS  HAYES,  MAY-JUNE  1977 

by 

Ulrich  Katz  and  Eugene  J.  Mack 
Calspan  Corporation,  Buffalo,  New  York  14221 

SUMMARY 

Under  contract  no.  N00173-77-C-0126  with  the  Naval  Research  Labora- 
tory (NRL),  Calspan  Corporation  provided  meteorological  and  aerosol  phys- 
ics support  during  the  NRL  cruise  77-16-04  aboard  the  USNS  HAYES.  Data 
were  acquired  during  the  period  from  14  May  to  7 June  1977  while  en  route 
from  Virginia  to  Athens,  Greece  via  the  Grand  Banks  and  Gibraltar.  The 
objective  of  Calspan' s participation  in  the  study  was  to  record  physical 
and  chemical  characteristics  of  the  aerosol  (including  fog)  in  the  lower 
atmospheric  marine  boundary  l^yer  as  well  as  the  pertinent  basic  meteoro- 
logical parameters  for  subsequent  interpretive  analyses. 

Calspan  instrumentation  deployed  aboard  the  USNS  HAYES  is  listed  in 
Table  B-II-6  along  with  information  on  parameters  measured,  range  of  vari- 
ables, frequency  of  measurements  and  measurement  heights  above  the  sea 
surface.  The  locations  of  the  various  instruments  on  the  ship  are  depic- 
ted schematically  in  Figure  B-1I-2  . The  first  five  items  in  Table  B-II-6 
were  installed  in  a shelter  which  formed  the  lower  half  of  a small  tower 
erected  atop  the  flying  bridge.  In  order  to  minimize  the  risk  of 
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contamination,  the  sample  air  tor  the  various  particle  monitors  was  as- 
pirated through  a 5 cm  diameter  (ID),  2.5  m long  tube  extending  to  2 m 
forward  of  and  2 m above  the  flying  bridge.  (Flow  rate  was  such  that 
aerosol  residence  within  the  tube  was  2 sec.)  Fog  sampling  equipment 
(items  6,  7,  and  8)  was  mounted  on  top  of  the  tower  to  minimize  the 
ship's  thermal  influence  when  the  wind  direction  is  within  the  forward 
quadrant.  A second  drop  sampler  was  positioned  on  the  top  of  one  of  the 
bows  (depending  on  wind  direction)  for  collecting  sea  spray  droplets  as 
close  to  the  water  surface  as  possible  while  escaping  any  spray  generated 
by  the  ship.  Continuous  temperature  recordings  were  obtained  at  four 
levels,  the  lowest  being  the  sea  surface.  (By  towing  an  armored  sea  sur- 
face temperature  probe  between  the  bows,  it  was  possible  to  make  the 
probe  follow  within  the  top  10  centimeters  of  the  water  surface.)  Con- 
tinuous humidity  measurements  were  obtained  at  two  levels  with  lithium 
chloride  based  dew  point  sensors  while  manual  psychrometry  provided  cali- 
bration and  backup  measurements.  The  main  purpose  of  the  wind  instrumen- 
tation mounted  on  the  antenna  mast  was  the  convenient  detection  of  poten- 
tially contaminating  wind  conditions.  While  the  last  two  items  in  Table 
B-II-6  provided  bulk  aerosol  samples  for  later  chemical  analysis,  the  cas- 
cade impactor  samples  are  most  suitable  for  combined  scanning  electron 
microscopy  and  energy  dispersive  X-ray  analysis  of  individual  particles, 
whereas  the  bulk  filter  samples  are  partly  earmarked  for  wet  chemical 
analyses  and  partly  designed  for  atomic  absorption  spectroscopy  or  X-ray 
fluorescence  analysis. 
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The  data  are  presently  available  in  the  form  of  hourly  logs,  strip 
chart  records  and  unprocessed  samples  from  the  thermal  diffusion  chamber, 
droplet  impactor,  cascade  impactor  and  hi-volume  filters.  The  logs  com- 
prise hourly  values  of  all  the  parameters  measured  except  CCN,  sea  spray 
and  chemical  samples  which  were  obtained  at  irregular  intervals  and  which 
require  further  processing  to  arrive  at  quantitative  values.  In  total, 
the  logs  comprise  approximately  400  successful  hourly  observations.  An 
additional  100  readings  were  contaminated  (aerosol  data  mainly)  by  ship's 
exhaust  resulting  from  unfavorable  winds.  Table  B-II-7  summarizes  the 
quantities  and  types  of  aerosol  measurements  which  were  obtained  at  dis- 
crete intervals.  Visibilities,  scattering  coefficient,  temperatures,  dew 
points,  and  winds  were  recorded  continuously  and  are  available  on  strip 
charts  as  well  as  in  the  logs. 

The  data  are  grouped  into  two  separate  hourly  logs,  one  containing 
the  meteorological  parameters  (Table  B-II-8),  the  other  listing  aerosol- 
related  readings  (Table  B-II-9).  It  should  be  noted  that  the  Royco  val- 
ues entered  in  the  log  are  counts  for  the  10-minute  period  centered  about 
the  indicated  time;  printouts  of  five  additional  counts  for  every  hour 
are  on  file  at  Calspan. 

Since  the  raw  EAA  data  required  some  processing,  the  aerosol  data 
archived  in  the  data  system  and  presented  in  Section  B-IV-d  of  this  re- 
port are  the  results  of  computer  analysis  done  at  Calspan. 
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TABLE  B-II-7 


SUMMARY  OF  CALSPAN  DISCRETE  AEROSOL  MEASUREMENTS 


Approximate 
Number  of 
Measurements 

Particles,  Instruments 

Period 

420 

Aitken,  Gardner 

Hourly 

390 

Small  p.,  EAA 

Hourly 

1600 

Large  p.,  Royco 

10  minute  averages. 

Some  gaps  (printer 
malfunction) 

90 

CCN,  Therm.  Diff.  Ch. 

3-4  hours 

20 

Sea  spray.  Drop  Impactor 

Irreg.,  25  May  - 6 June 

35 

Large  and  giant  p.. 

Cascade  Impactor 

Irreg.,  17  May  - 7 June 

12 

All  p. , Hi-Vol  filters 

Irreg.,  8-20  hours 

TABLE  B-II-8 

CALSPAN'S  METEOROLOGICAL  LOG 


Parameter 

Units 

Fil e/Column  No 

1) 

Relative  Wind  Direction 

(deg)  : 

CALVIS.DAT(2) 

2) 

Relative  Wind  Speed 

(mph)  : 

CALVIS.DAT(3) 

3) 

EG&G  Visibility 

(m) 

CALVIS.DAT (4) 

4) 

MRI  Visibility 

(mi)  : 

CALVIS.DAT(5) 

5) 

Psychrometer  Dry  Bulb 

(°F)  : 

CALVIS.DAT(6) 

6) 

Psychrometer  Wet  Bulb 

(°F)  : 

CALVIS.DAT(7) 

7) 

Air  Temperature  at  27  m 

(°C)  : 

CALTMP.DAT(2) 

8) 

Air  Temperature  at  15  m 

(°C)  : 

CALTMP.DAT(3) 

9) 

Air  Temperature  at  9 m 

(°C)  : 

CALTMP . DAT ( 4 ) 

10) 

Sea  Surface  Temperature 

(°C)  : 

CALTMP.DAT(5) 

ID 

Dewpoint  at  27  m 

(°C)  : 

CALTMP .DAT(6) 

12) 

Dewpoint  at  15  m 

(°C)  : 

CALTMP. DAT( 7) 

13) 

Ship’s  Magnetic  Heading 

(deg)  : 

CALTMP. DAT( 8) 
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TABLE  B-I 1-9 


CAL SPAN ‘ S AEROSOL  LOG 


Parameter 


Tilel Column  No.) 


1) 

Aitken  Nuclei  Concentration  (#/cc)  : 

OU  T . AER ( 2 ) 

2) 

Royco  Data  in  #/2. 

a) 

>0 . 3 n 

0UT.AER(3) 

b) 

>0.6y  I 

0UT.AER(4) 

c) 

>1 . 2 u 

0UT.AER(5) 

d) 

>3.  y i 

OUT . AER ( 6 ) 

e) 

>5 . y : 

OUT . AER ( 7 ) 

3) 

bscat  f'"'4"'1'  : 

OUT.AER(B) 

4) 

Thermo  Systems  Model  3030 

Electrical  Aerosol  Analyzer  N = #/cc 

a) 

Log  N (r>.0032y) 

WH I T2 -DAT ( 2 ) 

b) 

Log  N (r>.0056y) 

WH I T2 .DAT ( 3 ) 

c) 

Log  N ( r> -01 u) 

WHIT2.DAT(4) 

d) 

Log  N (r>.U17By) 

WH I T2 -DAT ( 5 ) 

e) 

Log  N (r>.u316y) 

WHIT2 .DAT(6) 

f) 

Log  N (r>.0562y) 

WHIT2.DAT(7) 

g) 

Log  N (r>.ly) 

WH1T2.DAT{8) 

h) 

Log  N (r>.178y) 

WHIT2.DAT( 9) 

i) 

Log  N (r>.360y) 

WHIT2.DAT( 10) 

j) 

Log  N (r>.560y) 

WHIT2 .DAT( 11 ) 

DSNS  HAYES  Marine  Boundary  Layer 
Research  Cruise: 

Preliminary  Evaluation  of  NPS  Data 


C.  W.  Fairall,  L.  F.  May,  K.  L.  Davidson 
T.  Houlihan  and  G.  E.  Schacher 

Environmental  Physics  Group 
Naval  Postgraduate  School 
Monterey,  CA  93940 

In  May  and  June  of  1977  the  Naval  Air  Systems  Command  (AIR  370)  spon- 
sored a marine  boundary  layer  research  cruise  aboard  the  Naval  Research 
Laboratory  (NRL)  ship,  USNS  HAYES,  in  the  North  Atlantic  and  Mediterra- 
nean. The  cruise  involved  scientific  personnel  from  several  laboratories 
including  the  Naval  Postgraduate  School  (NPS).  We  departed  Chetham  Annex, 
Virginia  on  May  15,  paralleled  the  Coast  to  Nova  Scotia,  then  took  a great 
circle  route  across  the  Atlantic  arriving  in  Rota,  Spain  on  May  27.  We 
departed  Rota,  Spain  on  May  30  and  arrived  in  Piraeus,  Greece  on  June  7. 

NPS  participated  in  the  cruise  primarily  to  gather  atmospheric  tur- 
bulence data  for  a wide  range  of  atmospheric  and  sea  state  conditions. 
These  data  are  used  to  determine  turbulent  fluxes  of  heat,  momentum,  and 
water  vapor  and  turbulent  mixing  of  other  atmospheric  scalars  such  as 
aerosols,  electric  change  and  pollutants.  The  data  are  also  important 
for  characterizing  the  optical  propagation  qualities  of  the  atmosphere. 


This  is  a preliminary  report  on  the  basic  mean  and  turbulence  data 
gathered  by  the  NPS  personnel. 


II.  USNS  HAYES  INSTALLATION 

The  USNS  HAYES  (T-AGOR  16)  is  a twin  hull  ship  that  was  constructed 
for  NRL  for  use  as  an  oceanographic  research  vessel.  Pertinent  dimensions 
are:  length  overall  246  ft,  beam  75  ft,  and  distance  between  hulls  27  ft. 
Due  to  the  catamaran  construction  the  ship  has  a very  broad  forward  wall 
beginning  approximately  50  ft  back  from  the  bows  of  the  ship,  which  ex- 
tends nearly  the  full  width  of  the  ship  and  to  a height  of  approximately 
56  ft  above  the  mean  water  line.  The  presence  of  this  wall  greatly  per- 
turbs the  natural  airflow  and  causes  problems  in  the  interpretation  of  our 
data,  as  will  be  described. 

Figure  B- 1 1-3  shows  the  profile  of  the  ship  and  the  location  of  the 
two  stations  at  which  NPS  installed  sensors.  We  welded  a 12  ft  tubular 
steel  tower  as  far  forward  on  the  port  bow  as  possible  (approximately  2 ft 
back  from  the  tip  of  the  bulwark).  The  mean  sensors  were  placed  at  the 
top  of  this  tower  and  the  turbulence  sensors  were  mounted  on  a wind  vane, 
which  was  mounted  on  a 4 ft  extension  at  the  top  of  the  tower.  The  vane 
is  used  to  keep  the  sensors  pointing  into  the  wind,  and  the  extension 
places  the  sensors  forward  of  the  tower  so  that  shipboard  influence  is 
eliminated.  The  second  station  was  located  on  the  first  level  of  the 
ship's  forward  mast.  Again  the  turbulence  sensors  were  mounted  on  a vane 
on  an  extension  to  keep  these  sensors  as  far  forward  of  the  platform  as 
was  feasible.  The  arrangement  at  this  station  was  such  that  the  mean  sen- 
sors were  approximately  2 ft  below  the  turbulence  sensors.  The  heights  of 
the  various  sensors  above  the  mean  water  level  are  found  in  Table  B-II-10. 
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TABLE  B- 1 1— 10  Height  of  NPS  Sensors  Above  Mean  Sea  Level 


STATION 

SENSOR  ht. 

EFFECT,  ht. 

(tunnel ) 

EFFECT. 

ht  (meas.) 

Mean 

Turb 

Mean 

Turb 

Mean 

Turb 

BOW 

12.1m 

12.1m 

11.3m 

11.3m 

11.0m 

11.0m 

MAST 

23.8m 

24.5m 

17.6m 

18.5m 

19.6m 

20.5m 

Figure  B-II-3  shows  a schematic  representation  of  the  air  flow  over 
the  ship  when  the  relative  wind  is  from  the  bow.  The  effect  of  the  flow 
is  that  the  air  sampled  by  the  sensors  at  the  mast  station  comes  from  a 
height  lower  than  their  actual  height  above  the  sea  surface.  The  location 
of  the  bow  station  is  such  that  the  air  sampled  comes  from  nearly  the  same 
height  as  the  sensors. 


Before  the  1975  fog  cruise,  NRL  had  a wind  tunnel  analysis  of  the  air 
flow  dround  the  Hayes  performed.  Also,  during  that  cruise.  Dr.  Richard 
Jeck  of  NRL  made  a number  of  measurements  of  the  wind  velocity  in  the  for- 
ward region  of  the  ship  with  a bivane  anemometer.  Evaluation  of  these  two 
sets  of  data  do  not  give  us  as  accurate  information  as  we  need,  but  they 
do  show  that  the  air  arriving  at  the  mast  rises  approximately  20  ft.  This 
correction  is  applied  to  the  height,  giving  the  air  height  listed  in  Table 
B-1I-10  as  effective  height  (tunnel).  Note  that  the  air  flow  information 
that  is  available  is  only  for  the  relative  wind  directly  from  the  bow  at 
about  10  knots.  We  have  made  a better  estimate  of  the  effective  heights 
from  the  ratio  of  turbulent  dissipation,  e,  at  the  two  heights.  In  the 
boundary  layer. 


ix/i2  « (e2/ei)  fc  (yu/fe  (Z2/L) 


where  f (Z/L)  is  a stability  correcting  factor  and  L is  the  Monin-Obukhov 
length  (a  typical  value  for  the  Northern  oceans  if  L = 100  meters).  For 
this  cruise  we  have  230  simultaneous  values  of  and  e2  giving 

< e]/ ^2  > = 
which  results  in 

= 1.86. 

Assuming  the  bow  level  is  much  less  affected  than  the  mast  we  have  as- 
signed the  values  of  Z j and  Z2  given  in  the  third  column  of  Table  B-II-10. 

III.  NPS  Equipment 

1)  General  Information 

We  installed  equipment  to  measure  both  mean  and  turbulence  param- 
eters at  both  the  bow  and  mast  stations.  The  parameters  measured  were: 

a.  humidity 

b.  temperature 

c.  horizontal  wind  speed 

d.  temperature  fluctuations 

e.  horizontal  wind  speed  fluctuations 

A sea  surface  temperature  sensor  was  installed  between  the  hulls  about 
75  ft  forward  of  the  stern.  Relative  wind  direction  was  measured  at  the 
port  bow  level. 
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Details  of  the  equipment  used  for  the  measurements  follows. 

The  sensors  for  both  mean  humidity  and  mean  temperature  were  mounted 
in  a single  aspirator  at  each  station.  They  were,  therefore,  protected 
from  the  weather  while  a steady  flow  of  air  insured  that  they  remained  at 
ambient  conditions. 

2)  Humidity 

A Hygrodynamics  Digital  II  system  employing  15-1818  sensors  was 
used.  The  sensors  contain  lithium  chloride  cells  and  thermistors  so  that 
humidity  and  temperature  were  both  measured  continuously.  This  tempera- 
ture sensor  was  used  only  as  a monitor,  not  for  our  mean  temperature  mea- 
surements. The  humidity  sensors  were  calibrated  in  a humidity  chamber, 
the  results  are  shown  in  Table  B-II-11.  Sensor  #2  was  used  on  the  bow 
level  and  Sensor  #6  was  used  on  the  mast. 

3)  Mean  Temperature 

Hewlett  Packard  Model  2801A.  In  this  system  the  sensor  is  a 
quartz  crystal  which  forms  the  capacitative  element  in  an  oscillator  tank 
circuit.  Temperature  changes  change  the  oscillator  frequency,  which  is 
sensed  for  the  temperature  readout.  This  system  is  quite  accurate,  allow- 
ing temperatures  to  be  determined  within  O.U3°C.  Sea  surface  temperature 
is  measured  with  the  same  system. 
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Table  B-II-11  NPS  Humidity  Sensor  Calibration 
Hygrodynamics  Digital  II  Calibration 


Wet- 

Dry 

Bulb 

%RH 

Ave 

Hygro 

%RH 

1 

Deviation  of 

2 3 

Sensor 

4 

from  Ave  Hygro  % RH 

5 6 7 

8 

33.0 

31.6 

1.4 

-1.4 

-0.6 

1.5 

-0.2 

-1.6 

-0.9 

1.8 

36.0 

37.2 

-0.5 

-0.8 

-0.3 

1.6 

0.5 

-1.4 

-0.1 

1.2 

44.0 

41.2 

-0.8 

-0.4 

-0.3 

0.6 

0.9 

0.4 

-0.1 

0.3 

NA 

41.5 

-0.5 

-0.4 

-0.3 

0.5 

0.9 

-0.4 

-0.1 

-0.1 

NA 

42.6 

-1.2 

-0.5 

-0.5 

-0.9 

1.8 

-0.4 

-0.2 

0.1 

NA 

65.8 

-0.2 

-1.4 

-2.5 

0.8 

1.8 

-0.3 

1.3 

0.3 

74.8 

73.6 

0.2 

0.2 

-1.3 

0.5 

-0.2 

0.4 

-0.4 

0.8 

88.3 

85.0 

-0.7 

0.9 

-0.4 

0.1 

-0.6 

0.0 

0.4 

0.2 

82.6 

86.7 

-1.4 

-0.6 

-0.5 

0.1 

-0.2 

0.3 

0.7 

0.3 

82.6 

87.0 

-1.0 

0.5 

-0.7 

0.0 

0.3 

0.7 

0.0 

0.2 

85.0 

88.0 

-0.4 

1.1 

-0.1 

0.1 

-1.8 

1.1 

-0.6 

-0.1 

93.9 

89.6 

0.1 

0.8 

-0.2 

0.3 

-2.2 

2.1 

-1.5 

0.4 

94.4 

92.3 

-0.6 

1.1 

-0.5 

-0.3 

-0.4 

1.2 

-0.5 

0.4 

94.4 

93.6 

-0.8 

1.0 

-0.8 

-0.3 

0.4 

-0.2 

0.8 

-0.4 

Ave  Devi  at 

-0.5 

-0.0 

-0.6 

0.3 

0.1 

0.1 

-0.1 

0.4 

a 

0.7 

0.9 

0.6 

0.7 

1.1 

1.0 

0.7 

0.6 

The  average  deviation  of  the  wet-dry  bulb  from  the  Ave  Hygro  reading  is 
0.4%  RH  with  a standard  deviation  of  2.9%  RH. 
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4)  Mean  Wind  Speed 


Thornthwaite  Model  101  Wind  Register  System.  This  is  a photo- 
electric system  with  very  light  plastic  cups  that  allow  measurements  of 
winds  as  slight  as  one  half  knot. 

5)  Temperature  Fluctuations 

A low  power  ac  Wheatstone  bridge  was  constructed  using  a GTE 
Sylvania,  Inc.  Model  140  Lightweight  Thermo-sonde  System  with  TSI  Model 
1210  probes  with  P 8 platinum  wire  used  as  sensing  elements.  The  low 
power  bridge  was  operated  at  100  u Amp,  the  resistances  of  the  probes 
were  approximately  40  Ohms,  thus,  the  energy  dissipation  in  the  sensor  is 
5 x 10  7 Watts.  Very  low  power  was  used  to  insure  that  the  sensor  is  not 
elevated  above  ambient  temperature  which  prevents  velocity  fluctuations 
from  influencing  the  wire  response.  The  wire  is  1.2mm  long  and  2.5p  in 
diameter,  and  the  probe  has  a time  constant  much  shorter  than  the  times 
encountered  in  atmospheric  turbulence. 

6)  Wind  Speed  Fluctuations 

TSI  Model  1054B  Anemometer,  TSI  Model  1210  probes  with  6 mil.  hot 
film  elements  were  used  as  the  sensing  elements.  An  overheat  ratio  of  1.2 
was  utilized  so  that  the  wire  temperature  was  approximately  150°C.  The 
overheat  greatly  diminishes  the  influence  of  temperature  fluctuations  on 
these  measurements. 

The  axes  of  the  films  were  aligned  in  the  vertical  direction  so  that 
the  wires  were  not  sensitive  to  air  flow  in  the  vertical  direction.  Thus, 
only  the  horizontal  component  of  wind  speed  fluctuations  is  detected. 
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7)  Data  Recording  Equipment 

Fluctuation  data  for  temperature  and  wind  speed  were  recorded  on 
1/4  inch  magnetic  tape. 

Mean  wind  data  was  recorded  by  hand.  The  mean  temperature  and  humid- 
ity were  printed  on  paper  tape  with  a Hewlett  Packard  562AR  printer. 

In  order  to  correctly  understand  the  mean  temperature  data,  it  is 
necessary  to  describe  the  data  gathering  sequence  in  some  detail.  There 
are  five  signals  to  record:  humidity  at  two  levels  and  temperature  at 
three  levels,  including  the  sea  surface.  These  signals  are  processed  by 
a homemade  sequencer  which  sorts  and  sends  the  signals  to  the  printer. 

The  sequencer  steps  cyclically  from  level  to  level  in  the  sequence:  sea 
surface  to  bow  to  mast,  and  the  temperature  and  humidity  at  a given  level 
are  printed  simultaneously.  When  the  sequencer  is  at  the  sea  surface 
step,  zeros  are  printed  for  the  humidity  and  this  allows  us  to  identify 
the  levels  on  the  print  tape. 

The  interfacing  between  the  sequencer  and  the  Hewlett  Packard  tem- 
perature readout  posed  some  problems.  Since  the  readout  is  a frequency 
counter,  it  has  its  own  count  cycle,  the  timing  of  which  is  set  by  a front 
panel  control.  It  was  not  possible  to  control  the  readout  timing  by  a 
command  from  the  sequencer  and  we  did  not  construct  the  sequencer  to  ac- 
cept a command  from  the  readout.  Thus,  it  is  possible  for  the  sequencer 
to  issue  a print  command  when  the  counter  is  in  a read  cycle,  leading  to 
an  error.  We  set  the  readout  timer  In  a way  so  that  such  errors  are  very 
infrequent  and  they  are  normally  easily  identified. 
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The  three  mean  temperatures  and  the  two  mean  humidities  were  read 
and  recorded  every  2.5  minutes  throughout  the  cruise  except  when  a mal- 


function occurred.  Fluctuation  data  was  recorded  at  hourly  intervals,  as 
conditions  warranted.  The  fluctuating  signals  were  also  continually  pro- 
cessed by  filtering  and  RMS  units  and  recorded  on  strip  chart  units  for 
real  time  averages  of  Cy  and  e. 

8)  Acoustic  Radar 

An  Aeroenvi ronment  Model  300  acoustic  radar  was  installed  on  the 
upper  aft  deck  of  the  ship.  The  acoustic  radar  provides  an  information 
cross  section  of  the  lower  atmosphere  by  emitting  a brief  pulse  of  sound 
upward  (f  = 1600  HZ),  listening  to  the  echoes  which  are  reflected  by  tem- 
perature variations  aloft,  and  displaying  the  height  of  the  echoes  on  a 
chart.  This  device  has  been  used  successfully  on  shipboard  in  the  Pacific 
to  measure  the  height  of  the  marine  inversion  up  to  1 km.  No  inversions 
were  observed  during  the  entire  trip. 

IV.  MEAN  DATA 

Temperature  and  Humidity 

Table  B- 1 1-12  contains  the  location  in  the  data  archive  system  of 
mean  values  of  temperatures  and  humidity  obtained  from  the  NPS  equipment 
during  the  trip.  Level  1 is  the  bow  station,  level  2 is  the  radar  mast 
station  and  T$  is  the  sea  "surface"  temperature.  Since  the  sea  tempera- 
ture sensor  was  between  the  hulls  in  a very  turbulent  region,  it  is  as- 
sumed that  it  is  more  of  a "bucket"  temperature  than  an  actual  surface 
temperature. 
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Table  B-II-12 


i 


I • 


NPS  Meteorological  Data  Files 


1) 

V 

Surface  Temperature 

(°C)  : 

NAVPGS.DAT(2) 

2) 

V 

Air  Temperature  level  1 

(°C)  : 

NAVPGS.DAT(3) 

3) 

t2» 

Air  Temperature  level  2 

(8C)  : 

NAVPGS.DAT(4) 

4) 

Hl* 

Relative  Humidity  level  1 

(%)  : 

NAVPGS.DAT (5) 

5) 

Ho, 

Relative  Humidity  level  2 

(%>  : 

NAVPGS.DAT (6) 

V.  TURBULENCE  DATA 
1)  Definitions 

2 

The  temperature  structure  function,  Cy,  is  defined  as 
Cy  = < [T ( x ) - T(x+d)]2  > d-2/3 

where  T(x)  is  the  temperature  at  position  x and  T(x+d)  is  tne  temperature 

2 

at  position  x+d.  In  the  inertial  subrange  of  isotropic  turbulence,  Cy 
is  independent  of  d and  can  be  related  to  the  one  dimensional  power  spec- 
trum of  temperature  fluctuations,  $y,  by 

*y(k)  = .25  C2  k’5/3 

where  k is  the  wave  number.  The  parameter  of  optical  interest  is  the  in- 

o 

dex  of  refraction  structure  function,  CN,  which  is  given  approximately 
by  (at  sea  level), 

C2  = (9.8  x ID'7  »C-Wy 
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The  rate  of  dissipation  of  turbulent  energy,  e,  is  related  to  the  one 
dimensional  power  spectrum  of  velocity  fluctuations,  *u>  by 

♦u(k)  = .5  e2/3k‘5/3 

in  the  inertial  subrange.  In  the  constant  stress  surface  layer,  the  fric- 
tion velocity,  U*,  is  related  to  e by 

eKZ  = U3  fc(Z/L), 


where  K is  von  Karmon's  constant  = .35  and  Z is  the  height  above  the  sur- 
face. In  near  neutral  conditions,  MZ/L)  = 1 and 

U*  = ( cKZ)1/3 

The  friction  velocity  is  related  to  the  mean  velocity  (u)  profile  by 


6u  _ 


U* 


f..  (Z/L), 


62  KZ  'u 


where  fu(Z/L)  is  a stability  correction.  The  eddy  diffusivity,  Km,  is 


given  by 


Km  = KZU*fu"A(Z/L). 


2)  Atlantic  Crossing 

2 

Daily  average  values  of  Cy  and  T are  shown  in  Fig.  B-II-4  for 
the  Atlantic  Ocean  crossing.  This  data  illustrates  the  large  values  of  Cy 
found  at  the  boundary  of  the  warm  waters  of  the  Gulf  Stream  around  May  22. 
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2 

3)  Diurnal  Variation  of 


2 

The  diurnal  variation  of  CT  is  of  interest  to  optical  propaga- 
tion modelers  as  well  as  micrometeorologists.  In  general,  the  marine 
boundary  layer  is  considerably  less  variable  than  its  overland  counter- 
part. We  have  divided  the  data  into  three  separate  geographical  loca- 
tions: the  Atlantic  coast  of  the  U.S.  and  Canada  (May  15-21),  the  warm 
waters  of  the  Mid-Atlantic  (May  22-27)  and  the  Mediterranean  Sea  (May  30- 
June  6).  Figure  B-II-5  is  the  bow  level  data  (Z  = 11  meters).  Figure 
B-II-6  is  the  radar  mast  level  (Z  = 20.5  meters)  and  Figure  B-II-7  shows 
both  levels  for  all  three  regions  combined. 

4)  Turbulence  Data 

2 

Table  B-II-13  is  a compilation  of  the  fluctuation  parameters  Cy 
and  e measured  during  the  cruise.  In  addition,  we  have  the  friction  ve- 
locity U*  which  is  related  to  the  eddy  diffusivity  and  iQ,  the  microscale 
of  turbulence.  The  atmospheric  stability  is  given  in  terms  of  Richardson 
number,  R.. . 

n • _ g ( 60  V/6Z) 

Ri  " “7  T~ 

1 (6U/6Zr 

where  g is  the  gravitational  acceleration  and  is  the  virtual  potential 
temperature. 


Table  B-II-13 


NPS  Data  Files 


1) 

CT12  (10'3oC/m2/3)  level  1 

: NAVPGS.DAT(7) 

2) 

Ct22  (10"3oC/m2/3)  level  2 

: NAVPGS.DAT(8) 

3) 

c1  (10”3m2/sec3)  level  1 

: NAVPGS.DAT(9) 

4) 

e2  (10”3m2/sec3)  level  2 

: NAVPGS.DAT ( 10) 

5) 

U*  (m/sec) 

: NAVPGS.DAT(ll) 

6) 

Richardson's  Number 

: NAVPGS.DAT(12) 
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B-III  METEOROLOGICAL  DATA 

i 

a)  Sea  Surface  Temperature 

Sea  surface  temperatures  were  measured  by  infrared  radiation  thermom- 
eters, fish  towed  by  the  ship  in  the  surface  waters,  and  from  sensors  lo- 
cated at  keel  level  in  a well  on  the  ship.  Each  of  these  measurements 
has  its  own  special  problems  and  meanings  to  investigators  and  files  of 
each  of  these  measurements  are  stored  in  the  data  archive  (see  Section  D). 

However  in  order  to  provide  one  sea  surface  temperature  for  general  use, 
the  following  compromise  is  used.  For  each  hour  an  average  of  all  of  the 
instruments  measuring  sea  surface  temperatures  at  that  particular  time  is 
used  for  the  following  listing  and  in  the  plot  of  SST  vs  time  in  Figure 
B-III-1.  Statistical  analysis  of  the  scatter  of  reading  in  this  data 
show  that  for  the  most  part  the  standard  deviation  of  readings  for  each 
hour  is  less  than  1 degree  C for  most  hours  although  higher  standard  devi- 
ations occur  on  May  IS  and  20,  1977. 

The  frequency  distribution  of  the  sea  surface  measurements  is  shown 
in  Figure  B-III-2  where  the  x axis  represents  temperature  in  degrees  C 
and  the  y axis  is  the  number  of  observations  of  sea  surface  temperature 
which  fall  within  the  limits  of  the  bar  graph  elements.  The  wide  spread 
of  this  data  indicate  that  the  cruise  did  represent  a good  sampling  of 
oceanic  sea  surface  temperature  conditions.  Table  B-III-1  is  a listing 
of  the  sea  surface  temperatures  on  an  hour  by  hour  basis. 
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OF  SEA  SURFACE  TEMPERATURE 
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Table  B-III-1.  Average  Sea  Surface  Temperature 


SEA  SURFACE  TEMPERATURE  (C> 
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Table  B-III-1.  Average  Sea  Surface  Temperature  (cont'd) 
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Table  B-III-1.  Average  Sea  Surface  Temperature  (cont'd) 


I 


b)  Air  Temperature 

Each  hour  of  the  data  from  which  the  air  temperature  plot  shown  in 
Figure  B-III-3  was  made  up  from  is  the  average  value  of  up  to  seven  simul- 
taneous measurements  of  air  temperature.  This  is  our  method  of  providing 
one  standard  air  temperature  value  for  various  calculation  purposes.  The 
data  is  made  up  of  data  taken  on  the  ship's  bow,  on  towers  mounted  above 
the  pilot  house  and  from  thermometers  located  on  the  forward  mast.  These 
instruments  were  operated  by  NRL,  the  Naval  Postgraduate  School  and  Cal- 
span  Corporation.  These  instruments  and  their  locations  are  described  in 
Section  B-II  of  this  report.  The  standard  deviation  of  the  measured  tem- 
perature for  a typical  hour  was  less  than  1 degree  Centigrade  but  on  occa- 
sion standard  deviations  on  the  order  of  4 degrees  C were  observed.  Such 
an  average  value  for  air  temperature  of  course  masks  all  real  effects  of 
sensor  location  and  of  observing  techniques.  All  of  the  individual  spe- 
cialized data  is  of  course  available  in  the  data  archive  system,  the  use 
of  which  is  described  in  Section  0 of  this  report.  Table  B-III-2  is  a 
printed  record  of  each  of  the  average  value  of  air  temperature  data  that 
is  available  for  the  cruise. 

Figure  B-III-4  is  the  frequency  distribution  of  this  data  in  which 
the  y axis  is  the  number  of  observations  of  average  air  temperature  within 
the  indicated  band  and  the  x axis  is  the  air  temperature  in  degrees  C. 

This  curve  like  that  of  the  sea  surface  temperature  in  Figure  B-III-2 
shows  that  indeed  a broad  sampling  of  thermal  conditions  were  encountered 
on  the  EOMET  77  cruise. 
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AUERAGE  AIR  TEMPERATURE 


Fig.  (B-III-3)  Cruise  plot  of  the  average  air  temperature  expres 
in  Celsius  vs.  time  expressed  in  GMT  Julian  days. 


FREQUENCY  OF  AIR  TEMPERATURE  MEASUREMENTS 


Fig.  (B-III-4)  Frequency  of  air  temperature  measurements;  Y 
shows  number  of  observations  made;  X axis  is 
temperature  in  degrees  C. 


AVERAGE  AIR  TEMPERATURE 
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Table  B-III-2.  Average  Air  Temperature 


AUERAGE  AIR  TEMPERATURE  <C> 
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Table  B-III-2.  Average  Air  Temperature  (cont'd) 


AUER AGE  AIR  TEMPERATURE  <C> 
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c)  Air-Sea  Temperature  Difference 


The  data  that  is  present  in  Figure  B-III-5  is  the  calculated  differ- 
ence between  the  air  temperature  and  the  sea  surface  temperature  plotted 
against  time  throughout  the  cruise.  Each  hour  of  the  data  from  which  this 
plot  was  obtained  represents  the  arithmetical  difference  between  the  aver- 
age air  temperature  and  the  average  sea  surface  temperature  for  that  hour. 
The  average  air  temperature  and  the  average  sea  surface  temperature  data 
are  discussed  in  Sections  B-III-a  and  B-IlI-b  of  this  report.  This  plot 
shows  in  a graphic  way  the  areas  where  the  atmospheric  boundary  layer  and 
the  surface  waters  are  not  in  thermal  equilibrium.  The  cold  Nova  Scotian 
coastal  waters  and  the  Gulf  Stream  stand  out  as  significantly  affecting 
the  chemical  stability  of  the  boundary  layer  as  is  seen  by  large  devia- 
tions both  positively  and  negatively  in  this  presentation.  This  same  be- 
havior is  observable  in  the  plot  of  Richardson's  number  in  Section  B-III-g 
of  this  report.  The  printed  version  of  this  calculated  data  for  each  hour 
is  presented  in  Table  B-III-3. 
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AIR  SEA  TEMPERATURE  DIFFERENCE  <C 


Fig.  (B-III-5)  Air-sea  temperature  difference. 

Time  in  Julian  day,  temperature  in  degrees 


AIR  SEA  TEMPERATURE  DIFFERENCE  CC> 
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Table  B-III-3.  Air-Sea  Temperature  Differences 


AIR  SEA  TEMPERATURE  DIFFERENCE  <C> 
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Table  B-III-3.  Air-Sea  Temperature  Differences  (cont’d) 
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Table  B-III-3.  Air-Sea  Temperature  Differences  (cont'd) 


d)  Relative  Humidity 


The  data  shown  in  Figure  B-III-6  represents  average  values  of  rela- 
tive humidity  obtained  by  the  various  methods  and  presented  here  as  the 
representative  relative  humidity  data  for  the  cruise.  This  data  is  calcu- 
lated from  7 different  measurements  that  were  made  throughout  the  cruise. 
The  various  measurements  of  relative  humidity  are  listed  below: 


ORGANIZATION 

METHOD 

DATA  FORMAT  IN  ARCHIVE 

1) 

NRL  (tower) 

Dewpoi nt 

Dewpoint,  temperature 

2) 

Calspan  (bow) 

Lithium  Cloride 

Dewpoint,  temperature 

3) 

Calspan  (mast) 

Lithium  Cloride 

Dewpoint,  temperature 

4) 

/ 

Calspan  (tower) 

Sling  Pschrometer 

Dry  bulb  temperature, 
wet  bulb  temperature 

5) 

NRL  (deck) 

Pschrometer 

Dry  bulb  temperature, 
wet  bulb  temperature 

6) 

NPS  (bow) 

Lithium  Cloride 

Relative  humidity 

7) 

NPS  (mast) 

Lithium  Cloride 

Relative  humidity 

Each  set  of  humidity  measurements  were  converted  into  relative  humid- 
ity values  for  comparison  purposes  and  then  for  each  hour,  up  to  seven 
calculations  of  relative  humidity  are  available.  The  average  value  of 
these  calculations  is  then  presented  as  the  representative  relative  humid- 
ity for  that  particular  hour  on  the  cruise.  The  error  between  the  simul- 
taneous measurements  of  relative  humidity  is  shown  in  Figure  B-III-7  by 
the  frequency  distribution  of  the  standard  deviation  of  these  values  for 
each  hour.  It  shows  that  for  most  cases,  the  calculations  had  a standard 
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1 


r 


deviation  of  about  6 to  7%  relative  humidity.  On  special  occasions  how- 
ever the  standard  deviation  went  up  to  20 % relative  humidity.  The  indi- 
vidual measurements  of  the  particular  instruments  are  available  in  the 
data  archive. 

Figure  B- 1 1 1-8  is  the  frequency  distribution  of  the  relative  humid- 
ity measurements  from  the  EOMET  77  cruise.  It  shows  that  the  relative 
humidities  are  quite  sharply  centered  about  85%.  The  printed  version  of 
the  hourly  representative  relative  humidity  data  is  presented  in  Table 
B-III-4. 
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AVERAGE  RELATIUE  HUMIDITY  MEASUREMENTS  <*) 


Fie  (B-III-6)  Average  relative  humidity  measurements. 
Julian  date  vs.  RH  (7o)  . 


FREQUENCY  OF  RELATIUE  HUMIDITY  <*> 


Fig.  (B-III-8)  Frequency  of  relative  humidity 


AUERAGE  RELATIUE  HUMIDITY  < PERCENT > 
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Table  B-III-4.  Average  Hourly  RH  Measurement  (7o) 


AUERAGE  RELATIUE  HUMIDITY  < PERCENT  > 
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Table  B-III-4.  Average  Hourly  RH  Measurement  (%) 


AUERAGE  RELATIUE  HUMIDITY  < PERCENT) 
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Table  B-III-4.  Average  Hourly  RH  Measurement  (7») 


e)  Real  Wind  (Speed  and  Direction) 


The  calculation  of  the  real  wind  speed  and  direction  from  measure- 
ments aboard  a moving  ship  requires  the  simultaneous  knowledge  of  four 
shipboard  measurements.  They  are  the  relative  wind  speed  as  measured  from 
an  anemometer  located  on  the  moving  ship,  the  direction  of  the  apparent 
wind  relative  to  the  bow  of  the  ship,  the  true  speed  of  the  ship  and  the 
heading  of  the  ship.  While  direct  measurements  of  the  relative  wind  speed 
and  direction  and  the  ships  heading  are  easily  and  accurately  obtained, 
the  speed  of  the  ship  relative  to  a fixed  reference  point  is  difficult  to 
make  because  of  the  problems  of  ocean  currents  and  other  effects  on  the 
ship.  Consequently  the  data  on  ships  speed  used  for  this  cruise  was  ob- 
tained by  means  of  calculating  it  from  accurately  obtained  ships  posi- 
tions. This  approach  assumes  that  the  average  speed  of  the  ship  between 
two  nearby  positions  as  obtained  by  dividing  the  distance  between  the  two 
points  by  the  elapsed  time  required  to  traverse  the  distance  between  the 
points  is  a good  approximation  of  what  the  ship  was  actually  doing  as  it 
traversed  the  track  between  the  two  points.  This  assumption  is  adequate 
if  the  ship  maintained  a constant  heading  and  propeller  rpm  between  the 
two  measurement  points. 

Figure  B-III-9  shows  in  two  plots  the  values  of  the  true  wind  and  its 
direction  as  a function  of  time  for  the  cruise.  The  upper  plot  shows  cal- 
culated values  of  real  wind  direction.  In  this  curve  winds  from  the  north 
have  a value  of  either  U or  36U  degrees  as  measured  along  the  ordinate  of 
the  plot  while  east  winds  have  a value  of  9U  degrees.  The  lower  curve  is 
* plot  of  wind  speed  as  measured  in  meters  per  second  along  the  ordinate 
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and  time  along  the  x axis  is  measured  in  Julian  days.  A periodic  fluctu- 
ation in  this  parameter  is  indicative  of  large  scale  air  mass  structures 
moving  across  the  North  Atlantic.  Figure  B-III-10  is  the  frequency  dis- 
tribution of  the  real  wind  speed  calculations.  Tables  B-III-5  and  B-III-6 
contain  the  hourly  calculations  of  real  wind  speed  and  real  wind  direction 
respectively. 


Fig.  (B-III-10)  Frequency  distribution  of  real  wind  speeds  (m/s). 
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Table  B-III-5.  Average  Wind  Speed  (hourly) 


REAL  WIND  SPEED  CALCULATIONS  < M'S 
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Table  B-III-5.  Average  Wind  Speed  - Hourly  (cont’d) 
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Table  B-III-6.  Average  Wind  Direction  - Hourly  (cont'd) 


Table  B-III-6.  Average  Wind  Direction  - Hourly  (cont'd) 


f ) Friction  Velocit: 


One  micrometeorological  parameter  of  interest  to  boundary  layer  stud- 
ies is  the  friction  velocity.  The  friction  velocity  is  a parameter  which 
is  related  to  the  mean  velocity  profile  end  to  the  eddy  diffusivity  of  mo- 
mentum. The  values  shown  plotted  in  Figure  B-III-11  are  those  calculated 
from  measurements  by  the  NPS  during  the  EOMET  cruise  of  the  USNS  Hayes.  A 
description  of  the  instrumentation  and  some  definition  of  terms  and  rela- 
tionships are  discussed  in  Section  B-II-c  of  this  report.  In  the  figure, 
the  ordinate  scale  of  the  friction  velocity  values  is  meters  per  second. 
The  X axis  is  of  course  the  time  of  the  measurement  expressed  in  Julian 
days.  Table  B-III-7. is  an  hourly  listing  of  the  NPS  calculations  of  fric- 
tion  velocity;'' 
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Table  B-III-7.  Friction  Velocity  Hourly  Listing 


A 

(0 

\ 

E 

v 


> 

»- 
t— i 

O 

O 

-I 

UI 

3 

z 

o 


CJ 

M 

a: 

u. 


-j 

o 

0 

1 
o 
(0 

Ui 

<r 

Z) 

9 

a: 

u 

h* 

W 

o 

Q_ 


<r 

z 


iv 

ft. 

\ 

8 

w> 

ft. 

is- 

\ 

(ft 

0J 

\ 

ID 

ft. 

ft- 

\ 

00 

(VI 

\ 

«n 

ft. 

ft. 

\ moo vo to 

ft-  H rt  ri 

OJ  ... 

\ 6)0 QO 

in 

ft. 

IV 

\ ft-  (ft  (V  ft. 

CJ  ...  . 

NQOiSO 

m 

Iv 

IV 

\ ^M-roro 

in  (Vi  (vi  (vi  (\i 

(VI  ...  . 

\ QOOCD 

in 


a 

® 


{MR  SfeSS 

QOQ  00(5)0 


CVIOO 
h CD 

(SO 


H tH  tH  tH  H r<  H H 
O Q 0 0 0 0 0 0 0 


®co!S?>!S,«:  tfrojn®  gft-N 


ft)  (VI  CVI<VI<VI<VJ  (VICVlMft) 
CDOOOOG)  (90(9  0 


ISOS 

(ft  ft- VO  VD  VO 
(SSMSSMS 


(ft 

KD 

*~4 

N 

H38JE3K3a 

<3 

CD 

OCDSQCEi 

!®'rHOjft),v»-in'a)r--co<fto^<vjft)’M‘invoft.oo(ft®^(vift) 

vh  »-i  th  w ^ ^ ^ v\|  CM  V Ni<M 


-90- 


Table  B-III-7.  Friction  Velocity  Hourly  Listing  (cont'd) 
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Table  B-III-7.  Friction  Velocity  Hourly  Listing  (cont'd) 


g)  Richardson's  Number 


The  dimensionless  Richardson's  number  serves  as  an  estimate  for 
energy  transfer  in  a turbulent  flow  under  non  adiabatic  conditions.  It 
indicates  the  energy  supplied  or  consumed  by  thermal  stratification  as 
compared  with  the  energy  furnished  by  eddy  stresses.  It  is  positive  for 
stable  stratifications  where  the  turbulent  stresses  have  to  work  against 
gravity  and  becomes  negative  in  super  adiabatic  or  unstable  cases. 

Figure  B-I 11-12  is  a cruise  plot  of  the  variation  of  Richardson's 
number  as  a function  of  time/place  during  the  EOMET  77  cruise.  These 
numbers  were  calculated  by  the  NPS  and  details  of  their  measurements  and 
their  calculation  are  found  in  Section  B-II-c  of  this  report.  Table 
B-III-8  is  an  hourly  listing  of  these  data  as  they  were  reported  by  NPS 
during  EOMET  77. 
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RICHARDSON  NUMBER  CALCULATIONS 


Fig,  (B-III-12)  NPS  Richardson  number  calculations 
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Table  B-III-8  Richardson's  Number 
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Table  B-III-8  Richardson's  Number  (cont'd) 
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Table  B-III-8  Richardson's  Number  (cont'd) 


h)  Atmospheric  Radon  Calculation 


222 

Radon  ( Rn)  has  prove"  to  be  an  excellent  indicator  of  the  con- 
tinental nature  of  over  ocean  air  and  air  mass  boundaries.  Radon  is  a 
radioactive  rare  gas  with  a 3.8  day  half  life  which  is  continuously  ema- 
nated from  all  land  areas.  Radon  emanates  from  the  oceans  but  at  a rate 
two  or  three  orders  of  magnitude  smaller  than  from  land  areas,  and  this 
enables  radon  to  be  used  as  an  indicator  of  any  recent  continental  contri- 
bution to  marine  air  masses. 

The  transition  from  pure  "continental"  to  pure  "maritime"  air  pro- 
duces an  order  of  magnitude  change  in  radon  concentration  from  a few  pico- 

-3  -3 

curies  per  cubic  meters  (pCi  m ) or  more  to  a tenth  of  a pCi  m or  less. 
In  many  areas  over  the  North  Atlantic  and  Pacific  Oceans,  two  or  three 
pCi  m radon  is  about  as  low  as  one  can  anticipate  encountering,  while 

_3 

radon  levels  of  a few  tenths  pCi  m can  be  expected  over  some  areas  of 
the  South  Pacific. 

Figure  B-III-13  is  a plot  of  radon  concentration  as  a function  of 
time  as  measured  during  EOMET  77.  The  Y axis  is  the  concentration  in 

_3 

pCi  m while  the  X axis  is  the  time  in  Julian  days.  Table  B-III-9  is 
an  hourly  listing  of  radon  data  as  it  was  obtained  and  recorded. 


Atmospheric  Radon  concentration 
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Table  B-III-9  Atmospheric  Radon  Concentration 
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Table  B-III-9  Atmospheric  Radon  Concentration  (cont’d) 
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Table  B-III-9  Atmospheric  Radon  Concentration  (cont'd) 


r - - - 

i ) Visibility 

Accurate  visibility  measurements  from  on  board  a ship  at  sea  are 
normally  difficult  to  obtain.  Estimates  of  visibility  are  made  routinely 
at  sea  by  observers  who  have  only  the  horizon  on  which  to  base  their 
estimates. 

The  visibility  values  given  here  were  calculated  from  the  well-known 

Koschmieder  formula.  According  to  this  formula  visibility  is  equal  to  the 

constant  3.91  divided  by  the  scattering  coefficient.  The  required  values 

of  the  scattering  coefficient  were  obtained  by  adding  a molecular  scatter- 

-5  -1 

ing  contribution  of  1.65x10  m to  the  values  of  the  aerosol  scatteri<iy 
coefficient  measured  with  the  MRI  model  1567  integrating  nephelometer. 

Figure  B-1II-14  is  the  cruise  plot  of  the  visibility  calculations. 
Figure  B-lll-15  is  the  frequency  distribution  of  these  calculations. 

Table  B-III-10  lists  all  of  the  valid  calculated  visibilities  obtained 
for  the  cruise. 
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UISIBILITY  CALCULATIONS  <KM>  FROM  NRL  SCAT 


Lsibility  - Julian  day  vs. 


THIS  IS  THE  NRL  FREQUENCY  OF  VISIBILITY  MEASURE 


Fig.  (B-III-15)  Frequence  of  visibility  measurements. 


UISIBILITY  CALCULATIONS  <KM>  FROM  NRL  SCAT.  COEF.  MEAS. 
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Table  B-III-10.  Visibility  - Hourly  Listing 
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Table  B-III-10  Visibility  - Hourly  Listing  (cont'd) 


VISIBILITY  CALCULATIONS  <KM>  FROM  NRL  SCAT.  COEF.  MEAS. 
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Table  B-III-10  Visibility  - Hourly  Listing  (cont’d) 


j)  EOMET  77  Meteorological  Soundings 

The  data  in  this  section  were  obtained  by  two  methods.  Dr.  Ralph 
Markson  of  Airborne  Research  Associates  supported  the  cruise  by  making 
aircraft  measurements  of  temperature  and  dewpoint  in  the  vicinity  of  the 
ship  as  it  proceeded  along  the  coast  of  North  America.  A series  of  kite 
balloon  soundings  were  also  made  with  the  NRL  microprocessor  controlled 
boundary  layer  sounding  system  described  in  detail  in  NRL  Report  #8271. 

Figure  B-1II-16  shows  the  location  of  all  of  the  kite  balloon  flights 
which  are  indicated  by  circles  containing  the  flight  numbers  adjacent  to 
the  cruise  track.  Each  flight  consists  of  one  ascent  of  the  kite  balloon 
system  to  a maximum  altitude  where  it  can  no  longer  rise  and  then  its 
descent  again  to  the  surface.  The  time  required  for  one  complete  sounding 
is  approximately  40  minutes.  Every  minute  during  the  flight  the  current 
measurement  of  dry  bulb  temperature,  wet  bulb  temperature,  and  pressure 
altitude  are  digitally  recorded.  Consecutive  flights  in  the  same  general 
area  are  shown  by  adjacent  circles  on  the  map.  In  actuality  the  flight 
numbers  show  the  order  of  the  flights.  There  are  several  missing  flight 
numbers  and  these  correspond  to  real  flight  attempts  but  which  for  one 
reason  or  other  the  flights  were  aborted  without  an  adequate  data 
accumulation. 

Figures  B-I 11-17,  B-III-18  and  B-I 11-19  are  the  aircraft  profiles  ob- 
tained on  15,  17  and  18  May  1977  in  the  vicinity  of  the  ship's  track.  The 
aircraft  used  in  these  measurements  was  a single  engine  turbocharged 
Bellanca  which  has  been  flown  to  as  low  as  3 meters  and  as  high  as  10,000 
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meters.  Temperature  and  dewpoint  data  are  recorded  along  with  other 
parameters  on  an  eight  channel  strip  chart  recorder. 

Figures  B-1II-20  through  B - 1 1 1 - 28  show  plots  of  the  kite  balloon 
soundings  for  all  of  the  successful  flights.  These  plots  show  the  kite 
balloon  data  expressed  in  terms  of  dry  air  temperature  and  its  dewpoint. 

On  the  plots  the  Y axes  show  altitude  of  from  0 to  800  meters  and  the  X 
axis  shows  temperature  of  from  -10  to  25  degrees  centigrade.  On  each  of 
these  flights  the  actual  data  points  have  been  approximated  by  a third 
order  polynomial  least  square  fits  of  the  discrete  flight  data.  Tempera- 
ture is  shown  by  the  solid  curves  and  dewpoint  is  shown  by  the  dashed 
curves. 

A calibration  of  the  kite  balloon  system  instruments  was  obtained  by 
flying  the  instrument  package  at  mast  height  and  comparing  the  measure- 
ments of  the  kite  balloon  package  with  the  precision  ships  instrumenta- 
tion. Such  a calibration  was  made  before  and  after  each  flight  to  watch 
for  a possibility  of  instrumentation  drift  which  might  have  occurred  dur- 
ing the  flight. 

Tables  B-III-11  and  B-III-12  contain  the  coefficients  of  the  fitting 
curves  for  the  air  temperature  data  and  for  the  relative  humidity  data. 

The  actual  data  that  was  used  in  the  calculation  of  the  coefficient  in 
Table  B-III-11  are  the  kite  balloon  dry  air  temperatures  and  the  altitude 
of  the  measurements  together  with  the  simultaneously  measured  shipboard 
mast  height  air  temperature  weighted  by  10  and  the  sea  surface  temperature 
weighted  by  10.  The  least  squares  analysis  produces  coefficients  A,  B,  C 

-109- 


F 


and  D which  show  the  relation  between  dry  air  temperature,  T and  altitude 
z by  the  formula: 

* * 2 * 3 

T(C°)  =A+Bz+Cz  +Dz  where  z is  in  meters. 

This  fitted  curve  is  only  good  for  altitudes  from  the  surface  to  the  max- 
imum altitude  obtained  on  that  particular  flight.  As  with  other  fitted 

polynomials  gross  errors  can  occur  outside  the  area  where  the  original 
points  were  taken. 

The  coefficients  in  Table  B-I 11-12  contain  numbers  which  fit  the 

relative  humidity  at  each  data  point  to  a similar  third  order  polynomial. 

At  each  discrete  data  point,  relative  humidity  in  terms  of  percent  was 
calculated  from  the  dry  and  wee  bulb  temperatures  using  the  usual  pschyo- 
meteric  formulas.  The  resultant  relative  humidity  profile  can  be  repre- 
sented between  the  surface  and  the  maximum  altitude  by  the  formula: 

* * 2 * 3 

Relative  Humidity  (5)  = A + B z + C z + D z (z  is  in  meters). 


-110- 


BELLANCA  N6586V  SOUNDING  1830  Z 17  MAY  1977 


Aircraft  meteorological  sounding  — 17  May 


BELLANCA  N6586V  SOUNDING 
19202  18  MAY  1977 
45  20  N 60  03W 


Fig.  (B— III— 19)  Aircraft  meteorological  sounding  — 18  May 
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Temperature  and  humidity  profile  - kite  balloon 
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Temperature  and  humidity  profile  - kite  balloon 


FLT  16 


-22)  Temperature  and  humidity  profile  - kite  balloon 


-23)  Temperature  and  humidity  profile  - kite  balloon. 


Temperature  and  humidity  profile  - kite  balloon. 


FLT  48 


Temperature  and  humidity  profile  - kite  balloon 


FLT  54 


Temperature  and  humidity  profile  - kite  balloon 


KITE  BALLOON  MEASUREMENT  OF  TEMPERATURE 
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Table  B-III-11  Temperature  Measurement  - Kite  Balloon 


KITE  BALLOON  MEASUREMENT  OF  TEMPERATURE 
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Table  B-III-11  Temperature  Measurement  - Kite  Balloon  (cont'd) 


KITE  BALLOON  MEASUREMENT  OF  RELATIVE  HUMIDITY 
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Table  B-III-12  Relative  Humidity  - Kite  Balloon 


KITE  BALLOON  MEASUREMENT  OF  RELATIUE  HUMIDITY 
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Table  B-III-12  Relative  Humidity  - Kite  Balloon  (cont'd) 


k)  Sea  State  Observations 


Table  B-III-13  lists  the  hourly  visual  observations  of  sea  state  that 
were  obtained  during  the  cruise.  The  coded  values  have  the  meanings  found 
in  Table  B-II-1  of  Section  B-II-a. 
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SEA  STATE  OBSERVATIONS 
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Table  B-III-13  Sea  State  Observations 


SEA  STATE  OBSERVATIONS 
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Table  B-III-13  Sea  state  Observation  (cont'd) 


Table  B-III-13  Sea  State  Observations  (Cont'd) 
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B-IV  AEROSOL  DATA 

One  of  the  most  important  features  of  this  cruise  was  our  ability  to 
measure  various  characteristics  of  the  natural  oceanic  aerosols.  In  this 
section  of  the  report  we  will  discuss  and  present  the  data  from  our  rou- 
tine measurements  of  condensation  nuclei  density,  cloud  condensation  nu- 
clei density,  scattering  coefficients,  the  Royco  and  P.M.S.  optical  aero- 
sol size  analyzers  and  the  Thermo  Systems  Electrical  Aerosol  Analyzer. 

The  ship  itself  is  a prolific  source  of  aerosols  because  of  effluents 

E 

I from  the  stacks  and  galley  therefore  great  care  must  be  taken  to  eliminate 

from  the  data  those  times  in  which  there  is  known  ship  influence.  The 
data  presented  in  this  section  of  the  report  has  been  prescreened  to  elim- 
inate as  much  as  possible  this  source  of  self  influence  on  the  aerosols 
characteristics. 

Even  though  we  experienced  true  winds  from  a westerly  direction  dur- 
ing a large  portion  of  the  cruise,  the  forward  motion  of  the  ship  helped 
to  counteract  this  influence.  The  frequency  distribution  of  the  relative 
wind  direction  is  shown  in  Figure  B-IV-1.  Here  we  see  that  the  most  fre- 
quent relative  wind  direction  was  from  a favorable  quarter.  However  dur- 
ing the  times  in  which  relative  winds  were  from  the  stern  of  the  ship 
(thus  causing  severe  sampling  problems)  the  ship  stopped  its  forward 
motion  and  turned  into  the  wind  from  time  to  time  in  order  to  obtain  accu- 
rate aerosol  measurements. 

| 

L 
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OF  RELATIVE  WIND  DIRECTION 


Frequency  distribution  of  relative  wind  direction 


a)  Condensation  Nuclei  Density  (CN) 


1 


Suspended  solid  and  liquid  particles  or  aerosols  of  such  visible  size 
can  be  present  in  the  atmosphere  in  enormous  numbers.  As  most  of  these 
particles  play  a crucial  role  in  condensation  processes  they  are  usually 
measured  by  a rapid  expansion  of  a chamber  containing  saturated  air.  The 

expansion  causes  an  adiabatic  cooling  to  the  point  where  a high  super- 

saturation causes  the  particles  to  act  as  nuclei  for  condensation  and  the 
droplets  formed  thereby  grow  almost  instantaneously  to  visible  size. 

Figure  B-1V-2  is  a cruise  plot  of  the  condensation  nuclei  density  as 
measured  with  a standard  NRL  Poliak  counter.  This  instrument  was  time 
shared  with  another  instrument  system  and  therefore  was  not  available  for 
continuous  observations.  Figure  B-IV-3  is  a cruise  plot  of  the  condensa- 
tion nuclei  density  as  measured  by  Calspan  with  a Gardner  counter.  It  can 
be  seen  that  these  two  instruments  agreed  quite  well  with  each  other.  The 
frequency  distribution  of  the  CN  data  as  it  was  observed  on  the  E0MET  77 

cruise  is  presented  in  Figure  B-IV-4.  This  figure  contains  the  data  ob- 

tained from  the  Gardner  counter. 

Table  B-IV-1  lists  the  CN  density  for  each  hour  on  the  cruise  where 
a Poliak  counter  measurement  was  taken. 
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JULIAN  DATE  135-13  MAY  1977 


Condensation  nuclei  density-Pollak  counter- Julian 
day  vs . number / cc . 


JULIAN  DATE  133-15  MAY  1977 


Condensation  nuclei  density-Cardner  counter-Julian 
days  vs.  number/cc. 


Fig.  (B-IV-4)  Frequency  distribution  of  CN-Gardner  counter. 


Table  B-IV-1  CN  Density  (Hourly)  - Poliak  Counter 


NRL  CONDENSATION  NUCLEI  OENSITY  <#/CC)  POLLAK 


Table  B-IV-1  CN  Density  (Hourly)  - Poliak  Counter  (Cant'd) 


Table  B-IV-1  CN  Density  (Hourly)  - Poliak  Counter  (Cont'd) 


b)  Cloud  Condensation  Nuclei 


It  is  generally  agreed  that  only  a small  fraction  of  the  condensation 
nuclei  in  the  air  are  activated  in  natural  cloud  or  fog  condensation  pro- 
cesses where  only  very  small  supersaturations  exist.  The  subset  of  atmo- 
spheric condensation  nuclei  which  can  be  activated  by  natural  conditions 
are  called  cloud  condensation  nuclei.  The  data  presented  in  this  report 
were  obtained  by  the  NRL  thermal  gradient  diffusion  chamber.  Figure 
B-IV-5  is  a cruise  plot  of  the  CCN  density  obtained  with  the  NRL  CCN  cham- 
ber for  a temperature  differential  across  the  thermal  gradient  plates  of 
4.3°C.  Table  B-IV-2  is  a listing  of  the  hourly  values  of  this  parameter 
during  the  cruise  as  they  were  obtained. 
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JULIAN  DATE  133-15  MAY  1977 


Fig.  (B-IV-5)  CCN  density  using  NRL  CCN  chamber  with  4.3°C 

temperature  differential;  number/cc  vs.  Julian  date 


Table  B-TV-2 


CCN  Density  (hourly  listing)  - from  NRL  CCN  chamber 
having  4,3®C  Temperature  Differential 


CLOUD  CONDENSATION  NUCLEI  DENSITY  < u/CC  4.3  DEG.  > 


HOUR  5/15/77  5/16/77  5/17/77  5/16/77  5 / 13/77  5/20/77  5/21/77  5/22/77 


£ 274 

3 4020  5469  932  311 
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c)  Scattering  Coefficient 


The  scattering  coefficient  was  measured  simultaneously  by  two  MRI 

integrating  nephelometers.  A high  degree  of  reliability  is  attributed  to 

these  measurements  because  of  the  good  correlation  between  these  aerosol 

measurements.  Figure  B-IV-6  is  a scatterplot  of  simultaneous  measurements 

of  the  scattering  coefficient  from  the  two  instruments.  Values  along  each 

4 

axis  are  the  measured  scattering  coefficient  multiplied  by  10  and  ex- 
pressed in  (meters-1).  The  x axis  represents  the  Calspan  MRI  model  2050 
data  and  the  y axis  is  the  NRL-MRI  model  1567  data. 

Figure  B-IV-7  is  the  cruise  plot  of  the  NRL  measured  scattering  co- 

4 _i 

efficient  multiplied  by  10  and  expressed  in  meters  along  the  y axis 
and  the  Julian  time  is  expressed  along  the  x axis. 

4 

The  frequency  distribution  of  the  NRL  scattering  coefficient  (x  1- 
(m-1))  is  shown  in  Figure  B-IV-8.  Table  B-IV-3  contains  the  hourly  values 
of  the  NRL  scattering  coefficient  obtained  during  the  cruise. 


SCATTERING  COEF.  <*  10-4>  EOMET  77 
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Table  B-IV-3.  NRL  Scattering  Coefficient  - Hourly  Log 
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Table  B-IV-3  NRL  Scattering  Coefficient  - Hourly  Log  (Cont'd) 
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Table  B-IV-3  NRL  Scattering  Coefficient  - Hourly  leg  (Cont'd) 


d)  Royco  Measurements 


There  were  two  Royco  aerosol  spectrometers  operating  in  a continuous 
monitoring  mode  throughout  the  cruise.  The  channel  sizes  of  the  two  in- 
struments are  slightly  different;  nevertheless  the  two  instruments  ap- 
peared to  track  very  well.  Both  sets  of  data  are  of  course  in  the  data 
bank  and  may  be  addressed  at  will  by  means  of  the  archiving  system 
described  in  Section  D. 

Figure  B-IV-9  is  a cumulative  plot  obtained  with  the  archive  system 
of  the  NRL  Royco  data  for  0500  GMT  on  18  May  1977  showing  the  measured 
number  of  particles/cc  greater  than  the  indicated  size.  Figure  B-IV-10 
is  the  plot  of  the  Cal  span  Royco  data  for  the  same  hour  showing  satisfac- 
tory agreement  for  experimental  field  data. 

Figures  B-IV-11  through  B-IV-15  are  cruise  plots  expressing  the  num- 
ber of  particles  per  cc  obtained  in  each  of  the  five  channels  of  the  NRL 
Royco  as  functions  of  time  throughout  the  cruise. 

Table  B-IV-4  is  a listing  of  the  NRL  Royco  data  where  the  numbers 
under  the  heading: 

CHI  refers  to  log1Q  [#/cc]  of  particles  with  0.45  < d _<  0.6  microns 

CH2  refers  to  log1Q  [#/cc]  of  particles  with  0.6  < d ^ 1.5  microns 

CH3  refers  to  log1Q  [#/cc]  of  particles  with  1.5  < d <_  2.5  microns 

CH4  refers  to  log10  [#/cc]  of  particles  with  2.5  < d < 4.0  microns 

and  CH5  refers  to  log1Q  [#/cc]  of  particles  with  d > 4.0  microns. 
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F 


Each  page  of  the  table  contains  two  columns  with  each  sample  referred  to 
by  Julian  day  and  GMT  hour  within  the  day. 

Table  B-IV-5  contains  Calspan's  Royco  data  that  is  described  in  Sec- 
tion B-II-b  of  this  report. 
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plot  of  NRL  Royco 
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Fig.  (B-IV-10)  Calspan  Royco  data 
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data  channel  4 


NRL  ROYCO  DATA  < DCMICRONSUM . 0 X.  #/CC > 


NRL  ROYCO  SIZE  DISTRIBUTIONS  EOMET  77 
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Table  B-IV-4.  NRL  Royco  Size  Distributions  - Hourly  Log 


NRL  ROVCO  SIZE  DISTRIBUTIONS  EOMET  77 
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Table  B-IV-4  NRL  Royco  Size  Distributions  - Hourly  log  (Cont'd) 
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Table  B-IV-4  NRL  Rayco  Size  Distributions  - Hourly  Log  (Cant’d) 
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CALSPAN  ROYCO  SIZE  DISTRIBUTIONS  EOMET  77 
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Table  B-IV-5.  Calspan  Royco  Size  Distribution  - Hourly  Log 
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Table  B-IV-5  Calspan  Royco  Size  Distributions  - Hourly  Log  (Cont'd) 


CALSPAN  ROYCO  SIZE  DISTRIBUTIONS  EOMET  77 
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Table  B-IV-5  Calspan  Boyco  Size  Distributions  - Hourly  Log  (Cont'd) 
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Table  B-IV-5  Calspan  Roycx)  Size  Distributions  - Hourly  Log  (Cont'd) 


CALSPAN  ROYCO  SIZE  DISTRIBUTIONS  EOMET  77 
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Table  B-IV-5  Calspan  Royco  Size  Distributions  - Hourly  Log  (Cont'd) 


Table  B-IV-5  Calspan  Royco  Size  Distributions  - Hourly  Log  (Cont'd) 


CALSPAN  ROYCO  SIZE  DISTRIBUTIONS  EOMET  77 


e)  Electrical  Aerosol  Analyzer  Data 


The  data  in  this  section  was  obtained  with  Calspan's  Thermo  Systems 
Model  3030  Electrical  Aerosol  Analyzer,  EAA,  with  its  sensor  mounted  ap- 
proximately 18  meters  above  the  sea  surface.  This  device  provides  the 
aerosol  size  distribution  for  sizes  between  0.003  micrometers  and  1.0  mi- 
crometers diameter.  Data  in  this  section  represents  the  actual  size  dis- 
tributions in  cumulative  form.  They  were  obtained  after  computer  process- 
ing at  Cal  span  of  the  raw  EAA  data.  The  archived  data  has  been  sorted  to 
eliminate  periods  of  adverse  relative  wind  direction.  Figure  B-IV-16  is 
a typical  plot  of  the  results  of  the  EAA  where  the  particle  diameter  in 
micrometers  is  plotted  along  the  x axis  and  the  cumulative  number  concen- 

_3 

trations  in  (cm  ) is  plotted  along  the  y axis. 

Table  B-IV-7  is  a listing  of  the  archived  data  where  for  each  day  and 
hour  the  numbers  listed  are  the  logarithms  to  the  base  of  10  of  the  cumu- 
lative  number  densities  (cm  ).  The  corresponding  channel  sizes  are 
listed  in  Table  B-IV-6. 


Table  B-IV-6 

Channel  Number  vs  Aerosol  Size 
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WHITBY  DATA  FOR  18  MAY  1977,  580  GMT 


Fig.  (B-IV-16)  Whitby  data 
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Table  B-IV-7  Whitby  Aerosol  Size  Distributions-Hourly  Log  (cont'd) 
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Table  B-IV-7  Whitby  Aerosol  Size  Distributions-Hourly  Log  (cont'd) 
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Table  B-IV-7  Whitby  Aerosol  Size  Distributions-Hourly  Log  (cont'd) 
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Table  B-IV-7  Whitby  Aerosol  Size  Distributions-Hourly  Log  (cont'd) 


Particulate  Aerosol  Size  Spectra 
as  Measured  by  the  PMS 
Axially  Scattering  Spectrometer  Probe 

Richard  K.  Jeck 
Code  8323 

Naval  Research  Laboratory 
I.  Measurement  Conditions  and  Procedure 

3 

Aerosol  particle  size  spectra  (number  of  particles  per  cm  per  micron 
radius)  over  the  range  0.7  to  46  um  diameter  were  obtained  from  measure- 
ments aboard  the  USNS  Hayes.  The  ambient  aerosol  particles  were  monitored 
from  an  elevated  platform  about  5m  above  the  forward,  STBD  corner  of  the 
FJying  Bridge  (04  level).  The  platform  was  about  22  rn  above  sea  level  and 
directly  above  the  leading  edge  of  the  superstructure.  This  location 
placed  the  probe  in  air  that  was  free  from  turbulence  generated  by  the 
ship's  superstructure,  for  relative  wind  directions  from  about  70°  on  the 
PORT' side  to  90°  or  more  on  the  STBD  side.  For  nearly  forward  relative 
winds  the  air  intercepted  by  the  probe  originated  at  a level  about  4 m 
below  the  height  of  the  probe  and  is  gently  deflected  up  to  the  level  of 
the  probe  as  it  nears  the  forward  superstructure.  The  air  trajectory  is 
inclined  about  10°  upward  at  the  location  of  the  probe.  This  location 
also  keeps  the  probe  above  the  levels  where  the  air  is  heated  by  proximity 
to  the  ship  surface  in  daytime  hours.  Measurements  showed  that  on  sunny 
days  the  air  temperature  was  as  much  as  2°C  wanner  at  a level  3 m below 
the  elevated  platform  than  at  the  platform  itself. 
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The  probe  referred  to  was  an  Axially  Scattering  Spectrometer  Probe 
(ASSP) , a single  particle  optical  counter  manufactured  by  Particle  Measur- 
ing Systems  (PMS),  Inc.  The  ASSP  has  four  selectable  size  ranges  covering 
the  interval  0.7  to  46  um  diameter  and  each  range  is  subdivided  into  fif- 
teen particle  size  categories,  or  channels,  as  listed  in  Table  B-IV-8. 

The  particle  counts  were  allowed  to  accumulate  in  the  size  channels  during 
successive,  one  minute  sampling  intervals  and  the  counts  were  then  re- 
corded at  the  end  of  each  minute.  The  probe  was  set  to  alternate  between 
Ranges  4 and  3,  each  range  being  in  effect  for  a full,  one  minute  sampling 
interval . 

The  PMS  system  was  not  operational  until  May  18,  the  third  day  under- 
way. After  that,  the  ASSP  was  on-line  continuously  except  for  periods  of 
maintenance  and  for  situations  where  the  wind  was  from  the  aft  direction 
of  the  ship.  The  off-line  time  was  about  30%  of  the  cruise,  mostly  due  to 
following  winds  and  the  resulting  poor  sampling  conditions. 

II.  Explanation  of  Data  Tables 

The  data  table,  B-IV-9,  contains  values  of  dn/dr  (number  of  particles 

3 

per  cm  per  micron  radius)  computed  from  particle  counts  in  thirteen  size 
channels  whose  threshold  radii  are  listed  in  the  left  hand  column  of  each 
page. 

The  data  came  from  the  15  size  channels  of  Range  4 and  the  15  from 
Range  3.  Since  the  first  eight  channels  of  Range  3 are  redundant  with 
Range  4,  only  the  results  from  channels  9 through  14  of  Range  3 have  been 
listed.  These  latter  results  are  the  last  six  entries  in  each  column. 
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The  first  seven  entries  are  from  alternate  size  channels  arbitrarily 
selected  from  the  15  available  channels  for  Range  4. 


For  the  purposes  of  this  data  report,  the  particle  counts  were  arbi- 
trarily averaged  over  one  hour  intervals  in  order  to  keep  the  page  count 
down  to  a tractable  number.  The  Julian  day  and  the  hour  at  which  the  one 
hour  average  begins  are  given  at  the  top  of  each  column.  The  entries  for 
dn/dr  are  given  in  computerized  scientific  notation  where,  for  example, 

_3 

. 77E-03  stands  for  0.77x10  . Some  data  exist  that  were  not  included  in 
this  report.  They  correspond  to  periods  when  valid  measurements  were  ob- 
tained for  only  portions  of  an  hour. 

III.  Additional  Information 

1.  Assumptions  Involved  in  the  Probe  Calibration. 

The  threshold  radii  assigned  to  the  various  size  channels  corre- 
sponds to  the  expected  probe  response  function  for  water  droplets  of  re- 
fractive index  n = 1.33  (Jeck,  1978).  It  is  assumed  that  most  of  the 
maritime  aerosol  particles  in  the  size  range  to  which  the  ASSP  is  sensi- 
tive exist  as  solution  droplets  of  sea  salt  or  other  soluble  material  at 
the  humidities  encountered  on  this  cruise.  It  is  further  assumed  that 
these  droplets  are  sufficiently  dilute  that  the  actual  refractive  index  is 
the  same  as  that  for  pure  water.  In  any  case,  realistic  deviations  from 
n = 1.33  are  likely  to  be  small  enough  that  an  error  of  no  more  than  + one 
size  channel  will  be  made  when  the  response  curve  for  water  droplets  is 
used  to  determine  the  size  of  these  particles.  This  amount  of  error  is 
within  the  manufacturer' s specifications  anyway. 
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2.  Method  Used  to  Compute  dn/dr 


| 

; 

i 


For  particle  spectrometers  with  sufficiently  narrow  size  channels, 
approximate  values  of  dn/dr  are  commonly  arrived  at  for  each  size  channel 

3 

by  dividing  the  particle  count  by  the  product  of  the  sampling  volume  (cm  ) 
and  the  channel  width  ( un) . The  resultant  value  is  then  taken  to  be  rep- 
resentative of  the  population  of  particles  at  some  discrete  radius  within 
the  size  channel.  This  radius  may  be  chosen  as  the  lower  threshold,  the 
center,  or  any  other  location  within  the  applicable  size  channel,  depend- 
ing on  the  slope  or  some  other  weighting  function  of  the  indicated  size 
distribution.  The  resulting  size  spectrum  consists  of  a plot  or  table  of 
discrete  values  of  dn/dr  each  associated  with  a selected,  fixed  value  of 
particle  radius.  Values  of  dn/dr  at  intermediate  radii  are  then  inferred 
by  the  use  of  some  interpolation  scheme. 

With  the  application  in  mind  of  using  the  particle  data  for  deriving 
scattering  coefficients  or  other  optical  properties  which  involve  integra- 
tions over  the  particle  size  spectrum,  an  alternate  method  has  been  devel- 
oped for  computing  dn/dr  from  the  raw  particle  counts  (Jeck  and  Ruhnke). 

In  this  method,  continuous  functions  of  the  familiar  "power  law"  form 

dn  _ Cr-(p+l)  (1) 

dr 

are  assumed  to  describe  the  size  distribution  within  each  size  channel. 

The  recorded  particle  count,  N,  in  a given  size  channel  is  then  used  to 
evaluate  the  constant,  C,  for  that  channel  from  the  relationship 


(2) 


i 


-188- 


where  6 is  evaluated  separately  from  an  iterative  scheme,  and  rfl  and  rb 
are  the  lower  and  upper  thresholds,  respectively,  for  the  channel. 

The  advantage  of  this  method  is  that  dn/dr  is  represented  within  each 
size  channel  by  a continuous,  analytic  function  rather  than  by  a point 
value  assigned  to  some  arbitrary  radius  within  the  channel.  These  func- 
tions, one  for  each  size  channel,  can  then  be  integrated  analytically  in 
a piecewise  fashion  along  with  the  scattering  efficiency  factor  or  other 
weighting  function  to  compute  the  desired  optical  property  or  other  result. 

For  the  purposes  of  this  data  report,  however,  it  was  decided  to  re- 
tain the  more  familiar  data  format  where  point  values  of  dn/dr  at  discrete 
radii  are  presented  rather  than  values  of  C and  6 for  each  channel.  The 
tabled  values  of  dn/dr  have  still  been  computed  from  Eq.  (1),  as  described 
above,  for  tne  lower  threshold  radius  (values  in  left  hand  column  of  each 
page)  of  the  size  channels  selected  for  use.  Since  no  values  for  C and  B 
are  included,  approximate  values  of  dn/dr  at  intermediate  radii  may  be  ob- 
tained by  interpolation. 


3.  Comparison  with  Other  Measurements. 

It  should  be  mentioned  that  a consistent  discrepancy  was  found 
between  the  values  of  dn/dr  indicated  by  the  ASSP  and  those  indicated  by 
two  other  PMS  probes  employed  by  Dr.  Gary  Trusty  on  this  cruise.  See  the 
following  section,  B-IV-g.  These  Classical  Scattering  Aerosol  Spectrom- 
eters (CSAS)  cover  the  range  1-15  wn  in  radius  and  thus  are  sensitive  to 
the  same  sizes  as  Ranges  4 and  3 of  the  ASSP.  The  CSAS  probes  generally 
indicated  values  of  dn/dr  that  were  smaller  by  an  order  of  magnitude  or 
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more  than  those  indicated  by  the  ASSP  for  the  same  particle  sizes.  The 
source  of  this  discrepancy  is  still  unclear  at  this  writing,  but  the  CSAS 
values  appear  to  be  unusually  low  compared  to  previous,  at-sea  measure- 
ments by  other  investigators. 
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Table  B-IV-8.  Calibration  Data  for  PMS  Probes 


Axially  Scattering  Spectrometer  Probe  (ASSP) 

Channel  Range  4 (0.37  - 4.4y)  Range  2 (0.4  - 7.7y) 


Number 

Radius  at 
ctr.  of  chan. 

Range 

chan. 

! within 
limits 

Radius  at 
ctr.  of  cb»n. 

Range 

chan. 

! within 
limits 

1 

0.42y 

0.37 

- 0.45y 

0 . 55y 

0.42 

- 0. 70y 

2 

0.52 

0.45 

- 0.60 

0.88 

0.70 

- 1.05 

3 

0.72 

0.60 

- 0.82 

1.25 

1.05 

- 1.45 

4 

0.95 

0.82 

- 1.05 

1.65 

1.45 

- 1.85 

5 

1.18 

1.05 

- 1.30 

2.10 

1.85 

- 2.38 

6 

1.40 

1.30 

- 1.52 

2.65 

2.38 

- 2.92 

7 

1.65 

1.52 

- 1.75 

3.20 

2.92 

- 3.48 

8 

1.88 

1.75 

- 2.05 

3.75 

3.48 

- 4.02 

9 

2.22 

2.05 

- 2.38 

4.30 

4.02 

- 4.58 

10 

2.55 

2.38 

- 2.72 

4.80 

4.58 

- 5.10 

11 

2.90 

2.72 

- 3.08 

5.35 

5.10 

- 5.60 

12 

3.22 

3.08 

- 3.40 

5.85 

5.60 

- 6.15 

13 

3.58 

3.40 

- 3.75 

6.40 

6.15 

- 6.70 

14 

3.90 

3.75 

- 4.08 

7.00 

6.70 

- 7.25 

15 

4.25 

4.08 

- 4.42 

7.50 

7.25 

- 7.75 

Range  3 

(0.7  - 15. 5p) 

Range  1 

(1.05  - 23. 2y) 

1 

1. 05y 

0.7  - 1 . 4y 

1 . 65y 

1.0  - 2 . 2y 

2 

1.8 

1.4  - 2.2 

3.0 

2.2  - 3.7 

3 

2.8 

2.2  - 3.3 

4.5 

3.7  - 5.2 

4 

3.8 

3.3  - 4.3 

6.0 

5.2  - 6.7 

5 

4.85 

4.3  - 5.3 

7.5 

6.7  - 8.3 

6 

5.9 

5.3  - 6.4 

9.0 

8.3  - 9.8 

7 

7.0 

6.4  - 7.5 

10.5 

9.8  -11.3 

8 

8.0 

7.5  - 8.5 

12.0 

11.3  -12.7 

9 

9.0 

8.5  - 9.5 

13.5 

12.7  -14.2 

10 

10.0 

9.5  -10.5 

15.0 

14.2  -15.7 

11 

11.0 

10.5  -11.5 

16.5 

15.7  -17.2 

12 

12.0 

11.5  -12.5 

18.0 

17.2  -18.7 

13 

13.0 

12.5  -13.5 

19.5 

18.7  -20.2 

14 

14.0 

13.5  -14.5 

21.0 

20.2  -21.7 

15 

15.0 

14.5  -15.5 

22.5 

21.7  -23.2 
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AEROSOL  MEASUREMENTS  FOR  OPTICAL  EXTINCTION  PREDICTION 
G.  L.  Trusty  and  T.  H.  Cosden 
Code  5568 

Naval  Research  Laboratory 
Washington,  D.C.  20375 

This  section  presents  the  results  of  experiments  carried  out  by  the 
Optical  Sciences  Division  of  NRL.  The  experiments  were  designed  to  mea- 
sure aerosol  size  distributions  to  be  used  for  calculating  extinction  co- 
efficients as  a function  of  wavelength. 

Particle  Measuring  Systems’  (PMS)  particle  spectrometers  measured  the 
aerosol  size  distributions.  Two  of  the  instruments  sat  just  forward  of 
the  bridge  at  a height  of  15  meters  from  the  ocean  surface.  One  was  an 
Active  Scattering  Aerosol  Spectrometer  Probe  (ASASP),  which  measures  par- 
ticles with  radii  in  the  size  range  of  0.1  ym  to  2.0  ym.  The  second  was 
a High-Volume  Classical  Scattering  Aerosol  Spectrometer  Probe  (CSASP), 
which  covers  a size  range  of  1.0  ym  to  15.0  ym. 

As  examples  of  marine  atmosphere  aerosol  data,  we  have  listed  the 

twenty-minute  averages  of  measurements  that  we  made  during  the  times  the 

ship  was  stopped  with  the  bow  into  the  wind.  Table  B-IV-10  gives  the  par- 

-3  -1 

tide  density  distribution  dN/dR  (cm  ym  ) as  a function  of  the  radius 
of  the  probe  bin  centers.  For  Probe  1 (ASASP)  we  give  results  from  only 
the  first  seven  bins.  Due  to  a double-valued  sensitivity  in  the  detection 
response,  the  data  for  the  larger  size  ranges  of  that  probe  have  proven  to 
be  inconsistent  in  many  instances.  For  the  purpose  of  calculating  extinc- 
tion coefficients,  we  fit  a line  between  the  value  for  the  seventh  bin  of 
Probe  1 and  the  first  bin  of  Probe  2 (CSASP). 
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Table  B-IV-11  gives  the  meteorological  parameters  and,  for  four  wave- 
lengths, calculated  extinction  coefficients.  {The  extinction  calculations 
do  not  include  molecular  absorption.) 

We  have  listed  the  particle  probe  bin  edges  in  Table  B-IV-12  to  aid 
those  who  may  wish  to  put  the  aerosol  data  into  a form  different  from  the 
one  we  have  chosen. 
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Table  B-IV-10  Twenty-minute  averages  of  dN/dR  (cm  ym“  ) versus  radius  (ym) • 
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Table  B-IV-11  Twenty-minute  averages  of  measured  and  calculated  parameter 


Marine  Observations  with  the  NRL  Mobility  Size  Spectrometer 


1 


W.  A.  Hoppel 

Naval  Research  Laboratory 
Washington,  D.C.  20375 

The  mobility  analyzer  recently  developed  at  NRL  was  installed  on  the 
USNS  Hayes  for  the  1977  EOMET  cruise.  The  mobility  analyzer  could  func- 
tion either  as  a size  spectrometer  or  to  strip  off  and  transmit  particles 
with  a nearly  monodisperse  size  to  the  TGDCC. 

When  operated  as  a size  spectrometer  it  could  measure  the  size  dis- 
tribution in  the  size  range  from  0.006  to  0.3  ym  which  is  well  below  the 
lower  detection  limit  of  optical  counters.  It  is  this  size  range  which  is 
most  important  in  the  formation  of  fogs  and  clouds. 

In  the  mobility  analyzer  the  particles  are  brought  to  charge  equilib- 
rium and  then  the  mobility  distribution  of  the  charge  fraction  is  measured. 
The  mobility  distribution  is  then  converted  to  a size  distribution  using 
the  equilibrium  (Boltzmann)  charge  distribution.  The  aerosol  detection 
system  used  with  the  mobility  analyzer  was  a Poliak  expansion  cloud  cham- 
ber. This  imposed  a rather  severe  lower  limit  on  the  sensitivity  of  the 
analyzer.  The  mobility  resolution  of  the  analyzer  is  determined  by  the 
ratio  of  the  flow  rate  of  the  sheath  air  to  that  of  the  sample  flow; 
whereas,  the  signal  is  inversely  proportional  to  that  ratio.  In  the  ma- 
rine environment  the  aerosol  count  is  so  low  that  the  desired  resolution 
was  sacrificed  to  increase  the  signal.  Even  then  in  the  very  clean  marine 
environment  the  number  density  was  too  low  to  give  reliable  results.  As 

_3 

a general  rule  the  total  aerosol  count  had  to  be  above  1000  cm  to  obtain 


-211- 


usuable  results.  (A  new  aerosol  detection  system  which  will  greatly  im- 
prove the  sensitivity  is  under  construction.) 


Figure  C-I-l  shows  three  size  distributions  obtained  with  the  mobil- 
ity analyzer.  The  Norfolk,  VA  Harbor  size  distribution  shows  many  parti- 
cles in  the  smaller  size  range  typical  of  urban  air.  One  hundred  miles 
off  the  East  Coast  the  very  small  particles  had  disappeared  presumably  by 
coagulation  and  the  peak  in  the  size  distribution  was  at  a radius  greater 
than  2 x 10"6  cm.  As  will  be  evident  from  the  CCN  data  given  later,  all 
of  these  aerosols  were  active  as  CCN  at  supersaturations  of  \%  or  less. 
The  spectrum  labeled  Mi d-Atl antic  is  the  average  of  several  spectra  taken 
in  the  Mid-Atlantic  where  the  number  of  particles  was  so  low  that  any 
single  spectrum  was  judged  to  be  unreliable. 


The  dotted  lines  represent  concurrent  data  taken  with  the  PMS  active 
scattering  probe. 

The  mobility  analyzer  can  also  be  used  to  strip  off  and  transmit  an 
aerosol  composed  of  uniform  particle  size.  This  mode  of  operation  was 
used  to  study  the  relationship  between  aerosol  size  and  critical  super- 
saturation necessary  for  nucleation  and  droplet  growth.  The  relationship 
between  critical  supersaturation  S£  and  dry  radius  rQ  is  given  by 

_A^\1/2  (1) 

B r 3 / 

where  A is  a constant  containing  the  surface  tension  and  Bq  is  given  by 
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the  density  of  the  dry  particle,  and  pw  the  density  of  water. 


Because  of  the  small  aerosol  concentrations  encountered  on  the  ship- 
board measurements  the  size  interval  of  the  transmitted  aerosol  was  in- 
creased to  obtain  a detectable  number  of  nucleated  droplets  in  the  TGOCC. 
This  of  course  decreased  the  resolution  with  which  the  relation  between 
dry  size  and  critical  supersaturation  could  be  measured.  Knowing  the  re- 
lationship between  dry  size  r and  critical  supersaturation  S makes  it 

U v 

possible  to  determine  BQ  in  Eq.  (1).  Values  of  BQ  for  three  periods  were 
determined.  These  were  taken  in  Norfolk,  VA  harbor,  100  miles  off  the  New 
Jersey  coast,  and  in  the  Mediterranean  and  had  the  values  of  0.28  ± 0.04, 
0.40  ± 0.09,  and  0.143  respectively.  The  first  two  were  taken  immediately 
after  the  corresponding  size  distributions  shown  in  Figure  C-I-l  and  the 
stated  uncertainty  in  these  values  is  a result  of  the  large  size  interval 
used.  The  treatment  of  the  Mediterranean  data  is  discussed  in  more  detail 
in  Hoppel,  Gerber  and  Wojciechowski  (1977)  and  is  the  average  of  six  runs 
during  a period  when  the  size  distribution  remained  nearly  constant.  The 
transmitted  size  was  chosen  to  be  near  the  peak  in  the  distribution  and 
hence  were  not  the  same  in  each  case.  The  transmitted  size  intervals  were 
0.022-0.0318  urn,  0.0318-0.0468  wn,  and  0.0468-0.070  urn  for  the  Norfolk, 

New  Jersey  Coast  and  Mediterranean  data,  respectively.  These  values  of 
Bq  should  not  be  taken  to  be  constants  at  any  given  location  but  variables 
which  also  depend  upon  air  mass  history. 
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The  equation  which  gives  the  relationship  between  the  saturation 


ratio  S and  equilibrium  size  can  be  written  in  terms  of  BQ  in  the  follow- 
ing form: 


It  is  clear  from  Eq.  (3)  that  in  addition  to  Bq,  e and  pd  must  also 
be  known  to  obtain  S as  a function  of  r.  However  if  p'  is  also  expressed 
in  terms  of  e,  pw,  and  pd  and  if  the  sensitivity  of  Eq.  (3)  to  various 
values  of  e and  pd  is  studied  numerically,  it  is  found  that  Eq.  (3)  is 
very  insensitive  to  e and  pd  over  ranges  of  these  variables  expected  in 
atmosphere  aerosols.  This  means  that  using  an  experimental  value  of  BQ  in 
Eq.  (3)  together  with  any  reasonable  value  of  e and  pd  will  give  the  rela- 
tive humidity  dependence  of  the  equilibrium  size.  This  dependence  for  the 
three  different  values  of  Bq  is  shown  in  Figure  C-I-2.  This  relationship 
should  hold  for  relative  humidities  above  the  deliquescent  point  of  the 
particular  solute  in  the  aerosol  and  thus  the  curves  in  Figure  C-I-2  are 
only  valid  for  R.H.  above  the  deliquescent  point. 
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Fig.  (C-I-l)  Size  distribution  of  aerosol  at  Norfolk,  NJ  coast, 
mid-Atlantic,  using  NRL  mobility  analyzer. 


L 
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RELATIVE  HUMIDITY  (%) 

I (7o)  at  Norfolk,  NJ  coast,  Mediterranean 


SHIPBOARD  MEASUREMENTS  WITH  A PORTABLE  TRANSMISSOMETER 


H.  Gerber  and  R.  Stilling 

Naval  Research  Laboratory 
Washington,  D.C.  20375 

1.  INTRODUCTION 

The  purpose  of  our  participation  in  the  1977  EOMET  cruise  of  the  USNS 
Hayes  was  to  test  the  feasibility  of  using  a portable  transmissometer  on- 
board a ship  at  sea.  If  proved  successful  this  transmissometer  can  become 
a valuable  new  tool  for  obtaining  maritime  transmittance  measurements 
which  are  presently  very  scarce.  This  lack  of  data  is  due  to  the  fact 
that  existing  transmi ssometers  require  a baseline  of  about  1 km  which  is 
obviously  incompatible  with  ships. 

Transmi ttances  witn  a resolution  of  0.1$  over  a 20  meter  folded  path 
have  been  measured  in  laboratory  aerosol  with  the  new  transmissometer. 

This  extraordinary  improvement  over  the  accuracy  of  long-baseline  instru- 
ments is  due  to  several  unique  features  which  will  be  discussed.  The  0.1% 
accuracy  reflects  a laboratory  environment  which  is  amenable  to  precise 
optical  measurements.  Since  such  an  environment  is  improbable  on  a ship, 
this  accuracy  will  deteriorate.  It  was  the  goal  of  our  participation  on 
the  1977  cruise  to  discover  the  influence  of  wave-induced  ship's  motion, 
vibration  due  to  the  ship's  machinery,  and  temperature  variations  on  the 
accuracy  of  the  new  transmissometer.  The  following  includes  a brief  dis- 
cussion of  the  new  technique,  a list  of  the  collected  data,  and  a summary 
of  the  results.  A detailed  account  of  the  new  system  and  its  performance 
on  the  cruise  are  given  elsewhere  (Gerber,  1978). 
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2.  INSTRUMENTATION 


A schematic  of  the  transmission  cell  is  shown  in  Figure  C-II-1.  A 
HeNe  laser  beam  is  passed  through  a beam  splitter.  One  portion  of  the 
beam  steps  in  a sawtooth  fashion  through  the  2 meter  cell  for  a total 
pathlength  of  about  20  meters,  upon  leaving  the  cell  it  passes  through  a 
chopper  before  its  intensity  is  measured  by  a sensor.  The  second  portion 
of  the  beam  is  routed  through  the  same  chopper  and  to  the  same  sensor 
without  passing  through  the  cell.  Holes  on  the  chopper  are  so  arranged 
that  the  sensor  is  illuminated  alternately  and  at  equal  intervals  by  the 
two  beams.  In  this  unique  arrangement  a relative  difference  measurement 
of  the  light  extinction  due  to  aerosols  in  the  cell  can  be  made  which  is 
independent  of  sensor-sensitivity  and  laser-intensity  drift. 

Figure  C- 1 1-2  shows  a block  diagram  of  the  phase-locked  demultiplexer 
and  pulse  conditioner  which  provides  analog  outputs  of  the  light  extinc- 
tion Al  due  to  the  aerosols,  and  the  initial  intensity  I of  the  beam. 

The  accuracy  of  this  circuit,  which  resolves  aI  to  0.01%  of  IQ,  is  limited 
by  A3. 

The  other  sensor  attached  to  the  middle  of  the  cell  has  a precise 
cosine  response  and  measures  the  flux  of  light  scattered  by  the  aerosol 
particles  in  the  cell.  Since  the  two  sensors  provide  extinction  and  scat- 
tering measurements  in  the  same  volume  of  aerosol,  an  in-si tu  measurement 
of  the  particles'  absorption  coefficient  is  provided. 

Under  most  atmospheric  conditions  the  extinction  of  light  by  aerosols 
in  the  cell  would  be  insufficient  to  give  an  accurate  measurement.  The 
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measurements  In  the  cell  only  became  possible  by  filling  the  cell  with 
concentrated  ambient  aerosol  particles.  A centrifugal  aerosol  concentra- 
tor  enhances  the  concentration  and  the  aerosol's  extinction  and  scattering 
coefficients  in  the  cell  by  a factor  of  about  20.  It  is,  of  course,  a re- 
quirement of  the  concentrator  to  enhance  the  particle  volume  density  with- 
out changing  the  relative  particle  size  distribution,  since  otherwise  the 
optical  parameters  of  the  aerosol  will  not  be  scaled  up  properly.  If  this 
is  done  accurately,  the  effective  pathlength  in  the  cell  increases  to  400 
meters.  Figure  C-II-3  shows  the  position  of  the  concentrator  with  respect 
to  the  cell. 

Systems  are  shown  in  Figure  C-II-3  which  flow  aerosol,  filtered  air, 
and  Freon  through  the  cell.  During  operation  the  cell  is  cycled  with  fil- 
tered air  and  concentrate  in  order  that  the  influence  of  the  beam- steering 
optics  in  the  system  can  be  subtracted  from  the  measured  aerosol  optical 
properties.  The  Freon  flow  is  used  to  calibrate  the  cosine  sensor  with  a 
known  scattering  coefficient.  The  cell  was  suspended  with  a shock-mounted 
universal  joint  from  a ceiling  beam  in  the  ship  in  order  that  the  influ- 
ence of  the  ship's  motions  on  the  measurements  were  minimized.  The  ver- 
tical positioning  of  the  cell  also  minimized  sedimentation  losses  in  the 
cel  1 . 

3.  DATA 

The  data  collected  during  the  cruise  is  summarized  in  Table  C-II-1. 

Each  of  the  values  of  AI  and  I consist  of  24  point  averages  of  5 minute 
measurement  intervals.  The  extinction  coefficient  calculated  from  ai  and 
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IQ  is  given  by  b,  and  the  coefficient  of  variation  of  b is  given  by  100  x 
s(b)/b  where  s(b)  is  the  sample  standard  deviation  of  b.  The  ratio  of 
F/Fc  is  the  average  of  24  point  values  of  the  flux  measured  by  the  cosine 
sensor  divided  by  the  Freon  calibration  flux,  and  the  aerosol  scattering 
coefficient  measured  in  the  cell  is  given  by  b$.  The  scattering  coeffi- 
cient b s ( MR I ) measured  with  the  MRI  integrating  nephelometer  is  also  given. 
Comparisons  of  b$  and  b s ( MR I ) must  consider  wavelength  compensation,  since 

o 

b^  is  for  x = 6328  A and  bs(MRI)  is  for  a quasi-photopic  source  cen- 

O 

tered  at  5250  A.  The  final  column  gives  values  of  the  aerosol's  absorp- 
tion coefficient  b3  divided  by  b.  The  occasions  during  the  cruise  when 

a 

aerosol  generated  by  the  ship  influenced  the  measurements  are  not  noted. 
Measurements  of  the  aerosol  size  distribution  in  the  ambient  air  and  in 
the  concentrate  were  made,  a particle-size  bandpass  of  somewhat  more  than 
one  order  of  magnitude  and  concentration  enhancement  of  about  a factor  of 
20  were  found. 

4.  RESULTS  AND  CONCLUSIONS 

O 

1)  Transmittance  measurements  at  6328  A were  made  with  a mean 
resolution  of  ± 0.16%  during  the  1977  EOMET  cruise.  This  corre- 
sponds to  a resolution  in  aerosol-optical-depth  measurements  in 
the  cell  of  ± 1.58  x 10"3. 

2)  The  deterioration  of  the  transmittance  measurements  during  the 
cruise  from  the  laboratory  measurements  were  due  to  mechanical 
weaknesses  in  the  cell.  Flexing  of  the  cell  due  to  ambient-cell 
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temperature  differences  and  longitudinal  stresses  from  wave- 
induced  ship's  motion  caused  beam  wander  on  the  transmitted-1 ight 
sensor. 

3)  An  improved  cell  could  better  the  resolution  of  the  transmittance 
measurements  by  a factor  of  10. 

4)  It  may  be  feasible  with  the  present  technique  to  measure  infrared 
extinction  by  aerosols  for  wavelengths  as  long  as  10.6  um  under 
hazy  atmospheric  conditions. 

‘ 6)  The  reciprocal  integrating  nephelometer  (cosine  sensor)  used  on 

the  cell  is  a practical  alternative  to  the  standard  nephelometer. 

6)  The  aerosol  concentrator  produces  a particle-size  bandpass  which 
scales  up  the  values  of  the  scattering  and  extinction  coeffi- 

o 

cients  (at  6328  A)  measured  in  the  cell  from  values  in  the 
atmosphere  with  an  uncertainty  of  about  10%. 

7)  The  combination  of  extinction  and  scattering  measurements  in  the 
same  aerosol  provides  a unique  means  of  measuring  the  particles' 
absorption  coefficient  in  real  time. 
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Table  C-II-1.  Aerosol  Optical  parameters  Measured  During  the 
1977  EOMET  Cruise  of  the  DSNS  HAYES 
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RAMAN  LIDAR  MEASUREMENT  OF  ATMOSPHERIC 


WATER  VAPOR  PROFILES  DURING  THE 
MAY-JUNE  1977  USNS  HAYES 
ATLANTIC/MEDITERRANEAN  CRUISE 
by 

Donald  A.  Leonard  and  Bernard  Caputo 

Computer  Genetics  Corporation,  Wakefield,  Massachusetts  01880 

Stuart  Gathman 

Naval  Research  Laboratory,  Washington,  D.C.  20375 
1.  Introduction 

The  Computer  Genetics  Corporation  (CGC)  participated  in  the  May-June 
1977  EOMET  cruise  of  the  USNS  Hayes  for  the  purpose  of  operating  a Raman 
lidar  system  to  experimentally  measure  vertical  profiles  of  atmospheric 
water  vapor,  temperature  and  visibility  in  the  marine  boundary  layer. 

This  note  reports  some  of  the  water  vapor  results. 

The  basic  Raman  method  for  atmospheric  humidity  profile  measurement 
has  been  well  explored  by  both  laboratory  research  and  field  tests.  Con- 
sequently this  technique  requires  no  further  feasibility  testing  and  has 
also  been  investigated  with  sufficient  thoroughness  so  that  successful 
field  systems  such  as  that  reported  on  here  can  be  designed.  The  princi- 
ples of  atmospheric  humidity  profile  measurements  by  Raman  scatter  have 
been  described  by  Melfi  (1969),  Derr  and  Little  (1970)  and  by  Strauch, 
et  al.  (1972). 

The  primary  result  of  this  experiment  was  to  further  advance  the  de- 
velopment of  remote  water  vapor  instrumentation  and  provide  for  the  first 
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time  a well  documented  field  demonstration  of  a shipboard,  semi-ruggedized, 
single  ended  remote  water  vapor  measurement  tool  under  several  weeks  of 
typical  marine  operational  environments. 

2.  Description  of  Raman  Lidar  Equipment 

The  atmospheric  Raman  spectra  were  obtained  using  a Raman  lidar  sys- 
tem originally  designed  and  used  as  a single  ended  transmi ssometer  (see 
Leonard  and  Caputo,  1974)  and  recently  modified  by  the  authors  to  measure 
atmospheric  profiles  of  water  vapor  and  temperature. 

A block  diagram  of  the  equipment  is  shown  in  Figure  C-1II-1.  The 
laser  source  was  a pulsed  nitrogen  laser  operating  at  337.1  nanometers  and 
producing  pulses  of  10U  kilowatts  peak  power  with  an  effective  pulse  dura- 
tion of  10  nanoseconds  at  a pulse  repetition  rate  of  500  Hertz.  The  ni- 
trogen laser  uses  an  unstable  resonator  cavity.  The  output  from  the  laser 
is  passed  through  an  interference  filter  which  passes  the  337.1  nm  laser 
line  with  high  efficiency  but  blocks  spontaneous  emission  occurring  in  the 
laser  gas  discharge  in  the  Raman  spectral  region. 

The  photons  collected  by  the  receiver  optics  are  passed  through  a 
double  1/4  meter  focal  length  scanning  spectrometer  having  a 0.5  nanometer 
spectral  resolution.  The  wavelength  scanning  of  the  spectrometer  is  con- 
trolled by  the  computer.  The  function  of  the  spectrometer  is  to  produce  a 
spectral  scan  of  the  Raman  air  spectrum  with  sufficient  resolution  so  that 
the  02,  and  H20  Raman  lines  can  be  resolved  and  the  background  measured. 
A 3-channel  filter  wheel  assembly  was  also  employed  so  that  the  peak  sig- 
nals of  the  02,  N2  and  H20  Raman  lines  could  be  detected  with  higher 
efficiency. 
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The  Raman  photon  output  from  the  spectrometer  is  detected  by  an  RCA 
8850  photomultiplier.  The  Raman  photoelectron  signal  and  a gate  pulse 
suitably  synchronized  and  delayed  with  respect  to  the  laser  firing  are 
combined  in  ten  (10)  separate  coincidence  gates,  the  output  from  which  are 
recorded  using  direct  memory  access  (DMA)  techniques  with  a Data  General 
NOVA  1220  computer.  The  number  of  pulses  occurring  during  the  same  data 
taking  interval  is  also  recorded.  The  ratio  of  the  Raman  photons  to  laser 
pulses  recorded  over  a given  time  interval  is  the  basic  data  obtained. 

The  ten  coincidence  gates  provide  Raman  data  at  ten  separate  locations. 

The  data  taking  is  fully  automated.  The  computer  can  be  directed 
through  the  teletype  (TTY)  to  specify  the  number  of  laser  pulses  over 
which  to  record  data  at  each  wavelength  and  the  maximum  and  minimum  wave- 
lengths and  the  wavelength  increments  of  each  spectral  scan.  The  data 
obtained  are  printed  on  the  TTY  in  tabular  form  and  also  recorded  on 
cassettes. 

3.  Data  Obtained 

Experiments  were  performed  with  both  the  1/4-m  scanning  spectrometer 
and  the  3-channel  filter  wheel.  Data  were  obtained  as  follows: 

3 . 1 Spectral  Scan  of  Op,  No  and  HpO  Raman  Lines. 

The  purpose  of  the  spectral  scans  with  the  1/4  meter  SPEX  double  spec- 
trometer of  the  02,  N2  and  H20  Raman  lines  is  to  determine  the  signal- 
to-noise  ratio  of  the  atmospheric  Raman  lines  with  respect  to  the  back- 
ground spectrum  both  ambient  and  laser  induced.  The  spectrum  plotted  on  a 
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logarithmic  scale  in  Figure  C-III-2  shows  typical  raw  laser  induced  data 
obtained  from  the  atmosphere  at  wavelengths  between  345  and  390  nanometers 
with  1.0  nanometer  resolution.  The  wavelength  scale  is  uncalibrated  and 
exhibits  a bias  of  about  3 nm,  i.e.  the  exact  wavelength  of  the  line  is 
365.8  nm.  For  this  data  set  the  range  resolution  was  10  meters  and  the 
altitude  was  100  meters.  A measurement  time  of  15  seconds  per  data  point 
was  used.  As  the  spectrum  of  Figure  C-III-2  shows  the  0^,  N2  and  H^O 
Raman  lines  appear  typically  as  well  resolved  spectral  lines  with  the 
basic  signal-to-noise  ratios  of  100,  250  and  20  respectively.  The  noise 
level  is  essentially  the  dark  current  of  the  photomultiplier. 

The  conclusion  to  be  drawn  is  that  the  Raman  spectrum  of  the  atmo- 
sphere in  the  cases  investigated  is  free  of  laser  induced  interference  and 
that  the  relative  concentrations  of  02,  K2  and  H20  may  be  determined  by 
simply  measuring  the  peak  intensity  at  each  spectral  line. 

3.2  Profiles  of  HqO/Nq  in  the  Marine  Boundary  Layer. 

Atmospheric  profiles  of  the  H20/N2  Raman  line  ratio  were  obtained  by 
using  interference  filters  in  place  of  the  1/4-meter  double  spectrometer. 
Due  to  the  larger  optical  efficiency  of  the  filters  relative  to  the  spec- 
trometer higher  signal  levels  were  obtained  and  thus  measurements  at 
higher  altitudes  were  possible.  The  results  of  the  spectral  scans,  dis- 
cussed in  Section  3.1  above,  justify  the  use  of  filters. 

Data  was  normally  obtained  simultaneously  at  10  distinct  and  separate 
altitudes.  The  Raman  return  was  sequentially  sampled  at  the  02,  N2  and  H20 
Raman  line  positions  using  the  filters.  The  ratio  of  the  H20  to  the  N2 
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Raman  returns  was  used  to  determine  the  H20  partial  pressure  assuming  that 
the  N2  partial  pressure  corresponded  to  standard  atmospheric  conditions. 
This  normalization  of  the  H20  signal  to  the  N2  signal  effectively  cancelled 
out  effects  due  to  variable  atmospheric  transmission,  laser  power  fluctu- 
ations and  changes  in  the  overall  photoelectric  efficiency  of  the  trans- 
mitter/receiver system  and  gave  the  H20  to  N2  mixing  ratio  directly.  The 
02  Raman  line  was  also  monitored  and  the  ratio  of  the  02  to  N2  Raman  sig- 
nals was  checked  as  an  internal  system  redundancy  procedure. 

As  typical  of  the  data  obtained.  Figures  C-III-3,  -4,  -5,  and  -6  show 
as  a function  of  altitude  the  partial  pressure  of  water  vapor  in  the  atmo- 
sphere in  millibars  (mb)  as  measured  by  the  CGC  Raman  lidar  on  the  Hayes. 
These  marine  boundary  layer  water  vapor  profiles  are  representative  of 
both  the  cold  Atlantic  waters  off  the  northeastern  coast  of  North  America 
as  well  as  the  warmer  eastern  Atlantic  and  Mediterranean  type  waters. 

In  situ  measurements  of  water  vapor  profiles  in  the  marine  boundary 
layer  throughout  the  entire  cruise  were  obtained  by  means  of  a series  of 
vertical  soundings  of  atmospheric  pressure  altitude,  dry  air  temperature 
and  wet  bulb  temperature  using  a kite-balloon  (see  Gathman,  1978). 

Error  bars  are  shown  for  the  Raman  data  in  Figures  C-III-3  through 
C-III-6  which  correspond  to  a + 10%  accuracy  level  based  on  the  photon 
statistics  obtained  in  a typical  measurement  interval  usually  the  order 
of  one  minute.  However,  each  typical  profile  plotted  represents  the 
average  Raman  data  over  a time  period  the  order  of  one  half  to  one  hour 
so  that  comparisons  could  be  made  with  the  balloon-borne  ground  truth 
instrumentation. 
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Very  favorable  ground  truth  comparisons  are  indicated  for  the  pro- 
files shown  in  Figures  C-III-3  through  C-Ill-6  where  simultaneous  balloon- 
borne  wet  and  dry  bulb  temperature  measurements  agree  well  with  the  remote 
sensing  technique.  Ground  truth  was  not  available  for  all  the  Raman  pro- 
files because  under  certain  conditions,  such  as  high  winds,  the  balloon- 
borne  instrumentation  was  not  operable. 

4.  Conclusions 

The  Raman  lidar  was  operated  aboard  the  Hayes  for  a period  of  three 
weeks  and  produced  range  resolved  02,  N2  and  H20  atmospheric  profile  data 
in  accordance  with  its  design  expectations  to  an  altitude  of  0.5  kilometer. 
Ground  truth  obtainable  for  the  H20  profiles  showed  excellent  agreement 
in  every  case  in  which  comparisons  were  available.  This  work  thus  demon- 
strated that  a Raman  lidar  can  successfully  measure  atmospheric  water 
vapor  profiles  for  an  extended  period  of  time  in  typical  operational 
marine  environments. 
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Fig.  C-III-6  Partial  Pressure  of  Water  Vapor  in 
Atmosphere  in  Millibars  (mb)  as  a 
Function  of  Altitude  in  Meters  as 
Measured  by  the  Raman  Lidar  on  the 
HAYES,  4 June  1977  (Mediterranean) 
Calibrated  by  Web  Bulb/Dry  Bulb 
Measurements. 
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MEASUREMENT  OF  CN  ON  BOARD  THE  USNS  HAYES  CRUISE 
OF  MAY-JUNE  1977 


John  N.  Hayes 
Atmospheric  Physics  Branch 
NRL 


The  purpose  of  this  investigation  was  to  determine  the  approximate 
distribution  of  ship  induced  condensation  nuclei,  CN,  about  the  ship  and 
their  concentration  as  a function  of  external  weather  conditions. 


The  CN  measurements  were  made  with  a commercially  available  G.E.  CN 
counter.  About  40  locations  on  the  ship  were  "staked"  as  standard  where 
the  measurements  were  repeated  as  often  as  time  and  the  endurance  of  the 
experimenter  would  allow.  Figure  C-IV-1  is  a schematic  layout  of  the  ship 
as  seen  from  the  top.  (More  detailed  diagrams,  including  good  side-views, 
can  be  obtained  from  NRL  Report  7370.*)  The  letter  L stands  for  level. 

Thus  POLW  stands  for  stakeouts  POUW,  PQ1W,  P03W,  P03W.  W stands  for  the 
well  of  the  ship.  Not  all  levels  are  equally  accessible  at  each  measure- 
ment location.  Forward  of  the  area  bounded  by  the  rectangle  P0LW--S0LW-- 
POL — SOL,  only  the  01  level  exists.  The  following  key  will  help  the  reader 
decipher  our  notation. 


PB  — Port  Bow 

SB  — Starboard  Bow 

POL  — Port  01,  02,  03,  04  levels 

SOL  --  Starboard  01,  02,  03,  04  levels 

POLW  — Port  Well  01,  02,  03,  04  levels 

SOLW  — Starboard  Well  00,  01,  02,  03,  04  levels 

FOLMW  — Middle  of  forward  section.  00,  02,  04  levels 


! 
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All  positions  with  "A"  in  them  imply  the  aft  section  of  the  ship. 

Not  all  levels  were  existent  there  either.  The  04A  did  not  exist  anywhere. 
S03AW  and  S02AW,  S03A,  and  S02  were  non  existent  as  well.  X00  and  YOU 
were  at  deck  level  with  nothing  above  them. 

Although  data  were  collected  every  day,  only  certain  representative 
data  sets  are  shown  here.  Not  all  were  usable  because  of  drastic  shifts 
in  relative  wind  direction  during  the  course  of  measurement,  instrument 
misbehavior,  or  rough  seas.  The  plotting  of  the  data  is  rather  difficult 
because  of  the  limited  number  of  data  points  and  because  the  plots  should 
be  represented  in  three  dimensions.  The  bow  data  are  from  the  01  level; 
the  forward  contours  are  from  the  04  level;  while  the  aft  data  are  from 
the  02  and  03  levels. 

The  contours  shown  are  CN  (On  Board)/CN  iSea).  A reasonable  amount 
of  interpolation  and  extrapolation  was  necessary  because  of  the  paucity 
of  data  stations,  but  will  not  significantly  affect  the  conclusions  to  be 
drawn  from  the  data. 

The  data  taken  in  the  Mediterranean  are  fewer  in  number  than  that 
taken  in  the  Atlantic.  Several  reasons  account  for  this:  less  time  was 
spent  in  the  Mediterranean;  there  was  a several  day  instrument  failure; 
and  the  operational  procedures  were  far  more  lengthy  there  due  to  the  loss 
of  help. 

Figure  C-lV-2a  shows  04  level  data  only  and  is  therefore  very  simple. 
The  day  (135)  and  time  of  day  are  indicated.  The  exact  location  can  be 
obtained  from  the  cruise  track  appearing  earlier  in  this  report.  The 
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contour  numbers  are  the  natural  airborne  condensation  nuclei  number  den- 


sity, divided  by  CN  density  at  the  measuring  point.  The  relative  wind  is 
coming  roughly  from  the  bow  with  a 20°  spread  over  the  time  of  measurement, 
accompanied  by  a 3 knot  spread  in  relative  wind  velocities.  The  particle 
densities  used  to  generate  the  contours  were  the  geometric  mean  of  four 
measurements  at  each  station.  We  see  from  Figure  C-IV-2a  that  in  the  vi- 
cinity of  the  smokestacks  the  CN  count  is  roughly  400  times  that  of  the 
background  level.  As  large  as  this  number  sounds,  its  associated  optical 
effect  was  not  discernible. 

Figure  C-IV-2b  shows  only  the  03  and  02  levels  in  the  aft  of  the  ship 
for  day  140.  In  spite  of  the  extremely  high  reading,  visibility,  while 
noticeably  reduced,  was  still  good;  however,  the  stack  gas  aerosols  were 
definitely  detectable  by  the  olfactory  senses. 

Figure  C-IV-2c  shows  the  contours  for  04  level  (forward)  and  03  and 
02  levels  a ft  combined  for  the  wind  conditions  indicated.  Because  of  the 
prevailing  winds  on  the  cruise  track  and  the  speed  of  the  ship,  these  were 
not  uncommon  conditions  on  the  trip.  In  order  to  get  the  relative  wind  to 
be  from  the  bow  to  obtain  natural  aerosol  samples,  it  was  frequently  neces- 
sary to  reverse  the  course  of  the  ship.  We  note  the  low  count  of  CN  (Sea) 
and  the  high  ratio  between  the  stacks.  Such  ratios  often  penetrated  deep 
into  the  well  in  which  the  S.R.I.  LIDAR  was  operating.  Even  though  the 
aerosol  sizes  were  of  the  same  order  as  the  wavelength  of  lasar  light 
used  in  the  LIDAR  system,  they  did  not  appear  to  affect  its  operation. 
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The  preceding  figures  as  well  as  C-IV-2d,  -2e,  and  -2f  show  what  was 
fully  expected  in  almost  all  respects.  The  stack  exhausts  fully  dominated 
the  data  collected  downwind,  as  did  the  galley  exhausts  which  on  occasion 
could  be  separated  from  the  stack  exhaust  in  our  data.  The  data  showed 
that: 

The  ship  should  sail  into  the  prevailing  winds  to  obtain  uncon- 
taminated aerosol  data  from  centers  located  at  the  bow  of  the  ship. 

The  port  and  starboard  bow  were  by  far  the  least  contaminated 
locations  on  the  ship. 

Stack  and  galley  exhaust  aerosols  dominated  the  submicron  aerosol 
region  so  totally  in  the  wind  shadow  that  the  natural  airborne  sea  salt 
particles  could  not  be  discerned. 

Turbulence  behind  the  wind  screen  on  the  flying  bridge  allowed 
for  contamination  of  the  natural  aerosols.  For  reliable  data,  all  gather- 
ing probes  should  be  stationed  forward  of  the  flying  bridge  wind  screen. 

In  addition,  we  sought  vertical  profile  data  which  are  shown  for  sev- 
eral levels  at  given  stations  in  the  remaining  figures  C-IV-3a,  -3b,  -3c, 
and  -3d  where  the  x axis  is  the  CN  measurement  in  #/cc  and  the  y axis  is 
the  deck  level.  They  show  clearly  that  uncontaminated  results  can  be  ob- 
tained, especially  in  high  winds,  if  the  collector  is  situated  on  the 
windward  side  of  the  ship;  but  if  data  are  collected  from  the  leeward  side, 
they  are  untrustworthy.  Figure  C-IV-4  tends  to  summarize  all  the  results 
together  with  deviations  from  the  average  effects.  The  open  ocean  gives 
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uncontaminated  aerosol  counts  of  the  order  of  100-1000  nuclei  per  cc.  In 
the  wind  shadow,  we  get  (depending,  in  part,  upon  how  the  engines  are  op- 
erating) 105  counts/cc  and  higher.  Intermediate  results  generally  involve 
turbulent  mixing  or  fresh  air  entrainment. 
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DAY  140  WIND  : 0 to  345 


fig.  (C-IV-3)  Vertical  profile,  CN  on  ship. 


DAY  142  12  00  2 - 1 40  02;  WIND:  10  to  27  KT,  0-2  20 


Fig.  (C-IV-4)  Summary  of  ship  effects  on  CN,  with  deviation  from 
average  effects. 


SALT  LOAD  DATA 


David  J.  Bressan 
Code  833D 


Salt  load  data  was  collected  at  different  points  on  the  HAYES,  desig- 


nated by  the  code  shown  below. 


90°  Rel 


0°  Rel- 


180*  Rel 


FIG.  C-V-l  SALT  SAMPLING  STATIONS 

The  locations  and  directions  were  chosen  to  provide  for  isokinetic  en- 
trance into  the  sampler. 

An  impactor  followed  by  a fine  screen  was  used  to  collect  the  salt 
samples.  It  had  previously  been  found  that  the  impactor  collects  parti- 
cles above  1-2  u diameter,  and  the  screen  collects  particles  greater  than 
0.1-0. 2 u.  Some  salt,  especially  the  larger  diameter  particles,  is  depos- 
ited in  the  8"  long  pipe  between  the  isokinetic  opening  and  the  impactor 
surface.  The  amount  deposited  there  tends  to  be  about  1/2  of  the  total 
collected  mass.  The  amount  of  salt  collected  was  determined  by  washing 
the  collectors  in  a known  amount  of  distilled  water  and  then  measuring  the 


conductivity  of  the  water  and  calculating  ug  salt/nr  air. 
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A comparison  between  the  impactor  system  and  three  of  the  fine  mesh 
screens  mounted  in  tandem  showed  that  the  total  salt  load  computed  from 


the  sum  of  that  on  the  three  screens  agreed  to  within  10-15%  with  the  total 
salt  retrieved  from  the  impactor  system. 

The  open  areas  of  the  screens  are  about  45%  of  the  total , so  that  the 
screen  retains  about  55%  of  the  material  0.1  j i diameter  which  reaches 
it.  It  should  be  then  noted  that  the  relative  humidity  of  the  air  will 
affect  the  size  of  the  salt  particles  collected,  as  will  passage  through 
an  impactor  nozzle,  since  Bernoulli  pressure  changes  will  also  affect  the 
relative  humidity. 

For  these  reasons  only  the  raw  data  is  presented.  It  is  shown  simply 

3 

as  that  collected  by  the  screen,  slide  and  inlet  system  in  pg  salt/M  of 
air.  It  is  expected  that  the  inlet  system  collects  mostly  particles  larger 
than  30-50  v diameter,  but  without  a sharp  cut  off  point.  The  remaining 
particles  above  ^ 1.0  u are  collected  on  the  slide,  and  the  screen  will 
collect  about  1/2  of  the  total  mass  reaching  it,  whose  range  should  be  0.1- 
1.0  y diameter. 

Interpretation  of  this  data  can  best  be  done  with  respect  to  changes 
in  wind  speed  and  white  cap  coverage,  and  changes  in  air  masses  with  ma- 
rine trajectories  and  turbulence  being  reflected  in  the  salt  load.  These 
observations  should  be  available  in  the  data  file. 
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Table  C-V- 

-1  Salt 

Load  Data 

(page  1 of 

2) 

pq  sa.It/M3 

Date 

Time 

Location 

Inlet 

Slide 

Screen 

Remarks 

05/16/77 

10-.30L 

l-C 

— 

43.5 

FW  170° 

05/16/77 

U:30L 

03-E 

— 

11.8 

05/17/77 

12:00L 

04-A 

— 

17.5 

7.9 

05/17/77 

13:00L 

1— D 

— 

30.6 

5.9 

05/17/77 

14:30L 

03-F 

— 

16.5 

9.6 

BW  240° 

05/17/77 

15-.45L 

04-A 

— 

16.3 

5.6 

FW  00 0° 

05/18/77 

11:50L 

04-A 

— 

1.8 

2.4 

No  white  caps 

05/18/77 

12:15L 

01-M 

— 

1.8 

4.3 

(Air  mass  change?) 

05/18/77 

16.-45L 

04-A 

— 

5.9 

3.0 

Many  white  caps 

5/18/77 

17-.10L 

01-M 

— 

9.0 

- 

05/19/77 

11:20L 

01-M 

— 

1.0 

2.4 

05/19/77 

13: SOL 

l-C 

— 

1.9 

2.3 

05/19/77 

16:10L 

04-A 

— 

2.5 

1.8 

05/19/77 

16:30L 

02-C 

— 

1.9 

3.5 

05/21/77 

15 :25L 

03-H 

— 

14.8 

4.1 

05/21/77 

23:152 

04-A 

— 

12.2 

1.1 

05/21/77 

23:30Z 

04-A 

— 

9.4 

1.2 

05/21/77 

23 :55Z 

04-A 

— 

11.2 

- 

05/22/77 

12:05Z 

04-A 

— 

15.1 

1.4 

Note  when  wind  speed  changes 

05/22/77 

12:20Z 

04-A 

— 

18.3 

2.1 

on  22nd 

05/22/77 

12:40Z 

04-A 

— 

15.9 

2.4 

05/22/77 

17 :15Z 

04-A 

— 

10.1 

1.3 

05/22/77 

17-.30Z 

04-A 

— 

10.3 

1.6 

05/22/77 

17:452 

04-A 

— 

10.3 

0.9 

05/24/77 

16:45Z 

04-A 

— 

3.2 

1.2 

05/25/77 

10:55Z 

04-A 

— 

3.3 

2.6 

05/25/77 

11:15Z 

04-A 

— 

6-8 

2.1 

Quite  a few  snail  white  caps 

5/25/77 

15:20Z 

04-A 

— 

3.3 

1.3 

5/25/77 

15:50Z 

04-A 

— 

3.8 

1.3 

5/26/77 

09:00Z 

03-G 

— 

2.5 

1.5 

5/26/77 

09:2CZ 

03-G 

— 

1.9 

1.3 

5/26/77 

14:40Z 

03-G 

— 

2.7 

1.5 

5/26/77 

14:50Z 

03-G 

— 

2.1 

1.8 

5/26/77 

15:15Z 

1-E 

— 

1.6 

2.2 
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Table 

C-V-l  Salt 

Load  Data 

(page  2 of 

ug  salt/M3 

2) 

Date 

Time 

Iecaticn 

Inlet 

Slide 

Screen 

05/27/77 

09:25Z 

03-H 

— 

10.9 

2.6 

05/27/77 

09:45Z 

03-H 

— 

11.2 

3.7 

05/31/77 

08:35Z 

03-G 

10.8 

8.0 

2.3 

5/31/77 

08:55Z 

03-G 

— 

7.9 

2.2 

05/31/77 

08:15Z 

03-G 

10.8 

8.2 

2.1 

06/01/77 

08:50Z 

03-G 

10.9 

11.7 

1.5 

06/01/77 

09:15Z 

03-G 

10.9 

11.9 

1.4 

06/01/77 

09:30Z 

03-G 

10.9 

13.1 

1.4 

06/02/77 

08:30Z 

03-G 

4.1 

1.5 

1.9 

06/02/77 

08 :35Z 

03-G 

4.1 

1.3 

2.5 

06/03/77 

U:35Z 

03-G 

3.6 

4.3 

2.0 

06/03/77 

11:55Z 

03-G 

3.6 

4.7 

2.0 

06/03/77 

12:10Z 

03-G 

3.6 

4.3 

1.8 

06/04/77 

06 :30Z 

03-H 

^lO.S 

M5.6 

3rd5.3 

06/04/77 

06:40Z 

03-H 

8.0 

8.7 

2.5 

06/04/77 

06:55Z 

03-H 

9.0 

10.2 

4.8 

06/04/77 

06:55Z 

03-H 

9.0 

10.2 

4.8 

06/04/77 

14:10Z 

03-G 

12.8 

11.1 

2.9 

06/04/77 

14:25Z 

03-G 

12.8 

9.9 

2.9 

06/05/77 

08.-10Z 

03-G 

7.8 

9.4 

1.5 

06/05/77 

08:30Z 

03-G 

7.8 

8.6 

2.3 

06/05/77 

20:30Z 

03-G 

5.2 

7.2 

1.8 

06/05/77 

20-.45Z 

03-G 

5.2 

6.3 

2.5 

06/16/77 

11:30Z 
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INTRODUCTION 

During  the  period  15  May  - 7 June  1977  the  SRI  Mark  IX  mobile  lidar 
system  was  operated  on  board  the  USNS  Hayes  on  a scientific  cruise  across 
the  Atlantic  Ocean  and  the  Mediterranean  Sea.  The  purpose  was  to  obtain  a 
continuous  record  of  the  vertical  aerosol  and  cloud  structure  up  to  an  al- 
titude of  3.5  km. 

MARK  IX  MOBILE  LIDAR  SYSTEM 

The  SRI  Mark  IX  ruby  lidar  system  is  the  product  of  14  years  of  ex- 
perience in  developing  and  applying  laser  radar  systems  to  the  observation 
of  atmospheric  aerosol,  cloud,  and  precipitation  distributions  over  remote 
areas  with  high  spatial  and  temporal  resolution.  The  lidar  system  is  op- 
erated from  a 6 m van  provided  with  power-generating  facilities  to  allow 
completely  mobile  operation  or  immediate  operation  at  any  stationary  loca- 
tion. One  of  the  Mark  IX' s principal  features  is  a digital  data  recording, 
processing,  and  display  system  for  real-time  TV-viewing  of  controlled  gray- 
scale presentations  of  atmospheric  structure.  These  gray-scale  displays 
are  generated  from  backscatter  data  obtained  from  multiple  laser  firings 
while  the  lidar  is  being  scanned  in  elevation  or  azimuth,  or  repetitively 
pulsed  with  time  or  distance.  In  addition,  the  minicomputer-based  system 
permits  the  determination  of  quantitative  estimates  of  aerosol  densities. 

Views  of  the  Mark  IX  system,  a block  diagram,  and  lidar  specifications  are 
presented  in  Figure  C-VI-1. 
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DATA  ACQUISITION  PROGRAM 


During  the  1977  USNS  Hayes  cruise  the  vertical  aerosol  structure  was 
sampled  at  a rate  of  approximately  three  profiles  per  minute  for  about 
twenty  hours  each  day.  Generally,  each  backscatter  signature  (profile) 

O 

was  digitized  in  real  time  with  a sampling  rate  of  5 x 10  sec  (range 
resolution  of  7.5  m)  from  the  lidar  height  aboard  ship  to  an  altitude  of 
3.6  km.  Periodically,  however,  profiles  were  digitized  with  a sampling 
rate  of  20  x 10  sec  (range  resolution  of  30  m)  up  to  an  altitude  of 
14.4  km  in  order  to  monitor  upper  level  clouds,  primarily  cirrus.  Concur- 
rent with  each  lidar  firing,  the  output  of  an  MRI  (Meteorology  Research, 
Inc.)  integrating  nephelometer  was  sampled  and  digitized  with  a low-speed 
analog-to-digital  converter.  The  lidar  and  nephelometer  data  along  with 
date  and  time  information  were  recorded  on  nine-track  magnetic  tape  for 
further  analysis. 

LIDAR  OBSERVED  ATMOSPHERIC  STRUCTURE 

A complete  set  of  Polaroid  photographs  depicting  intensity-modulated 
height/time  TV  displays  of  vertical  atmospheric  structure  along  the  cruise 
route  is  presented  in  the  S.R.I.  appendix.  Two  of  these  photographs  are 
shown  in  Figure  C-VI-2  in  order  to  explain  the  scale  labels  and  types  of 
backscatter  returns  observed.  Tick  marks  equally  spaced  at  ten  minute  in- 
tervals at  the  top  of  each  Polaroid  frame  delineate  time.  The  photographs 
are  arranged  in  rows  and  numbered  sequentially  from  the  beginning  of  the 
cruise.  The  date,  frame  number,  and  times  'GMT)  corresponding  to  the  first 
and  last  ticks,  respectively,  of  each  frame  in  a particular  row  are  tabu- 
lated to  the  left  of  each  row.  Periodic  time  gaps  due  to  missing  data  and 
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resulting  in  missing,  extra,  or  unequally-spaced  tick  marks  are  also  docu- 
mented. Height  above  the  lidar  is  indicated  in  each  frame  by  bright  hori- 
zontal lines  spaced  at  one  kilometer  intervals.  Gross  discontinuities  in 
the  vertical  aerosol  structure  appearing  within  a single  frame  are  artifi- 
cial: they  are  due  to  a reduction  in  the  height  scale  resulting  from  an 
increase  in  the  range  (altitude)  of  the  lidar  beam  in  order  to  detect  high 
altitude  clouds,  as  noted  previously.  In  these  cases,  the  vertical  scale 
has  been  reduced  by  a factor  of  4,  and  the  bright  horizontal  lines  appear 
at  4 km  intervals  (see  Figure  C-VI-2b). 

Picture  brightness  is  proportional  to  the  logarithm  of  observed  back- 
scatter.  The  TV  display  provides  16  gray  steps,  as  shown  by  the  gray  scale 
presented  beside  each  data  frame.  During  data  processing,  the  computer 
operator  selected  a signal  level  for  each  gray  step  to  emphasize  structure 
of  the  near-surface  aerosol  layer.  This  resulted  in  saturation  of  the  dis- 
play for  very  large  backscatter  values,  principally  related  to  clouds  and 
precipitation.  When  the  maximum  brightness  was  exceeded  for  these  data, 
the  computer  recirculated  the  data  to  the  lowest  brightness  gray  step. 

This  technique  provides  digital  isoecho  contouring  evident  in  the  cloud 
and  precipitation  returns  (see  Figure  C-VI-2b). 

The  0 to  5 volt  output  of  the  MR I integrating  nephelometer,  as  sampled 
at  the  time  of  each  lidar  firing,  is  plotted  at  the  top  of  the  height/time 
cross  section  (see  Figure  C-VI-2b).  The  0-V  level  corresponds  to  the  3-km 
range  mark  and  a 3.5-V  calibration  signal  is  normally  inserted  once  for 
each  experimental  run  (for  example,  at  1616  GMT  in  Frame  38).  Using  infor- 
mation provided  by  the  manufacturer,  the  3.5-V  level  corresponds  to  an 


3 

aerosol  concentration  of  265  pg/m  and  provides  a concentration  scale  for 
the  data  plot. 

RESEARCH  OBJECTIVES 


. 


Detailed  analyses  of  the  Rictorial  displays  and  digital  data  which  are 
presently  under  way  will  enable  us  to  characterize  the  marine  atmospheric 
aerosol  structure  in  terms  of  aerosol  mixing  layer  heights  and  relative 
columnar  aerosol  backscatter  values.  It  is  expected  that  planned  correl- 
ative analyses  between  these  Mark  IX  lidar-derived  aerosol  layer  heights 
and  backscatter  values  and  the  shipboard  observations  of  various  meteoro- 
logical parameters  will  also  contribute  to  our  understanding  of  the  atmo- 
spheric marine  boundary  layer. 
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Fig.  (C- VI-2) 


1910  GMT  Cirrus  cloud 


EXAMPLES  OF  LIDAR  OBSERVED  HEIGHT/TIME  SECTIONS 
OF  VERTICAL  AEROSOL,  CLOUD,  AND  PRECIPITATION  STRUCTURE. 

Frame  82  exemplifies  a height/time  display  containing  a time  gap.  In  addition 
to  detailed  scale  labeling,  the  nephelometer  data  trace,  a low  level  stratus  cloud 
deck,  and  the  gray  scale  are  all  denoted. 

Frame  123  contains  a height  scale  change  at  2000  GMT.  Note  the  cirrus  cloud 
return  between  8-10  km.  Both  the  bright  echoes  and  embedded  dark  echoes 
below  1 km  are  caused  by  rain.  These  embedded  dark  echoes,  which  also 
appear  in  the  cloud  returns,  are  examples  of  the  isoecho  contouring  effect 
explained  in  the  text.  Two  distinct  aerosol  layers  are  also  evident  below  1 km. 
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Appendix 

LIDAR-OBSERVED  AEROSOL  STRUCTURE 
ABOVE  EOMET  77  CRUISE  ROUTE 
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SATELLITE  OBSERVATIONS  OF  ATMOSPHERIC  AEROSOLS 


DURING  THE  EOMET  CRUISE 

M.  Griggs 

SCIENCE  APPLICATIONS,  INC. 

La  Jolla,  California  92038 

28  August  1978 

ABSTRACT 

Measurements  of  the  atmospheric  aerosol  optical  thickness  were  made 
during  the  1977  EOMET  cruise  across  the  Atlantic  Ocean  and  the  Mediter- 
ranean Sea.  These  data  were  obtained  at  the  same  time  as  NOSS-5  and  GOES 
visible  radiance  measurements  in  the  vicinity  of  the  ship.  Linear  rela- 
tionships between  the  upwelling  radiance  and  the  aerosol  optical  thickness 
were  found  for  each  satellite,  confirming  earlier  Landsat  results.  Dif- 
ferences in  the  relationships  for  each  satellite  are  attributed  to  differ- 
ences in  the  radiometric  calibrations  of  the  satellite  sensors. 


-291- 


r 

1.  Introduction 

The  atmospheric  marine  boundary  layer  is  of  considerable  importance 
to  Navy  operations.  Of  particular  interest,  with  the  increasing  use  of 
electro-optical  systems,  is  the  nature  and  distribution  of  aerosols  over 
the  oceans.  In  order  to  provide  further  information  about  the  aerosols, 
and  other  meteorological  parameters  in  the  marine  boundary  layer,  the 
Naval  Research  Laboratory  conducted  an  Electro-Optical  Meteorology  (EOMET) 
cruise.  May  15  to  June  6,  1977,  across  the  North  Atlantic  Ocean  and  Medi- 
terranean Sea.  This  cruise  offered  a good  opportunity  to  further  investi- 
gate the  satellite  technique,  developed  by  Griggs  (1975,  1977)  to  measure 
the  atmospheric  aerosol  content  over  oceans.  This  technique  relates  the 
r upwelling  visible  radiance  measured  by  the  satellite  to  the  atmospheric 

aerosol  optical  thickness.  Since  60%  of  the  aerosols  are  typically  in  the 
lowest  1 km,  and  90%  in  the  lowest  3 km,  it  is  clear  that  the  satellite 
measurements  can  provide  information  of  considerable  importance  to  Navy 
operations.  The  previous  studies  were  based  on  Landsat  data,  but  for  this 
investigation  the  NOAA  and  GOES  satellites  were  used. 

2.  Previous  Landsat  Results 

Previous  results  using  Landsat  1 data  (Griggs,  1975)  and  Landsat  2 
data  (Griggs,  1977)  have  shown  that  a linear  relationship  exists  between 
the  upwelling  radiance  in  the  visible  region  and  the  aerosol  content.  Data 
have  been  obtained  at  several  sites,  the  largest  data  set  being  for  the 
Pacific  Ocean  at  San  Diego  for  Landsat  2 overpasses.  These  results  are 
shown  in  Figure  C-VII-1.  The  radiances  are  determined  from  the  Landsat 
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digital  data  (densitometry  of  the  black  and  white  imagery  is  not  accurate 
enough  for  intercomparison  of  different  images),  and  the  aerosol  content 
values  are  determined  with  ground-based  Volz  sun  photometer  measurements 
at  the  time  of  the  Landsat  overpass.  The  aerosol  content  is  defined  in 
terms  of  the  Elterman  (1964)  model  vertical  aerosol  optical  thickness; 
i.e.,  the  aerosol  content  is  given  by  the  ratio  (measured  aerosol  optical 
thickness  at  wavelength  x to  the  model  aerosol  optical  thickness  at  wave- 
length x)  x N;  i.e.,  a value  of  2N  for  the  aerosol  content  indicates  that 
the  optical  thickness  is  twice  that  of  the  Elterman  model. 

The  relationships  appear  best  for  MSS  5 and  MSS  6;  this  is  due  to  the 
fact  that  the  radiance  in  MSS  4 is  affected  more  by  suspended  matter  in  the 
water.  Figure  C-VII-1  does  not  show  MSS  7 data,  since  the  digital  data  for 
this  channel  are  uncertain  owing  to  NASA  procedures  for  producing  Landsat 
2 computer  compatible  types.  The  measured  relationships  show  excellent 
agreement  with  theoretical  calculations  (Griggs,  1978). 

3.  Approach 

Arrangements  were  made  for  sun  photometer  measurements  of  the  atmo- 
spheric aerosol  optical  thickness  to  be  taken  daily,  weather  permitting, 
on  board  the  U.S.N.S.  Hayes  (EOMET  cruise  vessel)  at  times  as  nearly  coin- 
cident as  possible  with  the  overpasses  of  NOAA-5  (0800-1000  local  standard 
time)  and  at  1600  GMT,  when  GOES-1  digital  data  are  routinely  recorded  and 
stored.  The  radiances  measured  by  the  satellites  must  be  modified  in  or- 
der to  compare  them  with  the  Landsat  values  shown  in  Figure  C-VII-1.  The 
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Landsat  radiances  are  for  the  multispectral  scanner  (MSS)  spectral  band- 
passes  and  for  nadir  viewing.  The  NOAA-5  and  GOES-1  have  rlightly  differ- 
ent bandpasses,  and  in  general  the  radiances  of  interest  are  not  obtained 
in  the  nadir  direction.  These  radiances  are  normalized  to  the  Landsat 
viewing  and  sun  angle  conditions  by  means  of  theoretical  calculations  with 
the  Dave  (1972)  atmospheric  scattering  code,  which  was  previously  used  in 
support  of  the  Landsat  study. 

4.  Data  Analysis 

The  sun  photometer  data  had  to  be  carefully  reviewed  since  measure- 
ments are  difficult  to  make  on  board  ship  due  to  the  motion  of  the  ship. 
This  is  especially  true  for  the  measurement  of  the  airmass;  fortunately 
this  can  be  precisely  calculated  from  knowledge  of  the  location  of  the 
ship  and  the  time. 

The  location  of  the  ship  as  a function  of  time,  and  the  orbital  param- 
eters of  the  NOAA-5,  which  is  in  a polar  orbit,  are  used  to  determine  the 
respective  locations  of  the  ship  and  the  satellite  at  the  time  the  satel- 
lite radiance  measurement  is  acquired.  Then,  using  spherical  trigonometry, 
the  sun  zenith  (0Q)  and  azimuth  (<t>)  angles  and  the  viewing  angle  (0)  from 
the  satellite  to  the  ship  are  calculated.  For  the  geosynchronous  GOES-1, 
which  is  at  a fixed  location,  the  required  angles  are  similarly  calculated. 
These  angles  are  used  as  input  to  the  Dave  code  which  is  used  to  normalize 
the  measured  radiance  to  the  standard  conditions  used  for  the  Landsat  data, 
viz.  0Q  = 63.3°,  <t>  = 0,  0 = 0.  The  calculations  are  made  for  each  pair  of 
satellite  radiance  and  aerosol  content  observations,  using  the  measured 


-294- 


aerosol  content  as  input  to  the  code.  The  aerosol  model  used  in  the  code 
is  one  that  gave  good  agreement  with  the  Landsat-2  data,  i.e.,  a Junge 
size  distribution  with  v = 4.0  and  refractive  index  n = 1.5.  Thus,  the 
normalized  radiance  (I)  is  given  by 

I (0,63.3° ,0,N  ) 

1 = I x — — 

Ic(e*eo*+*Nm) 

when  I is  the  measured  radiance 

m 

I is  the  calculated  radiance 
c 

N is  the  measured  aerosol  content 
m 

The  calculations  are  performed  at  0.67  ym  for  the  N0AA-5  Scanning 
Radiometer  (SR)  and  at  0.64  ^ for  the  GOES  Visible  Infrared  Spin  Scan 
Radiometer  (VISSR). 

The  satellite  data  were  obtained  in  digital  form  from  the  National 
Environmental  Satellite  Service  (NESS)  of  NOAA.  The  SR  data  were  available 
only  in  the  mapped  format  (with  20  km  resolution)  for  this  investigation, 
and  did  not  provide  all  the  resolution  elements  actually  measured  during 
the  satellite  overpasses.  Another  shortcoming  of  the  SR  data  is  the  fact 
that  the  SR  output  is  subject  to  a non-random  noise  which  cannot  be  read- 
ily eliminated.  In  an  attempt  to  minimize  these  effects,  each  SR  radiance 
reported  here  is  the  mean  of  a 5 x 5 block  of  pixels  centered  on  the  cal- 
culated ship  location.  The  VISSR  radiances  (with  6 km  resolution)  re- 
ported here  are  also  the  means  of  5 x 5 blocks  of  pixels.  The  radiances 
are  given  in  units  of  mw/cm  /ym/sr,  and  are  based  on  calibration  data 
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obtained  from  NESS,  converting  digital  counts  to  foot-lamberts,  in  the  case 
of  the  SR,  and  directly  into  radiance  units  in  the  case  of  the  VISSR. 

5.  Results 

The  weather  during  the  cruise  was  generally  quite  good,  enabling  coin- 
cident sun  photometer  and  satellite  measurements  to  be  obtained  on  several 
occasions.  The  GOES  coverage  of  the  Atlantic  Ocean  is  good  only  to  about 
25°  W longitude,  which  was  reached  by  the  U.S.N.S.  Hayes  on  May  24,  1977. 

In  these  first  ten  days,  sets  of  data  were  obtained  on  six  occasions  for 
the  GOES.  During  the  same  time  period,  six  sets  of  data  were  also  ob- 
tained for  NOAA-5  overpasses.  However,  after  May  24,  1977,  ten  more  days 
were  spent  at  sea,  but  only  three  more  sets  of  data  were  obtained  for  the 
NOAA-5  overpasses. 

The  GOES  VISSR  data,  plotted  in  Figure  C-Vll-2  show  an  excellent 
linear  relationship,  as  anticipated  from  the  Landsat  results,  and  shows 
that  normalization  procedures  with  the  Dave  code  are  satisfactory. 

The  NOAA-5  data  were  not  expected  to  be  very  useful  due  to  their  for- 
mat and  noise  problems,  as  described  above.  However,  in  spite  of  these 
shortcomings,  the  relationship  shown  for  the  SR  radiances  in  Figure  C-V 1 1-3 
is  remarkably  good.  A linear  relationship  can  probably  be  inferred.  The 
x-points  show  an  enhanced  radiance  due  to  sun-glitter,  and  demonstrate  that 
observations  should  be  made  away  from  the  sun,  except  close  to  the  nadir 
as  illustrated  by  the  o-point. 


In  comparing  these  results  with  those  of  Landsat  2,  taking  into  ac- 
count the  wavelength  differences,  it  is  found  that  for  the  SR  data,  the 
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radiance  value  for  N = 0 is  as  expected  (this  value  is  independent  of  the 
aerosol  properties,  and  represents  a pure  molecular  atmosphere),  but  that 
the  other  radiances  are  lower  than  expected.  This  can  be  due  to  the  aero- 
sol properties  being  different  from  those  of  the  Landsat  San  Diego  data, 
or  can  be  due  to  uncertainties  in  the  radiometric  calibrations  in  each 
satellite.  However,  in  the  Landsat  Study,  data  obtained  at  Adrigole, 
Ireland,  for  Atlantic  Ocean  aerosols  showed  good  agreement  with  the  San 
Diego  data.  The  same  study  showed  that  differences  also  existed  between 
the  Landsat  1 and  Landsat  2 results  at  San  Diego,  and  it  was  concluded  that 
they  were  due  to  differences  in  the  radiometric  calibrations  of  the  two 
satellites.  It  is  believed  that  similar  calibration  problems  are  respon- 
sible for  the  SR  and  Landsat  differences.  Indeed,  in  examining  the  VISSR 
results  in  Figure  C-VII-2,  it  is  found  that  both  the  intercept  and  slope 
of  the  line  are  significantly  different  from  those  predicted  from  the 
Landsat  data,  suggesting  again  that  the  reason  is  due  to  the  radiometric 
cal ibrations. 

6.  Conclusions 

A linear  relationship  between  the  upwelling  visible  radiance,  as  ob- 
served by  the  Landsat  MSS,  and  the  atmospheric  aerosol  content,  has  also 
been  found  for  the  GOES-VISSR  and  the  NOAA-5  SR.  The  relationships  are 
slightly  different  for  each  satellite.  These  differences  are  attributed 
to  differences  in  the  radiometric  calibrations  of  the  satellites,  and 
points  to  the  necessity  of  precise  radiometric  calibrations  of  satellite 
radiometers  if  they  are  to  be  used  in  the  future  for  aerosol  measuremen  ;. 
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Without  precise  calibration  each  satellite  would  have  to  be  empirically 
calibrated  with  lengthy  periods  of  ground  truth  measurements. 
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Aerosol  Content 


Figure  C-VII-1  Landsat  2 Ocean  Radiances  Versus  Aerosol  Content.  The 
Radiances  are  for  Nadir  Viewing  Normalized  to  a Sun 
Zenith  Angle  of  63°. 
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Radiance  (mw/crn /pm/sr) 


Figure  C-VII-2  GOES  Ocean  Radiances  Versus  Aerosol  Content.  The  Radi- 
ances are  Normalized  to  Nadir  Viewing  with  a Sun  Zenith 
Angle  of  63°. 
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Radiance  (nw/cm  /Min/sr) 


Figure  C-VII-3  NOAA-5  Ocean  Radiances  Versus  Aerosol  Content.  The 
Radiances  are  Normalized  to  Nadir  Viewing  with  a Sun 
Zenith  Angle  of  63°. 
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D.  DATA  ARCHIVE 


All  of  the  data  which  was  obtained  on  a more  or  less  continuous  basis 
throughout  the  cruise  and  made  available  at  the  time  of  printing  of  this 
report  has  been  digitized  and  put  into  a form  compatible  with  the  data  ar- 
chive system  being  described  in  this  section  of  the  data  report.  The  pur- 
pose of  this  effort  is  to  allow  for  easy  access  to  any  combination  of  a 
large  number  of  different  variables  which  make  up  this  user  oriented  data 
retrieval  system.  This  system  was  designed  to  be  operated  in  conjunction 
with  the  Digital  Equipment  Corporation,  DEC-10  computer  at  NRL  and  thus 
utilizes  its  standard  software.  The  system  was  created  with  the  view  that 
it  will  be  used  mainly  by  the  scientist/programmer  who  has  the  ultimate 
goal  of  using  real  field  data  to  answer  real  scientific  questions  and  is 
not  interested  in  the  finer  points  of  computer  programming  by  itself. 

Six  major  uses  of  this  system  are  envisioned. 

1)  An  ability  to  obtain  a visual  overview  of  any  particular  variable 
throughout  a period  of  time.  Such  graphic  presentation  of  measured  points 
plotted  as  functions  of  time  provide  a quick  indication  of  trends,  data 
scatter,  data  ranges  as  well  as  areas  of  instrument  malfunction. 

2)  An  ability  to  obtain  scatter  plots  between  any  set  of  2 variables 
measured  simultaneously  to  look  for  correlations  and  significant  areas  of 
deviations  from  correlation. 

3)  An  ability  to  obtain  data  files  of  any  combination  of  variables 
to  be  used  for  processing  of  the  data  in  custom  programs,  and  to  make  it 
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available  to  other  users.  This  ability  of  combining  and  sorting  should 
eliminate  the  need  of  the  data  user  having  to  sort  through  and  work  with 
various  forms  of  data  presented  in  various  formats  and  found  in  various 
diverse  publications  and/or  data  sheets. 

4)  The  system  has  the  ability  to  take  a snap  shot  of  any  set  of  data 
at  a desired  period  of  time  and  it  should  be  helpful  if  this  system  would 
plot  in  a graphic  manner  these  various  parameters  in  the  appropriate  form. 

5)  The  system  is  able  to  make  calculations  from  stored  data  measure- 
ments and  is  able  to  operate  on  these  calculated  data  in  the  same  manner 
as  on  the  original  data  storage  files. 

6)  The  system  is  able  to  perform  statistical  analysis  on  its  multi- 
variate data  base  in  order  to  investigate  correlations  and  other  statisti- 
cal measures  of  relationships  between  the  various  types  of  data. 

The  DEC-10  system  utilizes  disc  resident  files  for  data  storage  and 
the  manipulation  of  these  disc  files  is  a simple  task  for  the  user.  All 
files  on  the  DEC-10  system  consist  of  a name  followed  by  a period  and  an- 
other group  of  up  to  three  letters  known  as  an  extension.  Thus  throughout 
this  report  we  describe  the  location  of  a particular  piece  of  data  as  be- 
ing in  a particular  file.  Examples  of  files  names  used  in  the  DEC-10  sys- 
tem are: 


ASC.DAT 
WHIT.DAT 
CAL. DAT 
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Furthermore  each  data  storage  file  used  in  this  retrieval  system  has  the 


format  shown  in  Table  D-l. 

TABLE  D-l 

(Standard  Format  Data  Storage  Files) 


(1 

1. 

2. 

3. 

• • • • H) 

1 

vi,i 

Vl,2 

Vl,3 

....  V^m 

3 

V31 

V32 

V33 

• • • • V 

*20 

1. 

2. 

3. 

0000 

24 

V24,l 

V24,2 

V24,3 

• • • • V 

Here  D represents  a date  code  such  as  051577  which  stands  for  15  May  1977. 
For  the  EOMET  77  cruise  the  relation  51577  <_  D <_  60777  holds.  The  number 
of  variables  held  in  a particular  data  file  is  represented  in  Table  D-l 
by  the  letter  m.  After  each  0 there  are  m numbers  representing  the  m col- 
umns of  data.  Hours  within  the  day  are  represented  by  the  T's  which  are 
the  GMT  hours  for  the  data  code  immediately  preceeding  it.  Each  date  may 
have  from  0 to  24  time  values  depending  on  how  often  the  variables  were 
measured  throughout  the  day.  Data  stored  in  the  matrix  are  represented  by 
Vi  m where  m 1s  the  column  or  type  parameter  of  the  variable  and  i is  the 
time  parameter.  Thus  all  Vs  with  the  same  i are  considered  simultaneous 
measurements.  This  representation  makes  it  easy  to  put  data  into  the  ar- 
ray and  to  search  for  a particular  piece  of  information  once  it  is  in  the 
array. 

A second  format  (Format  1)  is  used  by  this  system  where  each  line  of 
parameters  is  preceeded  by  a Julian  day  and  a decimal  fraction  of  that  day. 
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In  this  way  time  Is  a monitonically  Increasing  number  which  Is  easy  to 
plot  or  to  compare  with  other  numbers.  Table  D-2  Is  a representation  of 
this  form. 


TABLE  D-2 


Format  1 Type  File 


T1 

V11 

V12 

V13 

V14  * * * Vlm 

T2 

V21 

CM 

CM 

> 

V23 

V24  * * * V2m 

Tn  Vnl  Vn2  Vn3  Vn4  * * * Vnm 

where  the  Tn  are  the  Julian  time  numbers  for  which  there  is  at  least  some 
data  available  in  the  array  and  the  Vnm  are  as  before. 

Three  types  of  files  are  used  in  the  system.  They  are  storage' files, 
working  files,  and  output  files.  The  system  is  designed  so  that  all  of 
the  data  is  stored  in  standard  format  storage  files.  These  storage  files 
are  accessed  mainly  by  two  interactive  programs.  These  two  programs  known 
as  L00K9.BAS  and  FIND. BAS  are  the  means  by  which  the  user  tells  the  com- 
puter which  data  he  is  interested  in  and  what  he  wants  done  with  it. 

Figure  D-l  shows  a block  diagram  of  the  system  where  operational  pro- 
grams are  represented  by  blocks  and  files  which  carry  the  information  be- 
tween the  block  operators  are  represented  by  lines.  The  solid  lines  are 
standard  format  files,  the  dashed  lines  are  format  1 type  files  and  the 
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(lotted  lines  represent  files  which  are  of  various  types  and  used  for  out- 
put and/or  graphing.  Programs  containing  the  extensions  ".BAS"  are  oper- 
ated in  the  BASIC  mode  while  extensions  with  " .FTN"  are  Fortran  programs. 
Blocks  where  the  program  has  no  extension  represent  DEC-10  system  monitor 
commands. 

The  heart  of  the  system  is  a program  called  LQQK9.BAS  which  is  an  in- 
teractive program  written  in  BASIC  which  asks  the  user  various  questions 
about  the  working  file  which  the  program  is  making  for  the  user.  The  final 
product  of  L00K9.BAS  is  a format  1 type  file  which  is  itself  then  used  as 
an  input  file  for  all  of  the  different  programs  following  L00K9.BAS  in 
Figure  D-l.  The  contents  of  this  file  depend  on  the  user's  desires.  In 
essence  this  program  interfaces  between  the  user  and  the  data  bank.  It 
combines  in  one  format  1 type  file  up  to  8 data  columns  of  any  combina- 
tions of  the  simultaneously  measured  variables.  The  inputs  to  L00K9.BAS 
are  essentially  data  storage  files  but  a provision  has  been  added  to  use 
an  "EXTRA"  file.  This  is  the  feature  which  allows  the  system  to  test 
models,  to  manipulate  data  and  to  perform  calculations  on  the  measured 
data  and  to  have  this  data  available  for  further  use  on  the  system  such  as 
plotting,  etc. 

The  output  file  from  L00K9.BAS  has  four  main  uses: 

1)  Graphical  Presentation 

2)  Numerical  Output 

3)  Statistical  Analysis 

4)  Data  manipulation 
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Examples  of  these  various  types  of  programs  which  exist  in  the  system  are 
shown  in  Figure  D-l.  The  methods  of  obtaining  analytical  and  numerical 
output  from  the  data  bank  are  listed  in  Table  D-3.  Requests  for  this  data 
should  be  sent  to: 

Stuart  Gathman 

Naval  Research  Laboratory 

Code  8327 

Washington,  D.C.  20375 

Program  FIND. BAS  is  a snapshot  routine  which  is  used  to  obtain  dis- 
plays and  listing  of  multivariate  data  such  as  aerosol  size  spectra  and 
atmospheric  sounding  profiles.  The  program  is  interactive  and  asks  the 
user  what  data  he  wants  and  what  specific  date  and  time  he  would  like 
plotted  or  listed.  This  program  has  access  to  all  of  the  Data  Storage 
files  as  well  as  to  "EXTRA"  type  files  which  can  be  constructed  by  means 
of  the  L00K9.BAS  routine. 

Table  D-4  is  a catalogue  of  the  various  types  of  data  that  are  in  the 
archive.  The  physical  quantity  is  underlined  in  the  left  hand  column  and 
listed  directly  below  are  the  various  measurements  of  the  physical  quanti- 
ties. Each  is  identified  in  some  way  by  the  initials  of  the  performing 
organization.  Details  on  the  specifics  of  these  measurements,  i.e.,  the 
location  of  the  sensors,  type  of  instrument,  is  discussed  in  Section  B-II 
of  this  report.  In  the  right  hand  column  are  the  data  storage  files  names 
and  in  parenthesis  are  the  specific  data  columns  in  which  the  data  is 
kept.  Program  FIND. BAS  requires  just  the  name  of  the  file  as  it  will  work 
on  all  of  the  data  at  a particular  time  in  a particular  file  whereas  pro- 
gram L00K9.BAS  requires  both  a file  name  and  the  data  column  number. 
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TABLE  D-3 


Methods  of  Obtaining  Numerical  Data  From  the  Archive 


1)  Files  may  be  typed  by  a teletype  or  other  type  of  hard  copy  device. 

2)  Files  may  be  printed  on  the  system  line  printer  and  mailed  to  user. 

3)  Graphical  reproductions  are  available  in  Polaroid  reproduction  of 
plot-10  programs.  This  includes  a time  plot  of  any  variable  in  the  ar- 
chive or  a scatter  plot  of  any  two  variables  in  the  archive. 

4)  Facilities  are  also  available  to  obtain  hard  copies  of  all  graphics 
utilizing  the  Tektronix  4631  hard  copy  unit. 

5)  Mag  tapes  may  be  made  with  a variety  of  modes  using  the  program 
CHANGE  on  the  DEC-10.  These  will  be  9 track  800  bpi  tapes  in  the  follow- 
ing modes: 


ASCII 

6.  E.  ASCII 

H.  P.  ASCII 

HONBCD 

BOD 

BOL 

G.  E.  BOD 
FIXSIX 
EBCDIC 
SIXBIT 

6)  ASCII  paper  tapes  can  be  punched  on  the  system  paper  tape  punch. 

7)  Punched  card  decks  can  be  made  on  the  NRL  Advanced  Scientific 
Computer. 
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TABLE  D-4 


DATA  ARCHIVE 


Sea 

Surface  Temperature 

1) 

NRL  Keel  Level  Probe 

2) 

Cal  span  Fish 

3) 

NPS  Fish 

4) 

NRL  PRT-5( IR) 

5) 

NRL  PRT-IO(IR) 

6) 

Average  SST 

Air  Temperature 

1) 

NPS  Bow 

2) 

Cal  span  Bow  -9m 

3) 

NRL  8330  Rear  Deck 

4) 

NRL  8327  Tower 

5) 

Cal  span  Tower  - 15  m 

6) 

NRL  8323  Tower 

7) 

NPS  Mast 

8) 

Cal  span  Mast  - 27  m 

9) 

NRL  Kytoon  System  Profiles 

101 

Average  Air  Temp. 

Temperature  Fluctuations 

1) 

2) 

NPS  C2  . Bow 
?T 

NPS  Cr  _ Mast 

r 

3) 

NPS  e Bow 

4) 

NPS  e Mast 

Humi  di  ty 

1) 

NPS  Bow 

2) 

NRL  8330  Rear  Deck 

3) 

NRL  8323  Tower  (Dewpoint) 

4) 

Calspan  Tower  - 15  m 

5) 

Cal  span  Sling,  Tower 

6) 

NPS  Mast 

7) 

Calspan  Mast  - 27  m 

8) 

Kytoon  System  Profiles 

9) 

Best  Relative  Humidity  (%) 

File  Name 


THERM. DAT(2) 

CALTMP.DAT(5) 

NAVPGS.DAT(2) 

FILE.DATI4) 

PRT10.DAT(2) 

AVE.DAT(2) 


NAVPGS.DAT(3) 
CALTMP.DAT(4) 
THERM.DAT(5) 
THERM.DAT(3) 
CALTMP.DAT! 3) 

FILE -DAT( 2 ) 
NAVPGS.DAT (4) 
CALTMP.DAT! 2) 
DATA.DAT(2,3,4,5) 
AVE.DATI3) 


NAVPGS.DAT! 7) 

NAVPGS.DAT(8) 

NAVPGS.DAT! 9) 
NAVPGS.DAT(IO) 


NAVPGS.DAT! 5) 

THERM.DAT(4,5) 

FILE.DAT(3) 

CALTMP.DAT(7) 

CALVIS.DAT(6,7) 

NAVPGS.DAT(6) 

CALTMP.DAT! 6) 

DATA2.DAT! 2, 3, 4, 5) 

AVE.DAT(4) 
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File  Name 


r 


Relative  Wind  Speed 

1) 

Cal  span  Mast 

CALVIS.DAT(3) 

2) 

NRL  Mast 

ASC.DAT(2) 

Relative  Wind  Direction 

1) 

Calspan  Mast 

CALVIS.DAK2) 

2) 

NRL  Mast 

ASC.DAT(3) 

Real 

Wind  Speed  and 

Direction  Calculations 

1) 

Speed 

WDREAL.DATI3) 

2) 

Direction 

WDREAL.DAT(2) 

Ship  Speed  and  Heading 

1) 

Heading,  magnetic,  Calspan 

CALTMP.DAT(8) 

2) 

Course  Made  Good 

NAV.DATI6) 

3) 

Ship  Speed 

NAV . DAT ( 7 ) 

* * * * * AEROSOLS  * * * * * 


Calspan  Royco 

a) 

Royco  > .3/2.8  liter 

0UT.AER(3) 

b) 

Royco  > .6/2.8  liter 

OUT . AER ( 4 ) 

c) 

Royco  > 1.2/2. 8 liter 

0UT.AERI5) 

d) 

Royco  > 3. 0/2. 8 liter 

OUT .AER (6) 

e) 

Royco  > 5. 0/2. 8 liter 

0UT.AER(7) 

NRL  Royco  #/cc 

a) 

0.45  urn  < d < 0.60  ym 

R0Y.DAT{2) 

b) 

0.60  ym  < d < 1.50  ym 

ROY . DAT ( 3 ) 

c) 

1.50  ym  < d < 2.5  ym 

ROY . DAT ( 4 ) 

d) 

2.50  ym  < d < 4.0  ym 

R0Y.DAT(5) 

e) 

d 4.0  T ym 

R0Y.DAT(6) 
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File  Name 


Thermo  Systems  - E.A.A. 


a) 

login  N (r  > .0032  u) 

WHIT2.DAT{2) 

b) 

log{«  N (r  > .0056  u) 

WHIT2.DAT(3) 

c) 

log})}  N (r  > .01  v) 

WHIT2.DAT(4) 

d) 

log})}  N (r  > .0178  y) 

WHIT2.DAT(5) 

e) 

log})}  N (r  > .0316  y) 

WHIT2.DAK6) 

f) 

log})}  N (r  > .0562  y) 

WHIT2.DAK7) 

g) 

log})}  N (r  > .1  w) 
log}”  N (r  > .178  y) 
log}}}  N (r  > .360  u) 
log}”  N (r  > .560  y) 

WHIT2.DAT(8) 

h) 

WHIT2.DAT (9) 

i) 

WHIT2.DAT(10) 

j) 

WHIT2 .DAT (11) 

Visibility 

1) 

MRI  Nephelometer  (NRL) 

ASC.DAT(5) 

2) 

MRI  Nephelometer  (Cal span) 

CALV1S.DAT(5) 

3) 

EG&G  Foward  Scattering  (Cal span) 

CALVIS.DAT(4) 

4) 

MRI  Fog  Visibility  (NRL  8323) 

ASC.DAT(6) 

5) 

Observer  Estimates 

CLOUD. DAT(4) 

Inversion  Height 

1) 

SRI  Lidar 

SR I .DAT ( 2 ) 

Radon 

1) 

A.R.C.  (NRL  8326) 

RAD0N2.DAT(2) 

CN 

1) 

Calspan  - Gardner  Counter,  Aiken  Nuclei 

0UT.AER(2) 

2) 

NRL  Hoppel  - Poliak  Counter 

CN.DAT(2) 

3) 

E 1 Data  NRL  Hoppel  log^(#/cc) 

ASC.DAT(7) 

CCN 

(NRL  Code  8326) 

1) 

#/cc  at  2°C 

CCN.DAT(2) 

2) 

#/cc  at  3°C 

CCN.DAT(3) 

3) 

#/cc  at  4.3°C 

CCN.DAT(4) 

4) 

#/ cc  at  5°C 

CCN.DAT(5) 

5) 

#/cc  % Involatile  at  4.3°C 

CCN.DAT(6) 

Scattering  Coefficient 

1) 

MRI  Nephelometer  Calspan 

0UT.AER(8) 

2) 

MRI  Nephelometer  NRL 

ASC.DAT(4) 
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File  Name 


r 


Meteorological  Observations 


1) 

Sea  State 

CLOUD. DAT{2) 

2) 

% White  Caps 

CL0UD.DAT(3) 

3) 

% Cloud  Cover 

CLOUD. DAT(5) 

4) 

Cloud  Type  - Low  Cloud 

CLOUD. DAT (6) 

Middle  Cloud 

CL0UD.DAT(7) 

High  Cloud 

CLOUD. DAT(8) 

5) 

PRT-10 

i)  Sea  Surface  Temp. 

PRT10.DAT(2) 

i i ) Sky  Temp. 

PRT10.DAT(3) 

iii)  -10  deg. 

PRT10.DAT(4) 

iv)  +10  deg. 

PRT10.DAT( 5) 

v)  +30  deg. 

PRT10.DAT(6) 

vi ) +60  deg. 

PRT10.DAT( 7) 

vii)  Cloud  Base 

PRT10.DAT(8) 

Stability 

1) 

U*NPS 

NAVPGS.DAT(ll) 

2) 

R.NPS 

J 

NRL  Lapse  Rate  Kytoon  System 

NAVPGS.DAT (12) 

3) 

TKITE.DAT(3) 

Volz 

. Sun  Photometer  Readings 

1) 

Air  Mass 

V0LTZ.DAT(2) 

2) 

0.88  y VSP  Reading 

V0LTZ.DAT(3) 

3) 

0.50  y VSP  Reading 

V0LTZ.DAT(4) 

4) 

0.94  y VSP  Reading 

V0LTZ.DAT (5) 

5) 

Turbidity  (0.88  y) 

TURBID. DAT(2) 

6) 

Turbidity  (0.50  y) 

TURBID. DAT(3) 

7) 

Turbidity  (0.94  yj 

TURBID. DAT( 4) 

8323 

Tower  Positions  (Up  or  Down) 

UPD0WN.DAT(2) 

Ship  Position 

1) 

Latitude  (deg) 

NAV.DAT(2) 

2) 

Latitude  [min.  (decimal)] 

NAV.DAT(3) 

3) 

Longitude  (deg) 

NAV.DAT(4) 

4) 

Longitude  [min.  (decimal)] 

NAV.DAT(5) 
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4 Lakeside  Office  Park 
Wakefield,  Mass  01880 
617-246-2838 

Clamons,  Dean 
Code  8003 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-2384 

Cosden,  Thomas 
Code  5568 

Naval  Research  Laboratory 
Washington,  D.C.  20175 
202-767-3224 

Fairall,  Christopher  (Dr.) 
Naval  Postgraduate  School 
Monterey,  Calif  93940 
408-646-2219 


Fitzgerald,  James  (Dr.) 

Code  8326 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-2362 

Fraser,  Alistair  (Dr.) 

Penn  State  Univ. 

507  Deike  Blvd 

University  Park,  Penna  16802 
814-865-0478 

Gerber,  Hermann  (Dr.) 

Code  8322 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-2780 

Hayes,  John  N.  (Dr.) 

Code  8326 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-3589 

Hoppel,  William  A.  (Dr.) 

Code  8326 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-2362 

Jeck,  Richard  K.  (Dr.) 

Code  8323 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-2437 

Julian,  Ben  G. 

Code  8327 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-2022 

Katz,  Ulrich  (Dr.) 

Calspan  Corp. 

P.  0.  Box  235 
Buffalo,  N.  Y.  14221 
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Kidwell,  Raymond  Stilling,  Robert 

Code  8323  Code  8322 

Naval  Research  Laboratory  Naval  Research  Laboratory 

Washington,  D.C.  20375  Washington,  D.C.  20375 

202-767-2437  202-767-3317 


Larson,  Reginald  E. 

Code  8326 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-3589 

Leonard,  Donald  (Dr.) 

Computer  Genetics  Corp. 

4 Lakeside  Office  Park 
Wakefield,  Mass  01880 
617-246-2838 

Mack,  Eugene 
Cal  span  Corp 
Environmental  Sup.  Dept 
P.  0.  Box  235  (or) 

4455  Genessee  Street 
Buffalo,  N.  Y.  14221 
716-632-7500,  x643 

May,  Lenwood 

Naval  Postgraduate  School 
Monterey,  Calif  93940 
408-646-2219 

Nielson,  Norman 

Atmospheric  Sciences  Laboratory 
Stanford  Research  Institute 
Menlo  Park,  Calif  94025 
415-326-6200  x2838 

Rosenwasser,  Bruce 
Code  8320 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-2022 

Schefer,  Murray 
NavAir-Code  370 


Trusty,  Gary  (Dr.) 

Code  5568 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-3224 

Twomey,  Sean  (Dr.) 

Institute  of  Atmospheric  Physics 
University  of  Ariz. 

Tucson,  Ariz  85721 
602-626-1327 

Uthe,  Edward  (Dr.) 

Atmospheric  Sciences  Laboratory 
Stanford  Research  Institute 
Menlo  Park,  Calif  94025 
415-326-6200  x 4667 

Wojciechowski , T.  A. 

Code  8326 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-3589 

Zuccaro,  Anthony 
Code  8004 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
202-767-2731 


Spalding,  Glenn 
NavMat 


-316- 


